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Final Report on NAS8-37586 


Under this contract UAH was involved in a collaborative program with MSFC to 
carry out design and development activities for the Ultraviolet Imager to be flown on the 
Polar Spacecraft of the International Solar Terrestrial Physics {ISTP) Mission. In 
addition several other related tasks were carried out as per statement of work. These 
are described below. 


The following tasks were performed: 

1 . Design and Fabrication of Prototype/Engineering Model of the UVI 
Imager 

Following completion of the conceptual design and table top testing at UAH, 
modifications were identified in accordance with the test results and, a subcontract was 
awarded to Perkin Elmer Corporation to optimize the UAH design and to develop and 
fabricate an engineering model of the optical bench. 

The program was initiated with a statement of work to Perkin Elmer on 
2/15/89 This was followed by a meeting on March 31, 1989. Work commenced with 
UAH issuing a letter of intent on 4/17/89. However, to provide the funding required to 
support the overall effort at Perkin Elmer, additional funding had to be provided through 
a second contract supporting ISTP activities at UAH, namely NAS8-36955, an umbrella 
contract with MSFC that was used to provide the balance in the short turn-around needed 
to get the Perkin Elmer effort going on schedule. The funding was provided under 
Delivery Order Number 59. 

The work was divided into 2 phases. During phase 1 the following tasks were 
accomplished. 


The system and subsystem level requirements were defined 


The UAH optical design was optimized 

A sensitivity analysis of the final optical design was performed 

A straylight analysis was performed and the conceptual layout of the 
baffles was done with a rework done under D.O. 59 

The conceptual layout of the optical bench was generated, with final fine- 
tuning done under D.O. 59 

The design and layout of the optical elements was finalized and alignment 
tolerances and mounting interfaces defined 

During Phase II an engineering model of the optical bench was fabricated and 
detailed engineering drawings developed. Copies of the drawings which were 
preliminary are attached. 

The unit incorporated appropriate interfaces for the detector assemblies, filter 
wheel and other system interfaces. Critical components were lightweighted. A coarse 
NASTRAN model of the optical bench was produced and the delta temperature limits were 
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analyzed. A weight analysis was performed. System alignment tests were performed and 
an alignment procedure generated. 

The deliverable of a mid-term, final report and Engineering unit Ultraviolet 
Imager with associated engineering drawings were delivered to MSFC on June 11, 1990. 


Note: During the design development UAH identified the need for a second s y s J®™ °| 

reflective filter components in the optical train which were needed to provide additional 
filtering to meet the UVI science spectral purity requirements. To incorporate this 
filter requirement Perkin Elmer was requested to insert a 45° reflective surface at the 
entrance aperture to the system and to fold the optical beam through 90 . This effort 
was carried out under subcontract NAS8-36955 - Delivery Order Number 59 and will 
be reported on separately. This was the main task carried out on D.O. 59. 


Attached as Appendix A to this report are the Mid-Term and Final Reports 
received from Perkin Elmer (now Hughes Danbury). Note: Hughes Danbupr did not 
separate out the two tasks in their reports and their final design reflects the above- 
mentioned changes. The reports give details on the list of tasks discussed above including 
the performance evaluation of the design. 

It is noteworthy that the UAH design which was fine-tuned by Perkin Elmer 
constitutes a technological breakthrough in far ultraviolet imaging in that it represented 
the first high speed (f/# = 2.8) imaging camera with high spatial resolution developed 
to date for this wavelength regime. The sensitivity improvement of the instrument over 
existing systems exceeded a factor of 10, placing UAH at the forefront of this technology. 


2 . Preliminary Design Review 

UAH supported the UVI PDR at MSFC on October 17, 1989 in the following areas: 

Presentation of the optical system (not included in the attached report on 
PDR) 

Presentation of the data analysis software 
Report on the status of the filter design 

Note: The data analysis software flowchart shown in the PDR report was 
generated by Science and Engineering Associates. However, UAH provided the science 
algorithms needed to evaluate the dependence of intensity ratios on the characteristic 
energy of precipitating particles. UAH presented the UVI scientific goals and described 
the model development that has been used to support the UVI requirements definition. 

Relevant sections of the UAH contributions to PDR can be found in the copies of 
the PDR presentation included as Appendix B and in the Perkin Elmer reports (Appendix 
A). Because Perkin Elmer could not attend PDR, the optics review was conducted with 
the GSFC review team optics specialist in a splinter session telecon with Mr. Andreas 
Nonnenmacher of Perkin Elmer. Dr. D. G. Torr convened this session. 


3 . Vacuum Ultraviolet Filter Design 

Narrowband filters with the spectral performance specified for the UVI had 
never been fabricated before. In order to meet the spectral requirements of the filters, 
a novel design strategy was developed. First, it was discovered that suitable materials 
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for use in the VUV had never been identified, and second, the effects of absorption which 
occurs in the VUV had not yet been incorporated into thin film design theory. Third, 
accurate measurement techniques for determining the transmission and reflectance of 
thin films (needed for the determination of optical constants) had not yet been developed. 

All these shortcomings were addressed as follows: 

A technique was developed for the experimental determination of the 
reflectance and transmission of thin films and substrates in the VUV 

A mathematical algorithm was developed for the retrieval of optical 
constants of thin films and substrates from measurements of their reflectance and 
transmittance. 


Thin film multilayer theory was modified to include the effects of 
absorption. This essentially involved starting from scratch with Maxwell's equations. 

A detailed report of the above work was published in two papers entitled "Vacuum 
Ultraviolet Thin Films 1 : Optical Constants of BaF 2 , CaF 2 , LaF 2 , MgF 2 , AL 2 0 3 , HF0 2 and 
Si0 2 Thin Films" and "Vacuum Ultraviolet Thin Films 2: Vacuum Ultraviolet All- 

dielectric Narrowband Filters" attached as Appendix C. 


Reprints of both papers are attached. However, it must be noted that the work 
reported in these two papers covers the overall development of the program 
during the course of the ISTP program under contracts NAG8-086, NAG8-639, NAS8- 
37576, and D.O. 59. The final reports on the above four contracts describe the early 
phase of the work including establishing a VUV coating facility at UAH, debugging it, and 
developing the state-of-the -art filter technology in existence at UAH today. e 
fabrication of the engineering model filters being done under this contract. 


This work will have very broad applications in NASA and will open up the VUV to 
high resolution interferometry, photometry and high powered VUV lasers with the 
future development of Fabry-Perot etalons, beam splitters and high reflectance 
mirrors. 


At the time of PDR two prototype filters had been developed and their spectral 
characteristics are given in the PDR report attached as Appendix B. Also shown in the 
PDR report are the characteristics of commercially available VUV filters. A comparison 
of these with the UVI filters demonstrates that we have achieved an order of magnitude 
improvement in performance. The filter at 135.6 nm is a narrowband filter, and the 
filter labeled LBH l0N g is a specially designed broadband filter with a rectangular 
bandpass which meets the specific requirements of the ISTP science. 

Fabrication of the flight filters was carried out under the prime contract NAS8- 
38145 and is described separately in the reports on that contract. 


4 . Auroral Energy Deposition Code 

UAH was also responsible for the development of a computer code to calculate the 
auroral energy deposition rate from the images to be taken by the UVI instrument. This 
code is needed for two purposes: 
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1 . To determine what wavelengths and filter characteristics are needed for 
the UVI. 

2. To provide the auroral source function to a global ionospheric model 
which will be used by the ISTP community in the interpretation of the measurements 
taken. 


This task was a major undertaking and under NAS8-37586 the code was 
developed to a point where the filter design requirements could be finalized. Results 
achieved to date were published in the Journal of Geophysical Research under the title, 
"Auroral Modeling of the 3371 A Emission Rate: Dependence on Characteristic Electron 
Energy." In addition the application of this code to interpretation of the UVI images was 
published in the same journal under the title, "The Dependence of Model 01 1356 A and 
N 2 LBH Auroral Emissions on the Neutral Atmosphere". Other relevant associated papers 
partially supported under this contract are also attached as Appendix D. 

5 . Model of LBH Vehicle Glow 

This component of the work was carried out to study the effects of optical 
contamination due to interactions of the spacecraft with the natural environment. 
Previous observations of Vacuum Ultraviolet glow had been observed on Spacelab 1 and 
on the DoD satellite S3-4. A preliminary model of a mechanism that could generate the 
glow was developed under NAS8-057 "Assessment of Vehicle Contamination In Spacelab 
1 Spectroscopic Data." The model invoked surface production (and desorption) of 
metastable N 2 (A) via surface recombination of atomic nitrogen. The N 2 (A) is then 
collisionally excited to the a 1 rig state which radiates in the LBH bands. However, no 
adequate mechanism for the production of surface N could be found. In addition the 
mathematical development of the mechanism to account for the altitude dependence of the 
LBH glow was incorrect. Nevertheless, the basic concept of surface production of N 2 (A) 
leading to a^g was a major conceptual step forward. The calculations also utilized the 
concept of a gas build-up on vehicle surfaces which D. Torr had previously developed in 
conjunction with R. Rantanen of Science and Engineering Associates, Inc. This approach 
also proved to be crucial in the development of the understanding of the glow. 

Under the current contract the following new developments occurred: 

1 . A source for surface N was found. It was argued that ambient N 2 
impacting the vehicle surface would acquire vibrational excitation which could provide 
the exothermicity needed for the reaction 

N 2 (X 1 E + g) v > 13 + O -> NO + N 
to proceed, thereby providing an abundant source of surface N. 

2. The surface chemistry was updated to take account of the production and 
loss mechanisms for N 2 (X) v > 13, O, N and N 2 (A). 

3. The effects of radial outflow on the N 2 (A) gas cloud distribution was 
included. 

4 . An analytical formulation of the model was generated in which the scale 
height dependence of both the N 2 (A) and N 2 (a 1 ng) could be easily predicted. Previous 
calculations of Rantanen and Gordon of the re-emitted gas cloud around a 1 m disk were 
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scaled to provide an estimation of the contaminant cloud concentrations. The model 
correctly predicted both N 2 LBH and Vegard Kaplan (VK) emission intensities for the 
glow observed on Spacelab 1 and S3-4 as well as the height dependence observed on the 
S3-4 satellite. A paper reporting these results was submitted for publication to the 
Journal of Geophysical Research, and a copy of this version of the paper is attached as 
Appendix E. However, one reviewer argued that VK emissions were not observed on S3- 
4 and thus questioned the validity of the model. This was the status of the paper at the 
conclusion of the work under this contract. Since that time, the work has been continued 
by Science and Engineering Associates under DoD funding and it appears that the conflict 
with the S3-4 data will be resolved. 

6 . Laboratory Measurement Program of Collision Cross-Sections 

Dr. C. Keffer was supported for the period June 1 through September 30, 1989 
on this contract. During this period the following tasks were accomplished. 

Completed drawings for using the new large vacuum chamber as an EUV 
calibration facility. 

Implemented the EUV calibration facility and conducted tests on a CCD 
detector system. 

The cross-section chamber was assembled which was delivered mid July, 
1989. This task involved numerous minor time consuming activities. For example, the 
rotary platform did not turn freely when installed, because the mating flanges were not 
flat. All parts had to be cleaned before use of the 5 eV O beam. Assembly was completed 
by mid September. 

The system appeared to be working satisfactorily upon pump-down. 

Development of the system software was continued. 

Alignment tests commenced mid September, but it was found that the 
chamber goes out of alignment while under vacuum. Additional hardware was ordered to 
repair the problem. 

This effort represented all the work done on the establishment of the 
cross-section facility at MSFC under this contract. Continued work on this project was 
funded under other contracts. 

7 . Support of ISTP Meetings 

UAH personnel attended and participated in the ISTP design review meetings at 
GSFC and GE, Science Working Group meetings and meetings held with the NASA 
Headquarters Program Office personnel to resolve the issue of merging the UVI and 
Visible Imager teams of Iowa and Johns Hopkins Universities. 
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APPENDIX A 


Research Administration 



The University 
Of Alabama 
In Huntsville 

Huntsville, Alabama 35899 
(205) 895-6000 
Telefax: (205) 895-6677 


MEMORANDUM 


TO: 

FROM: 

SUBJECT: 

DATE: 


Dr. D. G. Torr 
Physics 

w- 

Kathy Niemi 
Contract Assistant 

Final Report for SUB89-117 amd SUB90-064 
August 6, 1990 


I have enclosed a copy of the final report on the above 
referenced subcontracts. If this report satisfactorily meets 
the reporting requirements of these subcontracts and all 
deliverables have been received and accepted^ please sign the 
concurrence line below signifying such and return this 
memorandum to me. 


If you have any questions, I can be reached at the above 
number . Thank-you . 



An Affirmative Acnon/Eauai Opportunity institution 



HUGHES 


Danbury Optical Systems. Inc 



RECE 


i\. uT~ 

i * ± 




CD-JMC-1611 
July 16, 1990 

University of Alabama 
Research Administration 
Huntsville, AL 35899 

Attention: R. McKinley 

Subject: Ultra-Violet Imager Final Report 

Reference: Subcontracts SUB89-117 and SUB90-064 

Dear Ms . McKinley: 

Enclosed is a Final Report for the Ultra-Violet Imager Optical 
System. This document is submitted in accordance with final report 
requirements and completes our contractual obligations under the 
referenced subcontracts . 


Very truly yours, 


HUGHES DANBURY OPTICAL SYSTEMS 
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530E-06 

-2.75E-06 

Housing 

3.1804 

Al 6061-T651 

1J6E-05 

-4X3E-04 

PM Flang to Housing Shim 

0.0500 

AISI 410 

530E-06 

-2.75E-06 

PM Body to Flange 

0.2499 

Al 6061-T651 

1J6E-05 

3.40E-05 


Al ■ ± 0.00034 in. 

(± 0.009 mm) 

Manufacturing and Assembly Errors 


Measurement of Optical Surfaces to Reference Surfaces : ±0X0014 in. (includes both PM and SM) 

Measurement of Spadng between Optics: ±0.0005 In. 

Total Errors: ±0.00052 in. (±0.0132 mm) 


A conservative total Primary Mirror w.r.t Secondary Mirror Despace Error is given by the rss of these 
values: 


£ * ± 0X014 in. (± 0.035 mm) 

Total Allowable: e» « ± 0X02 in. (±0X50 mm) Margin: £m - ± 0X014 In. (±0.035 mm) 
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The 3 element telescope system is flexible enough to allow for 
re-optimization of the focal plane location to compensate for 
fabrication and/or alignment errors. 
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1. Introduction 


This report summarizes the approach taken and analyses performed in 
deriving the optical performance error budget for the International Solar-Terrestrial 
Physics Ultra-Violet Imager (UVI). 

Presently under contract to the University of Alabama /Huntsville for the UVI 
engineering model (Subcontracts SUB89-117 and SUB90-064), Hughes DOS 
responsibilities formally encompass only the UVI telescope, i.e. the three off-axis 
aspheric mirrors and the structure required to maintain alignment between these 
optics, telescope baffling, interfaces to the spacecraft and MSFC hardware, as well as 
the associated systems engineering, and optical (encircled energy and straylight 
performance) and structural dynamics analyses. 

Additional analyses, primarily a detailed thermal analysis of the instrument to 
determine ground-to-orbit and on-orbit environments, are being performed by 
MSFC. As results from these analyses are released, appropriate revisions to the 
error budget will be made. 



2. Ultra-Violet Imager Mission Description 


The Ultra-Violet Imager mission objective is to spatially and temporally map the 
aurora oval in the vacuum ultra-violet region of the spectrum, thereby providing a 
description of the Earth's magnetic field boundaries. This in turn offers one insight 
into the behavior of the energetic particles associated with the solar wind. In 
particluar, total particle energy in-flux, particle acceleration processes, and 
characteristic energies of incoming particles can be determined. A major mission 
objective is the ability to image both the sunlit and nightside aurora, therefore 
requiring an instrument with extremely good straylight rejection and out-of-band 
rejection capabilites. 

UVI is to be mounted on the single-axis despun platform of the POLAR 
spacecraft which will orbit the earth in a highly elliptical polar orbit, providing 
viewing of the entire northern and southern hemispheres. Since observation of the 
entire aurora oval is desired at orbit apogee (approximately 9 earth radii), the UVI 
telescope has been designed with an 8° circular full field-of-view. 



3. UVI Telescope Description 


The Ultra-Violet Imager telescope being fabricated is a three mirror design with 
an effective focal length (EFL) of 123.97 mm and an entrance pupil diameter of 42.98 
mm implying a telescope speed of f/2.88. However, since the aperture of the system 
is 67.28 mm eccentric to the optical axis, the system is, from an optical design point 
of view, effectively an f/0.7 system, making it an extremely fast telescope. 

Figure 3.1 shows a functional plan view of the final optical layout. Depicted 
are optical element bodies (without mounting flange details), mirror body 
diameters, and nominal dimensions between optics showing axial locations and 
how far off-axis each element is located. Note that the filter mirror (refer to Figure 
3.5) is not depicted here. Table 3.1 provides some nominal system parameters. 


Table 3.1 UVI System Parameters 


Telescope EFL: 
Entrance Pupil o: 
System f/ #: 

Full Field of View: 
Field Bias: 

Plate Scale: 
Detector Pixel Size: 


123.968614 mm 

42.976800 mm (off-set 67.2846 mm from optical axis) 

2.884547 

8° (circular) 

- 6 ° 

2.1637 mm /degree 
54 pm x 32 pm (apparent) 


The original candidate three mirror optical system with a -6.25° field bias (Figure 
3.2) was provided to Hughes DOS by MSFC. Using MEXP, a HDOS proprietary 
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ENTRANCE PUP] L PLANE 



optical analysis code, preliminary analyses showed that the system had relatively 
poor on-axis encircled energy (EE) performance (48% in a 27 pm radius) which 
varied greatly across the total FOV from a minimum of 31% (at -4°,0°) to 87% (at 
+4°,0°) in the 27 |im radius. 

The three mirror system was subsequently optimized adjusting the spacings 
between the optics, changing the field bias to -6.0°, and altering the aspheric terms of 
the three mirrors (Figure 3.3), resulting in improved on-axis performance as well as 
minimizing variations across the field. The on-axis EE performance was improved 
to 71% in a 27 pm radius, and the EE variation across the FOV ranged from 71% (at 
0°, 0°) to 86% (at -4°, 0°) in the 27 pm radius. Figure 3.4 show spot diagrams for this 
system. Figure 3.5 show the final instrument configuration with the reflective filter 
mirror in front of the transmissive filters. This change was made because of filter 
fabrication constraints and provides for improved performance in terms of 
straylight rejection. 


Figure 3.2 

UVI - Original Optical Prescription 
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Figure 3.3 

UVI - Optimized Optical Prescription 
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Figure 3.4 

Spot Diagrams for Optimized System 








Figure 3.5 

Optimized Prescription with Filter Mirror 
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4. Development of the Error Budget 


The proper criterion for evaluating optical performance depends on a number of 
different factors: the nature of the source to be imaged, limitations of the optical 
system, the type of detector, and the goal of the application. Figure 4.1 shows the 
relationships between various commonly used image quality criteria. 

Since this telescope is dominated by geometric abberation as shown in the spot 
diagrams (the on-axis rms spot radius is =54 times larger than the diffraction-limited 
performance of an f/2.88 system at 1200 A), error budgeting in terms of wavefront 
error would not provide a practical figure of merit. Because spatial and temporal 
intensity variations in the scene are of primary interest to the instrument scientists, 
encircled energy and rms spot size, which are calculated from the point spread 
function, were chosen to provide a convenient means of assessing telescope 
performance. 

An optical transfer function (OTF) budget (or equivalent) might well be 
considered appropriate for this system given the desire to spatially resolve aurora 
features. OTF, however, can in turn be translated via a Fourier transform into a 
point spread function from which the original image quality criteria were 
calculated. 

Since no specific optical performance requirements were handed down to 
Hughes DOS by MSFC or the University of Alabama, the approach of "best effort" 
was taken giving the constraints of manufacturing and assembly tolerances, 
environmental contributors, and budget and schedule. Traditionally, a degradation 
of 15%-20% in the optical performance of a telescope due to environmental 
influences, and fabrication and assembly errors has been tolerated. With is in mind, 
a target of 60% encircled energy in a 27 pm radius circle, corresponding to the half- 
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some of the different image quality criteria. 


width of the long side of UVI's CCD pixels, was established. This requirement is 
contingent on a reasonable set of manufacturing and assembly tolerances for both 
the optics and the assembly of the telescope, i.e. the use of ultra-high precision (sub- 
micron) alignment techniques would not be used. 



5. UVI Error Budget 


Since cost is a major program consideration and given the speed of the parent 
optics, it was decided that the optics should be diamond-turned using nickel-coated 
aluminum substrates and subsequently post-polished. This provides one with very 
smooth (typically 10-12 A rms) surfaces and allows one to put high-order aspheric 
surfaces on the substrate. Based on input from the optical designers, it was 

determined that the optical surface should be < 0.5 waves rms (X = 6328 A) for each 
optic. This amount of surface figure error would not noticeably affect encircled 
energy or rms spot size (i.e. < 1%) and is comfortably achieved with current 
diamond-turning practices. 

The principle contributor to optical performance degradation for this optical 
system is the large bulk temperature change associated with the transition from 
ground to orbit. Based on MSFC's initial estimate of the on-orbit thermal 
environment, the nominal on-orbit bulk environment was specified to be 0 °C (32 
°F) with allowable orbit-to-orbit bulk temperature variations of ±5.5 °C (±10 °F) 
while the instrument was operating. Bulk temperature changes contribute to 
relative despace and decenter misalignments, as well as changes to the radius of 
curvatures of the three mirrors. Temperature gradients through the optical bench 
in each of the axes were estimated to be ±0.1 °C/cm (±0.5 °F/in.). These gradients 
were assumed to be superimposed over the bulk average, thereby contributing only 
to changes in the mirror radius of curvatures and relative tilts between the optics. 

Because UVI is compromised in terms of on-axis performance so that 
performance variations across the large FOV could be minimized, a simple error 
sensitivity analysis could not be performed. As a result, once error allocations were 
determined, each misalignment case had to be run through MEXP separately to 
determine its impact on performance. 



Appendix A details the calculations made to determine misalignment errors due 
to the various contributors. Calculated values are based on the current design 
dimensions as shown in detail piece part and assembly level drawings. Total 
misalignment (shown as total allowable error) due to environment and 
manufacturing errors plus margin were run and the impact on on-axis encircled 
energy and rms spot size were calculated. These values are illustrated in the 
following error budget (Figure 5.1) and the top level number is the rss of the 
contributors. This in turn is related back to encircled energy to determine margin. 

Note that misalignment terms were calculated with respect to the secondary 
mirror rather than the focal plane. This results from the physical location of the SM 
in the telescope and the approach taken in assembling and aligning the optics. 
During assembly and alignment both the primary and tertiary mirrors are to be 
aligned to the secondary mirror. 



Ultra-Violet Imager Error Budget 
(On-Axis with -6° bias) 



Figure 5^1 UVI Error Budget Tree showing image degradation contributors in terms of rms spot size radi 
Based on initial MSFC thermal environment with preliminary margin allocations (see Appendix A). 










6. Conclusions 


As is seen in the error budget, the target of 60% encircled energy in a 27 pm 
radius is meet with some margin. This 4% margin, together with the margin 
initially allocated to each error contributor, helps ensure that any subsequent 
changes in the thermal environment or drastic changes in assembly philosophy or 
procedures will not require a relaxation of the target budget requirement. 


During recent conversations with MSFC, preliminary indications were provided 
that the thermal environment would be significantly less severe in terms of ground- 
to-orbit bulk temperature changes. This information will be folded back into the 
error budget when available. 
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APPENDIX A 

UVI Error Budget Calculation Details 


A-l 



Primary Mirror w.r.t. Secondary Mirror Despace Contributors 


Al(in.) = a(in./in/°F) Kin.) AT(°F) 


Bulk Temperature Change - Ground to Orbit AT = -36 °F (-20 °C) 


SM Body to Flange 
Flang to Housing Shim 
Housing 

PM Flang to Housing Shim 
PM Body to Flange 


Bulk Temperature Change - Orbit 


SM Body to Flange 
Flang to Housing Shim 
Housing 

PM Flang to Housing Shim 
PM Body to Flange 


Length 

Material 

CTE(in/in/°F) 

A (in.) 

0.449 

Al 6061-T651 

1.36E-05 

2.20E-04 

0.05 

AISI 410 

5.50E-06 

-9.90E-06 

3.182 

Al 6061-T651 

1.36E-05 

-1.56E-03 

0.05 

AISI 410 

5.50E-06 

-9.90E-06 

0.25 

Al 6061-T651 

1.36E-05 

1.22E-04 


A1 = - 0.00124 in. 

(- 0.032 mm) 


to Orbit AT = ±10 °F (±5.5 °C) 


Length 

Material 

0.4488 

Al 6061-T651 

0.0500 

AISI 410 

3.1804 

Al 6061-T651 

0.0500 

AISI 410 

0.2499 

Al 6061-T651 


CTE(in/in/°F) ±a (in.) 
1.36E-05 6.10E-05 

5.50E-06 -2.75E-06 

1.36E-05 -4.33E-04 

5.50E-06 -2.75E-06 

1.36E-05 3.40E-05 


Al = ± 0.00034 in. 


(± 0.009 mm) 


Manufacturing and Assembly Errors 


Measurement of Optical Surfaces to Reference Surfaces : 
Measurement of Spacing between Optics: 

Total Errors: 


±0.00014 in. 
±0.0005 in. 
±0.00052 in. 


A conservative total Primary Mirror w.r.t Secondary Mirror Despace Error is 
values: 


(includes both PM and SM) 
(±0.0132 mm) 
given by the rss of these 


£ = ± 0.0014 in. (± 0.035 mm) 


Total Allowable: £t = ± 0.002 in. (±0.050 mm) Margin: Em = ± 0.0014 in. (±0.035 mm) 


A-2 


Tertiary Mirror w.r.t. Secondary Mirror Despace Contributors 


AKin.) = CC(in./in/°F) Kin.) AT(°F) 


Bulk Temperature Change - Ground to Orbit AT = -36 °F (-20 °C) 



Length 

Material 

CTE(in/in/°F) 

A (in.) 

SM Body to Flange 

0.449 

Al 6061-T651 

1.36E-05 

2.20E-04 

Flang to Housing Shim 

0.05 

AISI 410 

5.50E-06 

-9.90E-06 

Housing 

7.521 

Al 6061-T651 

1A6E-05 

-3.68E-03 

TM Flang to Housing Shim 

0.05 

AISI 410 

5.50E-06 

-9.90E-06 

TM Body to Flange 

0.25 

Al 6061-T651 

1.36E-05 

1.22E-04 


Al = - 0.0034 in. 

(- 0.085 mm) 

Bulk Temperature Change - Orbit to Orbit AT = ±10 °F (±5.5 °C) 



Length 

Material 

CTE(in/in/°F) 

±A (in.) 

SM Body to Flange 

0.4488 

Al 6061-T651 

1.36E-05 

6.10E-05 

Flang to Housing Shim 

0.0500 

AISI 410 

5.50E-06 

-2.75E-06 

Housing 

7.5173 

Al 6061-T651 

1.36E-05 

-1.02E-03 

TM Flang to Housing Shim 

0.0500 

AISI 410 

5.50E-06 

-2.75E-06 

TM Body to Flange 

0.2499 

Al 6061-T651 

1.36E-05 

3.40E-05 


Al = ± 0.00093 in. 

(± 0.024 mm) 

Manufacturing and Assembly Errors 


Measurement of Optical Surfaces to Reference Surfaces : 
Measurement of Spacing between Optics: 

Total Errors: 


±0.00014 in. (includes both PM and SM) 
±0.0005 in. 

±0.00052 in. (±0.0132 mm) 


A conservative total Tertiary Mirror w.r.t Secondary Mirror Despace Error is given by the rss of these 
values: 


£ = ± 0.0035 in. (± 0.091 mm) 


Total Allowable: £t = ± 0.0047 in. (±0.120 mm) 


Margin: £m = ± 0.0029 in. (±0.073 mm) 


A-3 


Focal Plane w.r.t. Secondary Mirror Despace Contributors 


Al(in.) = a(in./in/°F) Kin.) AT(°F) 


Bulk Temperature Change - Ground to Orbit AT = -36 °F (-20 °C) 



Length 

Material 

CTE(in/in/°F) A (in.) 

SM Body to Flange 

0.449 

A1 6061-T651 

1A6E-05 2.20E-04 

Flang to Housing Shim 

0.05 

AISI 410 

5.50E-06 -9.90E-06 

Housing (to pin) 

0.165 

A1 6061-T651 

1.36E-05 -8.08E-05 

Top Plate to Det. Plate 

1.28 

A1 6061-T651 

1.36E-05 -6.27E-04 

Pin C/L to Det. Pad 

0.97 

A1 6061-T651 

1.36E-05 -4.75E-04 

Det. Shim 

0.1 

AISI 410 

5.50E-O6 -1.98E-05 

Det. Foot to FP 

1.31 

A1 6061-T651 

1.36E-05 6.41E-04 

Al= ±0.0004 in. 
(± 0.010 mm) 

Bulk Temperature Change - Orbit 

to Orbit AT 

= ±10 °F (±5.5 °C) 



Lengths 

Material 

CTE(in/in/°F) ±A (in.) 

SM Body to Flange 

0.4488 

A1 6061-T651 

1.36E-05 1.22E-05 

Flang to Housing Shim 

0.0500 

AISI 410 

5.50E-06 5.50E-07 

Housing (to pin) 

0.1649 

A1 6061-T651 

1.36E-05 4.49E-06 

Top Plate to Det. Plate 

1.2794 

AI 6061-T651 

1.36E-05 3.48E-05 

Pin C/L to Det. Pad 

0.9695 

A1 6061-T651 

1.36E-05 2.64E-05 

Det. Shim 

0.1000 

AISI 410 

5.50E-06 1.10E-06 

Det. Foot to FP 

1.3094 

Al 6061-T651 

1.36E-05 3.56E-05 

Al= ± 0.0001 in. 
(± 0.0025 mm) 


Manufacturing and Assembly Errors (estimated) 


Assembly of Detector Assembly: ±0.001 in. 


Total Errors: ±0.001 in. (±0.0254 mm) 


A conservative totalFocal Plane w.r.t Secondary Mirror Despace Error is given by the rss of these 
values: 


£ = ± 0.0011 in. (± 0.028 mm) 

Total Allowable: £t = ±0.0012 in. (±0.031 mm) Margin: £ m = ± 0.0005 in. (±0.013 mm) 

This is only an estimate, since it will be dictated by MSFC tolerances and assembly procedures. 


A-4 



Primary Mirror w.r.t Secondary Mirror Decenter Contributors 


Al(in.) = a(in./in/°F) l(in.) AT(°F) 


Bulk Temperature Change - Ground to Orbit AT = -36 °F (-20 °C) 


Length Material CTE(in/in/°F) 

SM to PM Decenter(Y-axis) 2.165 A1 6061-T651 1.36E-05 

SM to PM Decenter (Z-axis) 


Bulk Temperature Change - Orbit to Orbit AT = ±10 °F (±5.5 °C) 


A (in.) 
-1.06E-03 
0.0 


SM to PM Decenter (Y-axis) 
SM to PM Decenter (Z-axis) 


Length Material 

2.1639 A1 6061-T651 


CTE(in/in/°F) ± A (in.) 

1.36E-05 5.89E-05 

0.0 


Manufacturing and Assembly Errors 

Measurement of Vertex to Reference Surfaces (Y-Axis): 
Measurement of Vertex to Reference Surfaces (Z-Axis): 
Measurement of Spacing between Optics (Y-Axis): 
Measurement of Spacing between Optics (Z-Axis) 

Total Errors: 


±0.00014 in. (includes both PM and SM) 
±0.00014 in. (includes both PM and SM) 
±0.0005 in. 

±0.0005 in. 

±0.0007 in. (±0.0132 mm) 


A conservative estimate of the total Primary Mirror w.r.t. Secondary Mirror Decenter error is given by 
the rss of the above values: 

£ = ± 0.0013 (in.) (± 0.032 mm) 

Total Allowable: £t = ± 0.0015 in. (±0.038 mm) Margin: Em = ± 0.0007 in. (±0.019 mm) 


A-5 


Tertiary Mirror w.r.t Secondary Mirror Decenter Contributors 


Al(in.) = a(in./in/°F) l(in.) AT(°F) 


Bulk Temperature Qiange - Ground to Orbit AT = -36 °F (-20 °C) 


Length Material CTE(in/in/°F) A (in.) 

SM to TM Decenter (Y-Axis) 3.161 A1 6061-T651 1.36E-05 -1.55E-03 

SM to TM Decenter (Z-Axis) n n 


Bulk Temperature Change - Orbit to Orbit AT = ±10 °F (±5.5 °C) 


Length Material CTE(in/in/°F) ±A (in ) 

SM to TM Decenter(Y-Axis) 3.1595 A1 6061-T651 1.36E-05 8.59E-05 

SM to TM Decenter (Z-Axis) q q 

Manufacturing and Assembly Errors 

Measurement of Vertex to Reference Surfaces (Y-Axis): ±0.00014 in. (includes both PM and SM) 
Measurement of Vertex to Reference Surfaces (Z-Axis): ±0.00014 in. (includes both PM and SM) 
Measurement of Spacing between Optics (Y-Axis): ±0.0005 in. 

Measurement of Spacing between Optics (Z-Axis) ±0.0005 in. 

Total Errors: ±0.0007 in. (±0.0132 mm) 

A conservative estimate of the total Tertiary Mirror w.r.t. Secondary Mirror Decenter error for all axes 
is given by the rss of the above values: 


£ = ± 0.0017 (in.) (± 0.043 mm) 


Total Allowable: £t = ± 0.002 in. (±0.038 mm) Margin: £ m = ± 0.001 in. (±0.026 mm) 


A-6 


Mirror A Radius of Curvature Contributors 


AR (in.) = 0C(in./in/°F) R(in.) ATbulk(°F) 


AR (in.) = a(in./in/°F) R 2 (in. 2 ) (Tfront - Tback)(°F)/Thickness(in.) 


Bulk Temperature Change - Ground to Orbit AT = -36 °F (-20 °C) 


Primary Mirror Radius 
Secondary Mirror Radius 
Tertiary Mirror Radius 


R (in.) 
639200 
3.28600 
6.85900 


Material 
A1 6061-T651 
A1 6061-T651 
A1 6061-T651 


CTE(in/in/°F) 

1.36E-05 

1.36E-05 

1.36E-05 


AR (in.) 
-3.13E-03 
-1.61E-03 
-3.36E-03 


Bulk Temperature Change - Orbit 

to Orbit AT 

= ±10 °F (±53 °C) 




R (in.) 

Material 

CTE(in/in/ 0 F) 

± AR (in.) 

Primary Mirror Radius 

6.3889 

A1 6061-T651 

1.36E-05 

±1.74E-04 

Secondary Mirror Radius 

3.2844 

A1 6061-T651 

1 .36E-05 

±8.93E-05 

Tertiary Mirror Radius 

6.8556 

A1 6061-T651 

1.36E-05 

±1 .86E-04 


Axial Temperature Gradient (along optical axis) - On-Orbit 


AT/in. = ±03 °F/in. (±0.1 °C/cm) 


Primary Mirror Radius 
Secondary Mirror Radius 
Tertiary Mirror Radius 


R (in.) 

Ave. Thick, (in.) 

CTE(in/in/°F) 

± AR (in.) 

6.3889 

0.8250 

1.36E-05 

±2.78E-04 

3.2844 

03745 

1.36E-05 

±7.34E-05 

6.8556 

1.2650 

1.36E-05 

±3.20E-04 


Radial Temperature Gradients effects are to small to account for. 


A conservative total delta radius of curvature for each mirror is given 

Total 


Primary Mirror Radius 
Secondary Mirror Radius 
Tertiary Mirror Radius 


± AR (in.) ±AR(mm) 

±0.0031 ±0.079 

±0.0016 ±0.041 

±0.0034 ±0.086 


by the rss of these values: 

Allow, (mm) 

±0.81 

±0.42 

±0.88 


A -7 


Primary Mirror w.r.t Secondary Mirror Tilt Contributors 


Thermal gradients in the Y and Z axes will contribute to a relative tilt as well as alignment errors, a 
gradient along the optical axis will not contribute. It is assumed that a 0.5 °F/in through the optical 
bench in each of these axes is present. Note that the numbers attributed to thermal misalignment is 
very conservative given the hardware configuration. 


PM to SM Tilt (Y-axis): 3 arc-sec 

PM to SM Tilt (Z-axis): 2 arc-sec 

rss Total: 3.5 arc-sec 

Manufacturing and Assembly Errors 

Measurement of relative tilt in 2-axes: 30 arc-sec 

Total: 30.2 arc-sec 

Total Allowable: 108 arc-sec (1.8 arc-min) 

This conservative total allowable was allocated prior to details of assembly and hardware. 


A-8 


Tertiary Mirror w.r.t. Secondary Mirror Tilt Contributors 


Thermal gradients in the Y and Z axes will contribute to a relative tilt as well as alignment errors, a 
gradient along the optical axis will not contribute. It is assumed that a 0.5 °F/in through the optical 
bench in each of these axes is present. Note that the numbers attributed to thermal misalignment is 
very conservative given the hardware configuration. 


TM to SM Tilt (Y-axis): 12 arc-sec 

TM to SM Tilt (Z-axis): 15 arc-sec 

rss Total: 20 arc-sec 

Manufacturing and Assembly Errors 

Measurement of relative tilt in 2-axes: 20 arc-sec 

Total: 28 arc-sec 

Total Allowable: 60 arc-sec (1 arc-min) 

This conservative total allowable was allocated prior to details of assembly and hardware. 
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Focal Plane w.r.t. Secondary Mirror Tilt Contributors 


This contributor again is driven by MSFC assembly procedures. It is estimated that alignment can be 
done as well as the mirrors. 
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Primary Mirror w.r.t. Secondary Mirror Rotation Contributors 


This error contributor is driven by how well one can measure height variations of the precision 
alignment surface on the flanges of each surface and subsequently determining the angular tilt of the 
alignment surface. 

PM Alignment Surface Length: =1.4 in. (35.56 mm) 

Height Measurement Sensitivity: 0.0002 in. (0.0025 mm) 

Angular Resolution: 30 arc-sec 

Total Allowable: 2 arc-min 
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Tertiary Mirror w.r.t. Secondary Mirror Rotation Contributors 


This error contributor is driven by how well one can measure height variations of the precision 
alignment surface on the flanges of each surface and subsequently determining the angular tilt of the 
alignment surface. 

TM Alignment Surface Length: =2 3 in. (58.4 mm) 

Height Measurement Sensitivity: 0.0002 in. (0.0025 mm) 

Angular Resolution: 18 arc-sec 

Total Allowable: 2 arc-min 
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PDR TITLE 


ULTRAVIOLET IMAGER INVESTIGATION/POLAR 


PRELIMINARY DESIGN REVIEW OCTOBER 17, 1989 

MSFC Building 4481/Room 369 

PRELIMINARY AGENDA 

8:30am INTRODUCTION M.TORR 

introductions 
agenda/objectives 
instrument overview 


9:30am 


10:25am 

10:40am 


1 1 :30am 


12:00noon 


MECHANICAL DESIGN 

optical bench, mechanisms, 
housing, s/c interface, radiator, 
electronics stack housing 
(45 min) 

mass/c. g. status 

ICD summary (lOmin) 

BREAK 

THERMAL DESIGN 

camera, radiator, 
electronics stack, 
reduced thermal model 
(30 min) 

ICD summary (20 min) 

MISC. ICD SECTIONS 
magnetic interface 
electrostatic cleanliness 
environmental interface 
(30min) 

LUNCH 


L.WALKER 
L. WALKER 
J. SPANN 


L. WALKER 
J. SPANN 
J. SPANN 



1 :15pm 


ELECTRONICS DESIGN 

system, processor, memory, 

detector, detector interface, 

s/c interface, TLM format, 

housekeeping, 

mechanism control, power 

supplies, health and safety 

data (60min) L. SAVAGE 

power status L.SAVAGE 

NSPARS L. SAVAGE 


ICD summary (30min) J. SPANN 


2:45pm 


3:00pm 
3:1 5pm 


4:30pm 


5:00pm 


DEVELOPMENT STATUS AND 
SCHEDULES (15min) M.TORR 

BREAK 

SOFTWARE/DATA MANAGEMENT 
flight software (15min) 

GSE software (15min) 

Data Management Plan 
(15 min) 

data analysis software 
(30min) 


L. SAVAGE 
L.SAVAGE 
G. GERMANY 

D. TORR 


ISSUES, CONCERNS 
ACTION ITEMS 

AJOURN 
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ULTRAVIOLET IMAGER 




AGREED TO RESOURCES SUMMARY 


MASS 

19.5 KGMS 

POWER 

20 WATTS 

DATARATE 

12 KBPS 


PDR/RESOURCES 





UV IMAGER SCATTERING LOSS 



RMS ROUGHNESS (A) 















COLDPUTE 

(PRELIMINARY) 









ULTRAVIOLET IMAGER 



FIGURE 9-4 - DETAIL "B” - ELECTRONICS STACK 







HIGHLIGHTS OF THE ULTRAVIOLET 
IMAGER INVESTIGATION 


OPTICS 

■Past (f/3) system with wide field of view <8 degrees) 
and excellent iaage quality « 1 pixel over most of wage), 
state-of-the-art low scatter design 

DETECTORS 

proven intensif ied-CCD system that prowdes compact 2-D 
sensor with low noise and large dynamic range <4000 instantaneous 
selectable over a range of six to seven orders of magnitude) 

DATA SYSTEM 

on-board microprocessor (88086) controller and data formatter 
allows programmable operations and substantial reattiae data 
compression that is further enhanced by data coepression 
algorithm 

FILTERS 

narrow bandpass filters with excellent out of band rejection 
for Uey emissions <e.g. 1384A, t356A) allows the first truKj 
quantitative auroral imaging 

CALIBRATION FACILITIES 

existinQ vacuum ultraviolet catenation facilities allowing absolute 
catenation of imager at all wavelengths 

ATMOSPHERIC/AURORAL MODELS 

models of the upper atmosphere and energetic particle lapact 
on the atmosphere developed by our group over the last 
28 years are essential to the interpretation of the auT'oral 
information content 
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ISTP-UVI CAMERA MACS AND CG 
08:02 PM lS-Oct-89 

ITEM NAME 


PRIMARY DETECTOR 
SECONDARY DETECTOR 
MIRROR MECrl. ASSY 
CABLE HARNESSES 
FILTER WHEEL ASSY 
BAFFLE DOOR ASSY 
BASEPLATE ASSY 
CAMERA COVER ASSY 
OPTICAL BENCH ASSY 
THERMAL LINKS 
COLDPLATE ASSY 

INSTRUMENT TOTAL 


ISTP-UVI CAMERA MASS AND CG 
ITEM NAME 


PRIMARY DETECTOR 
SECONDARY DETECTOR 
MIRROR MECH. ASSY. 
CABLE HARNESSES 
FILTER WHEEL ASSY 
RFILTER/APER ASSY 
BASEPLATE ASSY 
CAMERA COVER ASSY 
OPTICAL BENCH ASSY 
THERMAL LINKS 
COLDPLATE ASSY 

INSTRUMENT TOTAL 


MASS XCG YCG ZCG 

(KG) (MM) (MM) (MM) 


1 

.Ill 

167. 

64 

234. 

48 

144 . 

*7 Q 
i u 

1 

X . 

.111 

114 

. 3 

215 

. 9 

144 . 

73 

0 

. 572 

287. 

02 

223 

52 

76 

. 2 

1 

X . 

. 706 

101 

. 6 

317 

. 5 

101 

. 6 

1 

.075 

325. 

12 

223 

. 6 

144 . 

n o 

f u 

0 . 

610 

266 

. 7 

279 

. 4 


27 

i . 

. 737 

266 

*7 
. 1 

228 

. 6 

7 . 

62 

0 . 

642 

266 

. 7 

228 

. U 

152 

. 4 

4. 

. 950 

213. 

36 

243. 

84 

144 . 

73 

0. 

200 

1 

27 

241 

o 

. -J 

135 . 

53 

1 . 

.030 

139 

. 7 

190 

. 5 

232 


4 . 

303 

205. 

26 

246 . 

60 

1 <J x . 

A 


MASS 

XCG 

YCG 

ZCG 

(LB) 

(IN) 

(IN) 

(IN) 

2.450 

6 . 6 

11.2 

5.7 

2.450 

4 . 5 

8 . 5 

5.7 

1 . 260 

11 . 3 

8.3 

o 

3.760 

4 

12.5 

4 

2.370 

12.3 

3 

<— — T 

D . ; 

1 . 245 

10 . 5 

x x 

c 

3 . 360 

10.5 

3 

0.3 

1.415 

10 . 5 

9 

6 

10.910 

8.4 

9.6 

5.7 

0.44 

5 

9 . 5 

7 . 7 

2.270 

5.5 

7 . 5 

11.5 

32.630 

8 . 0814 

9.7086 

5. 1591 
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TP-OVI ELECTRONICS MASS AND CG 
3:13 PM 16 -Oct-39 


ITEM NAME 

MASS 

XCG 

YCG 

ZCG 


(KG) 

(MM) 

(MM) 

(MM) 

CPU/SCIF A&B 

0.422 

139.7 

88.9 

101.6 

MEMORIES A&B 

0.462 

133.35 

38.9 

101.6 

DETECTOR I/F A&B 

0.422 

142.24 

88 . 9 

101 . 6 

HOUSEKEEPING BOARD 

0. 185 

139 . 7 

21 . 59 

101.fi 

TEU CONTROL BOARD 

0.185 

139.7 

147 . 32 

101.6 

DUAL MOTOR CONTROLLE 

0 . 256 

142.24 

170.13 

iOl . 6 

MOTHERBOARD ASSY 

0. 138 

41.275 

33.9 

101.6 

POWER SUPPLIES A&B 

0.556 

i o *7 
c* . 

O O Q 
U U . O 

101.6 

CRATE STACK 

1.234 

121 . 92 

38.9 

101.6 

STACK EHDPLATES 

0.596 

127 

88 . 9 

101.6 

TENSION RODS 

0. 149 

127 

88.9 

101.6 

FASTENER ALLOWANCE 

0.054 

127 

88 . 9 

101 . 6 

INSTRUMENT TOTAL 

4.718 113.5604 

92.95792 

101.6 


ISTP-UVI ELECTRONICS MASS AND CG 


ITEM NAME 

MASS 

XCG 

YCG 

ZCG 


(LB) 

(IN) 

(IN) 

(IN) 

CPU/SCIF A&B 

0.930 

5 . 5 

3 . 5 

4 

MEMORIES A&B 

1.018 

5 . 25 

3 . 5 

4 

DETECTOR I/F A&B 

0.930 

5.6 

3.5 

4 

HOUSEKEEPING BOARD 

0.407 

5 . 5 

0.85 

4 

TEU CONTROL BOARD 

0.407 

5 . 5 

5.3 

4 

DUAL MOTOR CONTROLLE 

0.564 

5 . 6 

6 . 7 

4 

MOTHERBOARD ASSY 

0.436 

1.625 

3.5 

4 

POWER SUPPLIES A&B 

1 . 225 

0 . 5 

3 . 5 

4 

CRATE STACK 

2.720 

4 . 3 

3.5 

4 

STACK ENDPLATES 

1.314 

5 

3 . 5 

4 

TENSION RODS 

0.329 

5 

3.5 

4 

FASTENER ALLOWANCE 

0.120 

C 

u 

3.5 

4 

INSTRUMENT TOTAL 

10.401 4. 

.470883 3. 

. 359760 
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UVI ICD MECHANICAL SECTION 


STATUS : 

Signed by GE and UVI 


COMMENTS : 

Instrument Mass 

Total mass does not include 
MLI blankets and intra- 
instrument harness 

Instrument Frequency 

Primary resonance 
frequency for both boxes 
expected to be >100 Hz 


Values for CG are available 



UVI ICD THERMAL SECTION 


STATUS : 

Inputs provided to GE 


COMMENTS : 

GE providing MLI blankets 


CONCERNS : 

Camera Operating Temperature 

ICD does not reflect current 
I IRD temperature ranges 

Thermal Backload Analysis 

Values of heat load on 
instrument surfaces are 
preliminary 

Heat Flow Limits 

Values for heat flow limits 
based on current operating 
temperatures and mounting 
surface contact area are not 
provided 



3. 


3 Thermal Interface Requirements 

3.3.1 General. The interface conditions specified in this section are 
intended to maintain each UVI item within the temperature range shown in 
Table 3.3-1 when its internal heat dissipation is within the range indicated. The 
Camera and electronics stack are isolated form the despun platform and the space 
environment with thermal isolators and thermal blankets except portions of the +Z 
facing surfaces are covered with selected coatings to radiate internal dissipation to 
space. 


Table 3.3-1 UVI Thermal Characteristics 


Item 


Operating 

Tmin/Tmax.C 


Qmin/Qmax.w 


Non-Operate 

Tmin/Tmax.C 


Qmin/Qmax.w 


Camera 

Camera (Radiator) 
Digital Elect Box 
Digital Elect (Radiator) 


-20/+20 

0/0 

-20/+40 

0/0 

-20/0 

10.0A0.0 

-20/+40 

0/0 

-20/+40 

0/0 

-20/+40 

0/0 

-20/0 

10.0/10.0 

-20/+40 

0/0 


NOTE: 


1 ) 


Cold plate temperature should be +/- 1°C when operating, load dissipated by ra- 
diator is a function of temperature. See Figure 10. 


2) C am era box dissipation is assumed to be dissipated entirely through radiator. 


3) Electronics box dissipation is assumed to be dissipated entirely through radiator. 


Unless otherwise noted, all interface conditions specified in this section 
are to be interpreted as orbit-average and spatial-average values. 


Size Code dent No. 

A 49671 


IS-3282161 


Thermal interfaces are specified in terms of parameters depicted schematically 
in Appendix 1. These parameters define temperature, heat flow, and coupling 
limi ts which should result in adequate thermal balance for all instr um ent 
units. Appendix 1 explains the derivation and use of the thermal backload parame- 
ter for external surfaces. 

3.3.2 Laborato ry Exterior . 

3. 3. 2.1 Externally Mounted Instruments. The UVT instr um ent is mounted on the 
despun platform of the POLAR laboratory and is classified as an externally 
mounted instrument. The following paragraphs specify the thermal interface 

for this unit. 

3.3.2.2 Instrument Exterior Thermal Environment. Radiation couplings and 
thermal backloads for the exterior surfaces of UVI sensor are specified in 
Table 3.3-2. These values are based on the surface finishes defined in 
paragraph 3.3.2.3. 

3.3.2. 3 Instrument Thermal Surface Finishes. Finishes and thermal 
properties of the exterior surfaces of the UVI sensor are specified in 
Table 3.3-3. 

3-3.2. 4 Conductive Coupling. The UVI Camera is mounted to but is thermally iso- 
lated from the despun platform. The UVI Electronics box is mounted to but is ther- 
mally isolated from the despun platform with GE-Astro supplied thermal isolators. 
Physical constraints on this interface shall apply as written in paragraph 3. 3. 3.4. 2 
of the GIIS. Table 3.3-4a specifies the nominal mounting area, conductance, and 
limits of despun platform temperature at this interface. Limits of heat flow by 
conductive interchange across this interface are specified in Table 3.3-4b. 


Size 

Code Idem No ! 


A 

49671 t 

IS-3282161 


Shee 


: 35" 


A 177 -- ro 


Table 3.3-3 UVT Exterior Surface Definition (TBR) 


ER UV Alpha 
Emia. BOM EOM 


Area 

Item Surface sa.in. Typ e 


Camera Sides 

TBD 

(Black body) 

TBD 

TBD 

TBD 

+Z face (top) 

TBD 

(Black body) 

TBD 

TBD 

TBD 

-Z face 

TBD 

(Black body) 

TBD 

TBD 

TBD 

E-Box Sides 

TBD 

(Black body) 

TBD 

TBD 

TBD 

+Z face (top) 

TBD 

(Black body) 

TBD 

TBD 

TBD 

-Z face 

TBD 

(Black body) 

TBD 

TBD 

TBD 

Table 3.3-4a 

UVI Exterior Mounting Thermal 

Interface, 

Tif (TBR) 



Unit 


Filler 


Contact 

Area.sa.in 


Conduct- Temperature (Min/Max),C 

ance.w/C Operating Non-Qperate 


Camera 
E Box 


Isolators 5x0.75 0.36 

Isolators 6x0.25 0.12 


-20/30 -30/40 

-20/30 -30/40 


Notes: 


1) Temperatures shown are on the spacecraft side of the interface. 


Table 3.3-4b UVT Exterior Mounting Conductive Limits (TBR) 


Heat Flow Limits, Qif (Note 2), w 

Operating Non-operate 


Unit 

Item 

£M 

Hot 


Hot 

Camera 

Isolated mount 

-3/+3 

-3/+3 

-4/+4 

-4/+4 

Digital Box 

Isolated mount 

- 21+2 

-21+2 

-3/+3 

-3/+3 


Notes: 

1) Heat flow is defined as positive INTO the instrument. 

2) Cold and hot refer to laboratory conditions defined in Table 3.3-4a. 

The range shown for each case assumes the instrument unit at its 
coldest/hottest allowable temperature. 


3.3.3 Laboratory Interior. N/A 


Size | Code ident No 

A I 49671 


IS-3282161 


Sheet 


Table 3.3-2 UVI Exterior Surface Thermal Environment (TBR) 

PRELIMINARY 

Thermal Backload, Qbl (Min/Max),w 


Item 

Surface 

Rad Coup, 

Sun Ang 
5 deg 

Sim Ang 
90 deg 

Sun Ang 
160 deg 

Camera 

Sides 

36.0 

4.0/ 4.9 

9.5A1.5 

3.2/ 4.0 


+Z face (top) 

36.0 

27.0/35.8 

0.0/ 2.0 

0.0/ 1.0 


-Z face 

N/A 

N/A 

N/A 

N/A 

E-box 

Sides 

36.0 

4.0/ 4.9 

9.5/11.5 

3.2/ 4.0 


+Z face (top) 

36.0 

27.0/35.8 

0.0/ 2.0 

0.0/ 1.0 


-Z face 

N/A 

N/A 

N/A 

N/A 


Notes: 

1 ) Radiation couplings are defined as surface area times surface emissivity times a 
view factor of unity. When used with thermal backload, these values should 
couple the surface to a 0 deg K heat sink and no other external couplings should 
be used. 

2) Thermal backload is the total heat absorbed by the surface from its external 
environment. It includes solar, albedo, and earth flux contributions, as well as 
ER input from nearby surfaces. 

3) The instrument will be operating at sun angles between 90 and 160 deg. It will, 
however, be exposed to sun angles between 25 and 160 deg in a non-operate 
mode. 

4) Coupling has been normalized for a 36 sq. in. surface, actual load is dependant 
on box surface area and scales directly with area. 

5) Initial values based on e and alpha values of 1.00 (blackbody). 
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UVI I CD EMC / EMI / ESC SECTIONS 


STATUS : 

First inputs to GE 


COMMENTS : 

Aperture Window 

Dielectric MgF 2 

window inside baffle is 

non-conducting 

Baffle Surface 

Interior baffle surface will 
be conducting 


CONCERNS : 

Degaussing 

Degaussing of instrument must 
not affect operation of 
permanent magnet stepper 
motors 

Surface Coatings 
Value for minimum 
conductivity is not specified 


UVI ICD ENVIRONMENTAL SECTION 


STATUS : 

First inputs to GE 


CONCERNS : 

Instrument Purge 

Nitrogen should be 99.999% 
pure with hydrocarbon content 
< 2 ppm (boil-off purity) 

Transportation and Storage 

Temperature limits should be 
+10/+30 °C 

Hydrocarbon content of 
surrounding atmosphere 
should be < 10 ppm 

Vibration, Shock and 
Acceleration Tests 

Design is based on a time 
consistent environment 



3.4 


3.4.1 Spacecraft Generated Magnetic Fields. The maximum spacecraft induced 
magnetic flux the instrument will experience on the spacecraft will be TBS gauss. 

3.4.2 Instrument Generated Magnetic Fields. Instrument magnetic field 
characteristics shall be documented here. (TBS by instruement) 

3.4.3 Instrument Degaussing. The instrument, unless noted here, will be 
degaussed prior to mounting on the spacecraft. The maximum field strength used 
during degaussing will be TBS gauss. The degaussing frequency will be TBS Hz. 

3.5 Electromagnetic Interference (EMI) 

3.5.1 General. Instruments with electrical devices inherently susceptible to low 
level EMI shall indicate special procedures and requirements here. 

3.6 Electrostatic Cleanliness (ESC) 

3.6.1 General. Bonding and insulation requirements of the instrument(s): TBD 

3.6.2 Conductivity and Grounding of Conductive Finishes. 

List of instrument coatings: TBS 
resistivity:TBS (ohms/meter.) 

3.6.3 Bonding of Case Parts. Physical discontinuities (non-conductive) in the 
instrument case must be identified here. 

3.6.4 Instrument Aperture ESC Design. See Instrument aperture drawing TBD. 

3.6.5 ESC of Hinged Mechanisms. Hinged mechanisms must be identified here. 

3.6.6 Exposed Connectors and Harnessing. Applies as written the GHS. 

3.7 Environmental Interface. 

3.7.1 General Environment. Applies as written the GIIS unless specified here. 

3.7.2 Storage. Transportation, and Handling Environment. Applies as written the 
GHS unless specified here. 

3. 7.2.1 Instrument Environment Before Mounting on Spacecraft. Applies as 
written the GIIS unless specified here. 

3. 7.2. 2 Instrument Stored Mounted on Spacecraft. Applies as written the GHS 
unless specified here. 


Size Code Gent Mo. 
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Sheer 



3. 7. 2. 3 Instrument Transportation while Mounted on Spacecraft. Applies as 
written the GIIS unless specified here. 

3. 7.2. 4 Instrument Purging. 

Purging requirements:TBD 
Connector type: TBD 

Purge Interface: See Mechanical Drawing for location. 

The instrument team will be notified of interruptions longer than TBD. 

3.7.3 Flight Environment. The instrument, while integrated on the satellite, shall 
withstand the flight-induced environment detailed in the following paragraphs. 

3. 7. 3.1 Temperature . The instrument shall specify nonoperating survival 
temperature range and expected operating temperature ranges here. 

NOTE: To ensure that the Qualification Temperature limits are not exceeded 
during emergency conditions during spacecraft thermal-vacuum testing, a "Safety 
Heater” and a monitoring thermocouple will be attached to the instrument unless 
directed otherwise here. 

The "Safety Heater" will be approved for thermal-vacuum chamber operation and 
will be powered form a +28V supply not connected to the spacecraft power. The 
capacity shall be sufficient to maintain the instrument at +10degC when the 
chamber walls are at LN2 temperature. Astro will mount the heater to the 
instrument when it is put onto the spacecraft. Heater location TBD. See figure 
TBD. 

3. 7. 3. 2 Thermal -Vacuum. The instrument shall operate within specification over 
the temperature range specified here while subjected to a nominal pressure of one 
(1) atmosphere or to a vacuum pressure of .75 x 10-5 Torr or lower. 

Min (1 Atm): TBD 
Max (lAtm): TBD 
Min (Vac): TBD 
Max (Vac): TBD 


3. 7. 3. 3 Vibration, Shock, and Acceleration. Applies as written the GIIS. 

3. 7. 3. 3.1 Sinusoidal Vibration. Expected test levels: TBD 

3. 7. 3. 3. 2 Shock. Applies as written the GIIS unless specified here. 

3. 7. 3. 3. 3 Acceleration. Applies as written the GnS unless specified here. 

3. 7. 3.4 Acoustic/Random Vibration. Applies as written the GIIS unless specified 
here. 
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3-6.4 Instrument Aperture ESC Design. A continuous conductive surface shall 
exist between the solar array substrate, which will be conductively bonded to the 
pr imar y structure, and instrument sensors which protrude through solar array 
apertures (see Paragraph 3.2. 6.3). The instrument surface required to interface 
with the gasket will be developed jointly by Astro and the Instrumenter. 

3-6.5 ESC of Hinged Mechanisms. Hinged mechanisms which do not rotate 
through 360° shall utilize a conductive strap connected across the hinge connected 
to spacecraft ground to preclude charge buildup on the mechanism. 

3-6.6 E xposed Connectors and Harnessing All external, exposed electric terminals 
shall be .rounded and coated with insulating material and, finally,' overcoated with a" 
conductive material. Such encapsulation shall not preclude the removal, 
interchange, additions, and/or repair of any wire/pin connection. All external 
instrument connectors shall have connector covers with a conducting outer surface. 

All harnessing external to the laboratory surface will be shielded in accordance with 
the GGS Electromagnetic and Magnetic Compatibility Plan, (Astro TBD) with the 
outer surface of the shielding coated with conductive finish per paragraph 3.6.2. 

3.7 Environmental Interface. 

All instruments will be exposed to the environmental conditions specified in the 
following paragraphs. These environments represent conditions which will arise 
during the transportation, storage, handling, test, launch and orbital operation of 
the instrument when on the spacecraft. 

3.7.1 General Environment. At the spacecraft contractor’s facility, controlled 
environments will be provided when necessary as specified in the Instrument- 
Unique ICD to bring the temperature, humidity, shock and vibration to levels less 
severe than those pertai nin g to launch, ascent, and orbital operations. 

3.7.2 Storage, Trans portation, and Handling Environment. 

3.7.2.1 Instrument Environment Before Mounting on Spacecraft. The 
environments experienced by the instrument during fabrication, storage, and all 
modes of ha n dl in g and transportation should be controlled so as to be sig nifican tly 
less severe than worst-case flight conditions. Storage will be in vendor's shipping 
container in a normal factory environment unless the Instrument-Unique ICD 
specifies other requirements. 

3.7.2.2 Instrument Stored Mounted on Spacecraft 


(1) Storage of the instrument when mounted on the spacecraft will be at 
Class 100,000 or better, as defined in FED-STD-209. 

During transportation and handling, the laboratory may temporarily be in 
areas not meeting class 100,000. During these periods, the spacecraft will 
be in a protective tent. When this tent cannot be used, the instrument 
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will be bagged to prevent its being cont amin ated 


(2) Humidity Limits of Storage Area: Maximum Humidity 50%. Laboratory 
may be stored in a tent purged with LN2 boiloff. 

(3) Temperature Limits: -10 to +40°C 

3.7.2.3 Instrument Transportation while Mounted on Spacecraft. The 
environments experienced during spacecraft transportation and handling will be 
controlled to be significantly less severe than worst-case flight conditions. 

3.7.2.4 Instrument Purging. If required, the instrument will be provided -a nearly 
continuous dry filtered nitrogen gas purge system, distributed through teflon tubes. 
The nitrogen will have a greater than 99.998% purity and less than 1 ppm of 
hydrocarbons. The purge gas flow will be available throughout integration at Astro 
and during launch site integration. Purge will be interrupted during thermal-vac 
testing and during some ground handling procedures. Instrument te ams will be 
notified of interruptions longer than specified in the instrument-unique ICD. 
Purging requirements of the instruments shall be specified in the Instrument- 
Unique ICD. 

3.7.3 Flight Environment. All instruments integrated on the satellite shall be 
designed to withstand the flight-induced environment detailed in the following 
paragraphs. 

> 3.7.3. 1 Temperature. All instruments installed on the spacecraft shall specify 

nonoperating survival temperature range and expected operating temperature 
ranges in their respective Instrument-Unique ICD. 

NOTE: To ensure that the Qualification Temperature limi ts are not exceeded 

during emergency conditions during spacecraft thermal-vacu um 
testing, a "Safety Heater" and a monitoring thermocouple will be 
attached to the instrument 'unless directed otherwise in the 
Instrument Unique ICD. 

The "Safety Heater" will be approved for thermal-vacuum chamber operation and 
will be powered form a +28V supply not connected to the spacecraft power. The 
capacity shall be sufficient to maintain the instrument at +10°C when the chamber 
walls are at LN2 temperature. Astro will mount the heater to the instrument when 
it is put onto the spacecraft. 

The thermocouple shall be copper-constantan, mounted close to the location of the 
ins tr um ent's analog temperature telemetry sensor. The leads of this TC shall be 
independent of the S/C-instrument harness. Sensor leads inside the instrument 
case shall be properly secured so that only the flying leads need to be clipped during 
preflight preparations. 

3.7.3. 2 Therm al-Vec*:T!m. The instr um ent shall operate wi thin specification over 
the temperature range specified in the Instrument Unique ICD while subjected to a 
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nominal pressure of one (1) atmosphere or to a vacuum pressure of .75 x 10-5 Torr 
or lower. In addition, the instrument shall be capable of operating after a rate of 
pressure change of 1.27 psia/sec. The full pressure rate of change is shown ip,* 

Table 7. 

/£U3.3 Vibration. Shock, and Acceleration. All instruments shall be designed to 
withstand thewibration, shock, and acceleration environments experienced during 
test, launch and orbital operation and shall operate in accordance with the 
instrument performance specification requirements for those environments. 

3.7.3. 3.1 Sinusoidal Vibration. The sinusoidal vibration qualification levels 
experienced by the protoflight spacecraft in each of three orthogonal axes -are listed • 
in Table 8. The levels experienced by each instrument during these tests shall be 
specified in each Instrument-Unique ICD. 

3.7.3. 3.2 Shock. The spacecraft will experience a shock impulse when the separa- 
tion band is released during test as well as during separation in orbit. A test of this 
event will be performed on the integrated spacecraft. The maximum expected flight 
shock levels are shown in Figure 31 . 

3. 7.3. 3. 3 Acceleration. The instruments will experience worst case expected steady- 
state accelerations during launch as delineated in Table 9 for platform and body 
mounted WIND and POLAR instruments. 

3. 7.3. 4 Acnustic/Random Vibration. The instrument, while mounted on the 
spacecraft, shall be capable of withstanding and shall operate within specification 
after a single exposure of a 1.0 min. duration to the acoustic excitation levels 
defined in Figure 32. An acoustic test may also be required if evaluation indicates 
sensitivity to direct acoustic energy. The instrument shall survive and operate 
within specification after exposure to the random vibration levels specified in 
Table 10. 

3.7.4 Radiation Environment. Preliminary radiation levels for the orbital 
environments of the WIND and POLAR laboratories are shown in Figure 33 and 
Figure 34. Instruments located within the laboratory will experience protective 
shielding provided by the laboratory outer shell, as indicated by the effective 
shielding values in Figure 33 and Figure 34. Instruments, or portions of 
instruments adjacent to openings in the laboratory shell, or otherwise exposed to 
the external environment, will receive radiation dosages dependent upon the degree 
of exposure and the equivalent thickness of protection provided by the instrument 
housing. 


Code (dene No. 


A 49671 


32S206S 


Sheet 



> 

TABLE 8. SINUSOIDAL VIBRATION CRITERIA PROTOFLIGHT LEVELS 



FREQUENCY (Hz) ACCELERATION (g's zero-to-peak) 


Thrust Axis 5 to 6.8 

6.8 to 30 
30 to 40 
40 to 100 

Lateral Axis 5 to 6.2 0.5 in. double amplitude 

6.2 to 100 1.0 

Sweep Rate = 4 oct./min. 


0.5 in. double amplitude 
1.2 

1.5 - 

1.0 



X 


TABLE 9: INSTRUMENT DESIGN LOADS (g's) 1 - 2 


Instrument Mounting 
Location 

X 

Y 

Z t > 

Despun Platform 3 
(or +Z Exterior Sensors) 

15.2 

14.7 

13. 

Body Mounted on Equipment 
Decks 

5.9 

5.9 

13.0 

1. Referenced to S/C Coordinate System (See Figure 20). 



2. Pre liminar y results based on Delta II Inputs and S/C Structure Coupling 
Effects. 


3. Apulies to: POLAR - UVI, VIS, PIXIE, SEPS, DRT 

WIND - TGRS, KONUS, 3-D PLASMA 
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UVI ICD GSE SECTION 


STATUS : 

First inputs to GE 


CONCERNS : 

S/C Interface Emulator 

Verification of S/C 
interface emulator hardware 
and software 

Synchronization Signals 

Data rate signal lacking 

Environment 

Temperature limits should be 
20°C to 30°C 

Relative humidity limits 
should be 40% to 55% 



4.0 INSTRUMENT GROUND SUPPORT EQUIPMENT (GSE). 

Instrument GSE designed and fabricated primarily for use in tests performed at 
Astro facilities or a test site shall be capable of: 

a) Verifying that the instrument has survived shipment; and 

b) Demonstrating successful completion of spacecraft systems test 
requirements. 

All Instrument GSE required to electrically interface with the GGS spacecraft 
checkout station (SCS) shall be provided the interfaces specified in Section 4.1. 

4JL Instrument GSE Interfaces 

The Instrument GSE test sets shall interface with the laboratory GSE as illustrated 
in Figure 35, and described in the sections that follow. In general, the Instrument 
GSE shall monitor instrument performance via laboratory science data supplied in 
real time, or as played back from a laboratory GSE tape recorder. Instrument GSE 
processing of this data shall not require on-line support from any laboratory GSE 
computer. All Instrument GSE interface requirements shall be documented in the 
Instrument-Unique ICD. 


4.1.1 Power. The instrument GSE test set shall operate on single phase 120 -fl2/-S 
volt ac power supplied at 57 to 63 Hz and shall draw a marrimum of 20 amperes 
current. 


4.1.2 Instrument Data. The instrument GSE will be provided both low- rate 
engineering and science laboratory downlink telemetry and high-rate instrument 
playback downlink telemetry. The telemetry will be decoded biphase-L or NRZ-L 
PCM serial data, at bit rates equivalent to laboratory downlink rates. The data 
shall be received by the instrument GSE on a differential line driver/receiver 
interface per Figure 36. 

yy' 4J..3 Synchronization Signals. The instrument GSE shall be provided major and 
f minor frame sync pulses and a data bit sync clock for processing the engineering 
Vjmd science data of Section 4.1.2. Low- rate data shall be clocked at a bit rate equal 
to the real time rate for each laboratory. Kigh-rate instrument playback data shall 
be synchronized with a 500 kbps clock. All clock and frame sync signals shall be 
received on a differential line driver/receiver interface per Figure 36. 


4.1.4 Local Area Network Interface ('Optional! At the option of the instrument 
contractor, a local area network interface, e.g. Ethernet, also will be provided to 
•instrument GSE test sets to allow the SCS to update its command load database 
from the instrument GSE. The protocol for this command file retrieval function will 
be established at System Requirements Review. 
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Figure 35. Instrument GSE Interfaces at Astro 
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Safety engineering principles should be applied in accordance with MIL-S-39130. 
Particular attention shall be paid to fail-safe provisions in the equipment design to 
avoid damage to other laboratory GSE components in the event of catastrophic 
failure. 

4.3 Instrument GSE Environmental Requirements 

The Instrument GSE shall be capable of operation over a temperature range of 55° 
to 95°F and a relative humidity range of 20 to 70%. 

4.4 Instrument GSE Complement 

The instrument contractor shall supply Astro with a detailed list of the Instrument 
GSE to be used at Astro at the Instrument PDR. The list shall include: 


(a) Type of Instrument GSE 

(b) Function of Instrument GSE 

(c) Power requirements 

(d) Approximate sire and weight of racks, targets, etc. 
5.0 NOTES 


5.1 Acronyms 

► 

In addition to common acronyms, the special acronyms listed below are used 
throughout this document. 


A/D 

ATM2 

C&DH 

CDR 

CEI 

CDU 

COSMOS 

D/A 

EMI 

FAT 

FCC 

FSP 

GFE 

.GSE 

GSFC 

GTM 

PC 

PDR 

RFEM 

R&QAE 


Analog to Digital 
Ampere-Turn Meter Squared 
Command and Data Handling 
Critical Design Review 
Contract End Item 
Control Distribution Unit 

Complementary SymmetryMetal Oxide Semiconductor 

Digital to Analog 

Electromagnetic Interference 

Flight Acceptance Test 

Flat Conductor Cable 

Frame Sync Pulse 

Government Furnished Equipment Instrument 

Ground Support Equipment 

Goddard Space Flight Center 

GGS Telemetry Module 

Power Converter 

Preliminary Design Review 

Reduced Finite Element Model 

Reliability and Quality Assurance Engineering 
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GGS/UVI STACK CPU & S.C. INTERFACE DIAGRAM 
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GGS/UVI STACK MOTOR CONTROL DIAGRAM 13 
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Schaeffer Magnetics. Inc. 


TYPE 1 

JZ y* V ' ( - « -I . *W i «l w — •, , V ^ 


Specificallv developed for nigh reliability space flignt applications: Schceffer Mcgnetics Rotary 
incremental Actuators employ a small cngie cermcnenr magnet stepper motor to airectly drive an 
integral harmonic drive speed reducer. 
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DIAGRAM 
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DEPLOYMENT MECHANISM 
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UVI ICD ELECTRICAL SECTION 


STATUS : 

First inputs to GE 


COMMENTS : 

Cable Shielding 

Clarification of shielding 
in Electronics Stack for 
Camera S/C interface signals 

Command Timming 

On which edge of clock is 
command clocked in? 

Major Frame Telemetry Format 

Clarification of S/C 
housekeeping and instrument 
housekeeping formats 

Power Supply Frequency 

125 KHz 

Time Code in Command Data 

Required every 10 min 



UVI ICD ELECTRICAL SECTION 


CONCERNS : 

Major Frame Telemetry Format 

Contigious word enable in 
instrument houskeeping 
telemetry 



UVI S/C INTERFACE 


S/C 



+28 Volt (regulated) 

+28 Volt (pulse) 

Command Enable 

Command Clock 

Command Data 

Telemetry Clock 

Telemetry Enable 

Major Frame Sync 

Minor Frame Sync 

< Telemetry Data 

Platform Status 


< Camera Radiator Temp 

<-- Elec Stack Radiator Temp — 

< 0° Aperture Temp 

< 90° Aperture Temp 

< 180° Aperture Temp 

< 270° Aperture Temp 

< Aperture Door Status 

< — Filter Wheel Home Status — 



Note: Interface is redundant 





X 


All of the above grounds will be returned as twisted pairs with their power lines 
through separate leads (in the spacecraft harness) to a single ground tie-point 
herein after referred to as "laboratory" or "star" ground. All command, data,. and 
telemetry voltage levels shall be referenced to "laboratory" ground. For co mm and 
and data signals, the ground level offset from instrument to spacecraft "star" 
ground shall be assumed to be 0.1 volt maximum (0.04 V typical). See Figure 6 for 
the GGS grounding philosophy. The telemetry ground return length will be no 
longer than 20 feet of AWG 24 stranded wire, (28.1 ohms per thousand feet). 

3.1.1.2 -i-28 Volt Main Power Ground. The +28 Volt Main Power ground shall 

be used as the return for all current drawn from the spacecraft +28 volt main 
regulated bus. '* " 

Current in the power ground shall be limited to 500 mA per #22 conductor or 1 .Amp 
per allocated connector pin (see paragraph 3.1 .2.3.2) excluding transient current 
drawn during instrument turn on in order to minimize the effects of harness drops. 

Feedthrough filters must be grounded to chasses to maintain box shielding effec- 
tiveness and must have a capacitance value less than 0.1 pfd to limit AC currents 
flowing to structure from the power bus. 

3.1 .1.3 -^28 Volt Pulse Load Ground. The Pulse Load ground shall be used as the 
power return line for all pulse loads (steppers, heaters, etc.) which do not comply 
with the main bus ripple specification (see paragraph 3.1.3.2.6). Current in this 
power ground shall be limi ted to 500 mA per #22 conductor or 1A per allocated con- 
nector pin. 

Feedthrough filters must be grounded to chasses to maintain box shielding effec- 
tiveness and must have a capacitance value less than 0.1 jifd to limit AC currents 
flowing to structure from the pulse power bus. 

3.1 .1.4 -+5V Keep-Alive Power Ground. The +5 volt memory keep-alive power shall 
be returned over the signal ground pin. 

3.1.1.5 Signal Ground. Signal ground shall be the power return line for the 
secondary side of the instrument dc/dc converters. Current in the signal ground 
output of the experiment shall not be more than 75 mA per connector pin. 

The signal ground may be connected to chassis by audio frequency/radio frequency 
(AF/RF) bypass capacitors (less then 0.5 pfd total) with short leads to minimize the 
effects of impedance of the signal ground lead at high frequencies. 

3.1 .1.6 Shield Grounding. For each preassigned shielded cable that carries a signal 
into the instrument, a connect.. - pin shall be dedicated to that cable's shield. The 
connector pins so dedicated shall not be connected, but held in reserve. No other 
signal returns or grounds shall be connected to the shield ground pins. Shield 
ground shall not be used with triax cables. Except for coaxial cable, cable shields 
shall not be used as a return path. 
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For non-RF signal cables, the shield shall be grounded at the signal source end only. 
The shield will be terminated at the connector backshell (chassis) and also carried 
through a pin into the ins tr um ent. Each signal connector must provide at least one 
chassis grounded pin. The chassis ground pin must be bonded internally to the 
equipment chassis 

3.1.1. 7 Chassis Ground. Chassis (case) ground shall be provided on the same 
connector (see Section 3.1 .2.3.6) used to interface the spacecraft power input with 
the instrument. Chassis ground shall be dc connected to the case (including outside 
wrapper and mo untin g surface) of the instrument. Bond straps (jumpers) must be - 
used across all hinges and connections that would not otherwise meet the -bonding 
requirement. 

The chassis ground connection to the spacecraft will be made through the mounting 
surfaces and mounting hardware of the instrument. The instrument mount shall be 
designed to provide a total contact resistance of less than 0.0025 ohm from the 
chassis to the spacecraft. Provisions shall be made for the installation of a ground 
strap if required. 

The shells of all connectors interfacing with the spacecraft harness shall be 
grounded to the chassis of the instrument. Conductive harness overwrap will be 
connected to the connector shell on the harness side. 

3.1.1 .8 Pyrotechnic Return. A separate pyrotechnic return will be supplied for each 
pyrotechnic device. It shall be separated from all other grounds within the unit. 

3.1 .1.9 Passive Analog Telemetry Ground. Passive analog telemetry points (ther- 
mistors and multi- turn potentiometers) shall have a common return line to the 
GTM. 

3.1.2 Wiring and Connectors 

3.1 .2.1 Harness Design Philosophy. The laboratory harness will be designed in 
accordance with the following guidelines: 

(a) Harness runs will be selected to minimize overall length and weight.' 

(b) Power returns will be routed in the same bundle as the associated power 
lines and returned to the source in order to minimize current loops. 

(c) Telemetry (GTM) signals will be bundled together and in a separate bundle 
from power and digital command lines in order to minimize crosstalk. 

(d) Signal returns will be routed in the same bundle as the associated signal 
lines in order to minimize electromagnetic emission and susceptibility. 
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Figure 17. Serial Command Phasing Relationship 
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UVI PDR/ DESIGN STATUS 


SUBSYSTEM 

ELEMENT 

DESIGN 

COMPLETE 

ENG MOD 
FABRICATION 

FLIGHT 

FAB 

OPTICS 

MIRRORS 

100% 

AWAITING 

MACHINING 

7/1/90 


FILTERS 

70% 

IN PROCESS 

9/1/90 

MACHINING 

DETECTORS 

COMPLETE 

IN PROCESS 

1 0/30/90 


FILTER WHEEL 

50% 

1 /5 

1 0/30/90 


FOLDING MIRROR 

10% 

3/5 

1 0/30/90 


MIRROR/DOOR 

10% 

2/1 

1 0/30/90 


BAFFLES 

80% 

2/1 

1 0/30/90 


OPTICAL BENCH 

90% 

2/1 

1 0/30/90 


BASE PLATE 

70% 

2/1 

1 0/30/90 


ELECTRONICS BOX 

30% 

5/1 0 

1 0/30/90 


PASSIVE 

RADIATOR 

15% 

4/5 

1 0/30/90 

ELECTRONICS 

DETECTOR A 

95% 

10/30 

1 0/30/89 


DETECTOR B 

95% 

10/30 

1 0/30/90 


DETECTOR 

INTERFACE 

40% 

1/5/90 

1 0/30/90 


S/C INTERFACE 

50% 

12/1 5/89 

1 0/30/90 


TEC 

50% 

1 / 5 / 9 0 

1 0/30/90 


HOUSEKEEPING 

50% 

2/1/90 

10/30/90 


CPU 

70% 

1 2/1/89 

1 0/30/90 


MEMORY 

70% 

1/5/89 

1 0/30/90 


MOTOR DRIVERS 

70% 

12/15/89 

1 0/30/90 


POWER SUPPLIES 

80% 

1 1/1 /89 

1 0/30/89 


BACK PLANE 

50% 

12/1 5/89 

1 0/30/90 


ENCODER BOARDS 

10% 

1/15/90 

1 0/30/90 

SOFTWARE 

FLIGHT 

10% 

2/1/90 

1/1/91 


GSE 

10% 

2/1/90 

1/1/91 


KEY PARAM/TLM 

0% 

7/1/90 

7/1/91 


ALGORITHMS 

10% 


1 0/1/91 




PDR TITLE/006 


vil FLIGHT SOFTWARE DEV. DIAGRAM 19 

ight software development configuration 



SIMULATOR I MW 



GGS/UVI FLIGHT SOFTWARE REQUIREMENTS 


DIAGRAM 2 


1. Program code shall not utilize more than 8k words of PROM 
space . 

2. Program code shall not utilize more than 256k of RAM memory. 

3. A means shall be provided by the flight program for uploading 
program "patches" in the event of a software or hardware 
malfunction or failure. This implies that the program code 
shall allow portions of the main flight program to reside 

in RAM memory. 

4. The code shall be structured in such a manner as to circumvent 
and report minor faults where possible. 

5. The flight program shall implement a power on self test (POST). 

The POST routine shall verify instrument functional 1 ity . 

6. The software program architecture shall be a table driven state 
machine. This software structure optimizes code size and process 
execution time. 
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PROGRAM 

CODE 



uv imager 



uu imager 








uu imager 



uv imager 



integer 


uu imager 







- ORBIT 


CCD READOUT 


SUM OR SAMPLE 
Reduce numoer 0 / pixels by following 
sequence of forming neighborhood sums 
or of decimating image by sampling 


ENCODE VALUES OF 
REMAINING PIXELS 
Reduce number of bits/pixel 



GROUND PROCESSING 


TRANSMIT 


GSFC RECEIVES 
RAW DATA 


CDHF PRODUCES 
LEVEL 1 DATA 


DECODE 

Restore neignoorhood or sampiea 
values from received coded numoers 



ARCHIVE 

LEVEL 


CORRECTIONS 
CCD response corrections 
Motion comoensaDon if necessary 
Bacxground subtraction 
Aspect rano adjustment 
NS division by 4 


QUICK 

LOOK 


DETERMINE 
TYPE OF 
IMAGE 


CREATE KPF 


DERIVE 

HIGH RESOLUTION 
IMAGES 


EXPAND 

NEIGHBORHOOD SUM 
IMAGES 


ARCHIVE KPF 
LEVEL 2 


I wwwwwwxwww M 


ARCHIVE 
LEVEL 2 




SELECT IMAGES AND 


PERFORM DESIRED 


SCIENTIFIC 


— hi m 11 —1 


CDHF a Central Data Handling Facility, GSFC a Goddard Space Right Center, KPF a Kay Parameter File 

NS 3 Neighborhood Sum, SG a Sampled Grid 


OVERVIEW OF DATAFLOW FROM ULTRAVIOLET IMAGER IN NOMINAL OBSERVING SEQUENCE 

















uv imager 




uv imager 



uv imager 
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uu imager 




DATA ANALYSIS SOFTWARE 


THREE MAIN CATEGORIES: 

1 . KEY PARAMETERS 

THIS PROVIDES SUMMARY DATA USEFUL TO THE 
ENTIRE ISTP COMMUNITY 

2. FULL SCALE DATA REDUCTION 

REDUCTION OF SPECTRAL IMAGES TO GEOPHYSICAL 
QUANTITIES 

3. SCIENTIFIC MODELING OF OBSERVATIONS 

INTERPRETATION OF MEASURMEENTS BY GLOBAL 
MODELING OF EMISSIONS 

INPUT CALIBRATED RECONSTRUCTED SPECTRAL IMAGES 


BRIGHTNESS GIVEN IN RAYLEIGHS 



KEY PARAMETERS 

1. TOTAL GLOBAL POWER INPUT AND AVERAGE ENERGY FLUX 

2. SIZE AND SHAPE OF THE POLAR CAP 

3. SIZE AND SHAPE OF THE AURORAL OVAL 

4 . TOTAL POWER INPUT AND AVERAGE ENERGY FLUX TO 
THE POLAR CAP 

5. TOTAL POWER INFLUX AND AVERAGE ENERGY FLUX TO THE 
AURORAL OVAL 

6. ACTIVITY INDICES 

7. CHARACTERISTIC ENERGY ESTIMATE AT SELECT LOCATIONS 


8 . 


SELECTED AURORAL IMAGES 



INVESTIGATIVE METHOD 


THE UVI IMAGER WILL OBTAIN GLOBAL IMAGES OF THE EARTH 
AT SEVERAL TTT.TRAVIOLET WAVELENGTHS . THESE WAVELENGTHS 
WERE SELECTED SPECIFICALLY TO PROVIDE THE INFORMATION 
NEEDED TO EXTRACT THE KEY PARAMETERS 

THE WAVELENGTHS SELECTED DEFINE EMISSIONS THAT ARE 
GENERATED: 

- BY DIFFERENT ATMOSPHERIC CONSTITUENTS 

- AT DIFFERENT ALTITUDES 

- WITH DIFFERENT LEVELS OF ABSORPTION 


TASK: OBTAIN INTENSITY AT THE 

REQUIRED WAVELENGTHS 


* 


SPECTRAL DECONVOLUTION SOFTWARE 



UVI FILTER WAVELENGTHS AND FUNCTION 


X (nm) 

SOURCE SPECIES 

FUNCTION 

121.6 

H L a 


PROTONS & GEOCORONA 

130.4 

0 


0 ABUNDANCE 

135.6 

0 


ELECTRON ENERGY SPECTRUM 

150 (BROAD) 

N 2 (LBH 

SHORT) 

ELECTRON ENERGY SPECTRUM 

150 (NARROW) 

N/N 2 


SPECTRAL DECONVOLUTION 

165 (BROAD) 

N 2 (LBH 

MID) 

SPECTRAL DECONVOLUTION 
ELECTRON ENERGY SPECTRUM 

185 (BROAD) 

N j (LBH 

LONG) 

TOTAL ENERGY INPUT 
ELECTRON ENERGY SPECTRUM 

215 

NO 


NITRIC OXIDE ABUNDANCE 


IMAGE RECOGNITION PACKAGE 


IRP 


THIS SOFTWARE PACKAGE IS DESIGNED TO RECOGNIZE THE INNER 
AND OUTER BOUNDARIES OF THE AURORAL OVAL. 

TECHNIQUE INTENSITY THRESHOLD/ GRADIENT DETECTION 

RMS FITTING OF BOUNDARIES WITH SNAKE FUNCTIONS OR 
SIMILAR TECHNIQUES 

THE IRP PACKAGE WILL PROVIDE THE KEY PARAMETERS: 

- SIZE (AND SHAPE) OF THE AURORAL OVAL 

- SIZE (AND SHAPE) OF THE POLAR CAP 

FROM THE SIZE AND SHAPE, THE TOTAL ENERGY INFLUX 
PACKAGE (TEP) WILL GENERATE THE KEY PARAMETERS: 

- TOTAL ENERGY INFLUX TO THE AURORAL OVAL 

- TOTAL ENERGY INFLUX TO THE POLAR CAP 




TOTAL ENERGY PACKAGE 


TEP 


THE TEP PACKAGE USES OUTPUT FROM IRP (IMAGE RECOGNITION 
PACKAGE) TO OBTAIN A MAP OF THOSE PIXELS COMPRISING THE 
AURORAL OVAL AND THE POLAR CAP. 

TEP THEN CALLS EFF (ELECTRON FLUX FUNCTION) TO OBTAIN THE 
ENERGY FLUX AT EACH PIXEL. 

TEP SUMS THE ENERGY FLUX PER PIXEL WITHIN THE DEFINED 
PIXEL MAP. 

FROM THE FOLLOWING MAPS TEP WILL PROVIDE 

GLOBAL MAP: TOTAL GLOBAL POWER INPUT DUE TO ELECTRONS 

AURORAL OVAL MAP: TOTAL AURORAL POWER INPUT DUE TO 

ELECTRONS 

POLAR CAP MAP: TOTAL POLAR CAP POWER INPUT DUE TO ELECTRONS 

ENERGY FLUX FUNCTION [EFF] 

THE EFF ALGORITHM WILL COMPUTE THE ENERGY FLUX GIVEN THE 


MEASURED LBH-LONG INTENSITY 




SPECTRAL DECONVOLUTION PACKAGE 


SPECDEC 


UNDER THE UVI PROGRAM AN ENTIRE NEW TECHNOLOGY WAS 
DEVELOPED TO BUILD FILTERS WITH SUFFICIENTLY NARROW 
BANDBASS AS TO BE ABLE TO ISOLATE KEY SPECTRAL LINES. 
BECAUSE OF THE PROXIMITY OF SOME OF THE SPECTRAL FEATURES 
SOME DECONVOLUTION IS REQUIRED. THIS PROCESS UTILIZES 
THE SPECTRAL INFORMATION FROM SEVERAL FILTERS TO ISOLATE 


INDIVIDUAL LINES 




GRIP POINTS ROUTINE 


GP 


THIS ROUTINE IS USED TO DEFINE THE LOCATIONS AT WHICH 
THE CHARACTERISTIC ENERGIES OF THE PRECIPITATING ELECTRONS 
ARE TO BE COMPUTED. THE DEFAULT GRID WILL BE AXES 
CENTERED ON THE DAWN-DUSK AND NOON-MIGNIGHT GEOMAGNETIC 
MERIDIANS . 

THE GRID SPACING WILL BE 1 # TO 5 4 IN LATITUDE WITH 
SPECIFIED HIGH-LOW BOUNDARIES 



CHARACTERISTIC ELECTRON ENERGY 


CHEE 


THE CHEE PACKAGE WILL COMPUTE THE CHARACTERISTIC ENERGY OF 
PRECIPITATING ELECTRONS AT A SPECIFIED SET OF GRID POINTS. 
DEFINED BY THE GP ROUTINE 

THE DEFAULT GRID WILL BE ALONG AXES SAMPLING FOUR 
SECTORS COVERING THE DAWN-DUSK AND NOON-MIDNIGHT MERIDIANS 
AT 'I TO 5 0 IN LATITUDE 

THIS ROUTINE WILL COMPUTE THE RATIOS: 

I1356/LBH-LONG LBH-SHORT/LBH-LONG 

IT WILL USE LOOK UP TABLES ILLUSTRATED IN FIGURES 1 
AND 2 TO ESTIMATE THE CHARACTERISTIC ENERGY AT THE 


PRESCRIBED GRID POINTS 




Energy (kev) 


Intensity Ratio (01 1356/LBH 1838) 



Figure 1 


100 




Ratio (LBH 1838/LBH 1464) 


o 

» 

"* -*■ o 



Figure 2 




SELECTED AURORAL IMAGE PACKAGE 


SAI 


THE SAI PACKAGE WILL GENERATE A FULL INTENSITY COLOR 
PLOT OF THE LBH-LONG EMISSION ( a TOTAL ENERGY FLUX) 
APPROXIMATELY ONCE EVERY 10 MINUTES. 


ACTIVITY INDICES 


AI 


THE AI PACKAGE WILL SCAN THE LBH-LONG IMAGE AT INTERVALS 
OF 10 MINUTES AND COMPUTE THE PERCENTAGE AREA IN FOUR 
SECTORS FOR WHICH THE IMAGE BRIGHTNESS EXCEEDS GIVEN 
THRESHOLDS, E.G. 1KR, 10KR, 100 KR. 



KEY PARAMETER FLOW CHART 


READ INCOMING DATA STREAM 

l 

STRIP OUT DATA FOR KEY PARAMETERS 
CYCLE TIME 10 MINUTES 
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CALL SAI 


SELECTED AURORAL 
IMAGE (LBH-LONG) 
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CALL EFF 
ENERGY FLUX 
PER PIXEL 





CALL TEP 
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CALL A I 
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CALL IRP 
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AND POLAR CAP 
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SPECTRAL INTENSITIES 
1135.6, LBH-SHORT 
LBH-LONG 
















FULL-SCALE DATA REDUCTION 


THE UVI WILL PRODUCE ONE IMAGE EVERY -30 SECONDS. 

THE NOMINAL OPERATIONAL MODE WILL BE TO CYCLE THROUGH THE 
FILTERS . 


THUS THE MAXIMUM TEMPORAL RESOLUTION POSSIBLE IN THIS MODE IS -4 
MINUTES FOR 8 FILTERS. 

THE FULL-SCA LE DATA REDUCTION WILL INVOLVE THE FOLLOWING STEPS: 

1. GENERATE GLOBAL GRID OF THE ENERGY FLUX USING LBH-LONG . 

2. DECONVOLUTE THE SPECTRAL FEATURES. 

3. COMPUTE SPECTRAL RATIOS NEEDED FOR CHARACTERISTIC 

ENERGIES. 

4. COMPUTE THE GLOBAL GRID OF CHARACTERISTIC ENERGIES. 

1. USING 1135 . 6/LBH-LONG (SEE FIG. 1) 

2. USING LBH-SHORT/ LBH-LONG (SEE FIG. 2) 

5. COMPUTE THE GLOBAL NO COLUMN ABUNDANCE RUN 1215. 

6. COMPUTE THE GEOCORONAL BRIGHTNESS FROM H Ly a. 

7. COMPUTE THE GLOBAL COLUMN 0 ABUNDANCE FROM 1130.4 



SCIENTIFIC MODELING 


THE IMAGES TO BE TAKEN BY THE VISIBLE AND ULTRAVIOLET 
IMAGERS ON THE ISTP MISSION REPRESENT THE END PRODUCT OF MULTI- 
FARIOUS PROCESSES WHICH RESULT IN THE TRANSPORT OF ENERGETIC 
PARTICLES FROM THE SUN TO THE EARTH. THE IMAGES THEREFORE CAN BE 
MAPPED BACK INTO VARIOUS DISTANT REGIONS OF THE SOLAR TERRESTRIAL 
ENVIRONMENT. THUS MEMBERS OF THE ISTP IIWG WILL HAVE THEIR OWN 
PARTICULAR INTEREST IN VARIOUS ASPECTS OF THE IMAGING DATABASE. 

WE CONCENTRATE HERE ON THE SPECIFIC INTERESTS OF THE UVI 


TEAM. 


UVI SCIENTIFIC GOALS: 


PRIMARY GOAL: GLOBAL MODELING OF THE TERRESTRIAL IONOSPHERE, 

AIRGLOW AND AURORA. 


MODEL: FIELD LINE I NTERHEMMI SPHERIC PLASMA MODEL (FLIP) (TO BE 

UPGRADED FOR ISTP) * GLOBAL EMISSIONS MODEL 

MAJOR PREVIOUS LIMITATION: GLOBAL CORPUSCULAR IONIZATION AND 

EXCITATION SOURCES WERE UNKNOWN. 

NOTE: NO IMAGING SYSTEM OF THE PAST HAS PROVIDED 

SUFFICIENT INFORMATION TO GENERATE GIOBAL IONIZATION 
AND EXCITATION RATES. 

THE UVI WILL DO THIS. 

NEEDED DATA: ENERGY FLUX AND CHARACTERISTIC ENERGY PER PIXEL. 

FLIP * GLOBAL ATLTITUDE PROFILES OF 

IONIZATION AND EXTRACTION RATES. 


GLOBEL MODELING 



MID AND LOW LATITUDE MODEL OF THERMOSPHERIC EMISSIONS 
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Auroral Modelling of the 3371 A Emission Rate: Dependence 
on Characteristic Electron Energy 


P. G. Richards 

Computer Science Department and Center for Space Plasma and 
Aeronomic Research., University of Alabama in Huntsville 

D. G. Torr 

Physics Department and Center for Space Plasma and Aeronomic 
Research, Untverstty of Alabama in Huntsville 

We have developed an efficient two-stream auroral electron model to study the deposition of 
auroral energy and the dependence of auroraJ emission rates on characteristic energy. This model 
incorporates the concept of average energy loss to reduce the computation time. Our simple two- 
stream model produces integrated emission rates that are in excellent agreement with the much 
more complex multistream model of S(rtcHan<i fl al. [1983] but disagrees with a recent study 
by Rees and Lummerzheim [1989] that indicates that the N 2 second positive emission rate is a 
strongly decreasing function of the characteristic energy of the precipitating flux. Our calculations 
reveal that a 10 keV electron will undergo approximately 160 ionising collisions with an average 
energy loss per collision of 62 eV before thermalixing The secondary electrons are created with 
an average energy of 42 eV. When all processes including the baclcscatlered escape fluxes are 
taken into account, the average energy loss per electron-ion pair is 35 eV in good agreement with 
laboratory results. 
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Vacuum ultraviolet thin films. 

1: Optical constants of BaF 2 , 
CaF 2 , LaF 3 , MqF 2j Al 2 0 3 , Hf0 2j 
and Si0 2 thin films 


Muamer Zukic, Douglas G. Torr, James F. Spann, 
and Marsha R. Torr 


a reprint from Applied Optics 


a reprint from Applied Optics 


Vacuum ultraviolet thin films. 1 : Optical constants of 
BaF 2 , CaF 2 , LaF 3 , MqF 2 , Al 2 03 , Hf0 2 , and Si0 2 
thin films 


Muamer Zukic, Douglas G. Torr, James F. Spann, and Marsha R. Torr 


The optical constants of MgF2 (bulk) and BaF2, CaF?, LaFi, MgF2, AI2O3, Hf02, and SiO? films deposited on 
MgF 2 substrates are determined from photometric measurements through an iteration process of matching 
calculated and measured values of the reflectance and transmittance in the 120-230-nm vacuum ultraviolet 
wavelength region. The potential use of the listed fluorides and oxides as vacuum ultraviolet coating 
materials is discussed in part 2 of this paper. 


I. Introduction 

The optical constants of materials in the vacuum 
ultraviolet ( VUV) region of the spectrum are of inter- 
est to several areas of technology. Most applications 
such as high reflectivity mirrors and bandpass and 
narrow bandpass interference filters involve thin 
films. The design of optical instrumentation for space 
astronomy, space aeronomy, spectroscopy, and the de- 
velopment of electrooptic devices are critically depen- 
dent on experimental data for optical constants of thin 
films. 

We report measurements of optical constants of 
BaFo, CaF 2 , LaF^, MgF 2 , AI 0 O 3 , Hf0 2 , and Si0 2 thin 
films, and MgF 2 bulk material over the 120-230 nm 
VUV spectral range. The optical properties of MgF 2 
film and bulk materials in the vacuum ultraviolet spec- 
tral range have been reported in a number of papers , 1 " 9 
while optical constants of other fluoride and oxide 
films presented in this paper are reported we believe 
for the first time for the 120-230 nm wavelength range. 
It is found that BaF 2 and LaFi may be used as high 
refractive index film materials for constructing a high- 
low index pair with MgFo being the most useful low 
index material in the VUV. 

Optical properties of lanthanide trifluoride films 
and their potential use as high refractive index film 
materials in the vacuum ultraviolet were investigated 
by Lingg et al. 10 Although optical constants of lantha- 


James Spann and M. R. Torr are with NASA George C. Marshall 
Space Flight Center. Huntsville, Alabama 35812; the other authors 
are with University of Alabama in Huntsville. Physics Department, 
Huntsville, Alabama 35899. 
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nide trifluorides for the VUV range are not presented 
in the paper, plots of the measured transmittance sug- 
gest that some of these trifluorides may be useful coat- 
ing materials for wavelengths as low as 140 nm. Mea- 
sured transmittances of films obtained by the ion- 
assisted and conventional depositions are compared 
and discussed within the paper. 

Optical constants of an isotropic material — refrac- 
tive index n and extinction coefficient k — can be in- 
ferred from photometric or polarimetric measure- 
ments. Since compensators or transmission 
polarizers are not generally available in the VUV from 
120 to 230 nm, photometric methods are almost exclu- 
sively used for obtaining n and k . A number of meth- 
ods exist for extracting optical constants from reflec- 
tance ( R ) and transmittance (T) measurements at both 
normal and oblique angles of incidence. Because R 
and T are very complicated functions of optical con- 
stants 1112 it is generally impossible to express optical 
constants as explicit functions of R and T. Most ap- 
proaches to solving this complicated dependence in- 
volve either graphic or numerical techniques . 13-20 
The expressions are nonlinear so the problem may be 
considered as a numerical exercise in which n and k are 
found through an iteration process of matching calcu- 
lated and measured values of reflectance and transmit- 
tance. The retrieval of optical constants from mea- 
surements of R and T is a well-established approach, 
but it has not been fully exploited for the VUV regime. 
The numerical method employed by us is based on the 
use of a damped least-squares fit approach which is 
incorporated into a thin film design computer pro- 
gram. The damped least-squares fit technique pro- 
vides a rapid and reliable retrieval mechanism for n 
and k in matching measured and calculated quantities. 
It is surprising that to our knowledge the method has 
not been fully utilized previously. This fitting method 
is described in Sec. II. 
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I 

"I"m 

Fig. 1. Measurement of reflectance R\ from the wedged substrate 
to avoid contribution from the back side reflectance. Transmit- 
tance T\ t was measured on a 2-mm thick parallel substrate. 


In Sec. Ill the experimental techniques are nresent- 
ed, including sample preparation and handling, reflec- 
tance and transmittance measurements, and deposi- 
tion of thin films. The optical constants of MgFo 
(bulk) and BaF . JaF 2 » LaF i, MgF*?) A1?0^, HfOo, and 
SiC >2 thin films (deposited on MgF> substrates) are 
given in Sec. IV. Conclusions and a summary are 
given in Sec. V. 

II. Determination of Optical Constants 

Beam diagrams for the measurement and calcula- 
tion of the optical constants of the substrate and sub- 
strate with a six: e thin film are shown in Figs. 1-4. 
Reflectance R i turn the semi-infinite media is mea- 
sured bv means of a wedged substrate, and transmit- 
tance T\t is measured using a plane-parallel substrate 
as shown in Fig. 1. From the beam diagram for theo- 
retical derivation of transmittance of tiu nonabsorb- 
ing slab To, shown in Fig. 2, it follows that 

^ = 7,(1-/?,) + 7^(1 -K.) 

+ T 1 fft(l-« 1 ) + T l fl?(l-R 1 ) + . ... U> 



Fig. 2. Beam diagram for theoretical derivation of the transmit- 
tance of the nonabsorbing thick slab. 



T 



/ * 

7 R£ 

> 

f 

R T «= 


T 

T F 
1 M 


Fig. 3. Measurement of reflectance of the single film deposited 
, he wedged substrate. R F is the calculated reflectance of the 
film with the substrate as an incident medium and air as an 
eme ng medium. 7> is the calculated transmittance of the single 
film sandwiched between two semi-infinite media— air and sub- 
strate. 


and after multiplying we get 
T„ = Td + -«? + «!-«■? + «? + .. . ). 12) 

Using the binomial expression, the transmittance of a 
thick nonabsorbing slab can be written as 


1 +R, “ l + «/ 

where R ] is the measured reflectance from the single 
side of the substrate. If absorption is present in the 
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Fig. 4. Beam diagram for theoretical derivation of transmittance 
T {) of the single film deposited on the nonabsorbing substrate. 

slab, the measured transmittance is smaller than that 
calculated theoretically using Eq. (3). If the measured 
transmittance is denoted as T\f t the total loss due to 
the substrate absorption is given by 



and the correction factor for any other transmittance 
measurement on this substrate at the particular wave- 
length by 



Since the imaginary part of the refractive index of the 
MgF ‘2 substrate is very low (of the order of 10 -7 ), it does 
not affect significantly the measured reflectance R i on 
the wedged substrate. Thus, R { can be used in the 
calculation of refractive index n of the substrate, i.e., 


The ratio I(z)/I( 0) is the so-called intrinsic transmis- 
sion of the medium (no reflection occurs). Using our 
correction factor C defined in Eq. (6) we have 



where D is the thickness of the substrate. Thus, from 
measured values R\ and Tv/ using Eqs. (7) and (11), 
optical constants of the weakly absorbing substrate 
(MgF *2 in our measurements) can be obtained. As an 
example, consider measurements done at X = 135 nm 
on the MgF 2 substrate with thickness D = 2 mm. The 
measured reflectance of the wedged substrate was R\ = 
5.5% and the measured transmittance Tm = 77.4% 
givingn = 1.61 and& = 6.2 X 10“ ' which agrees perfect- 
ly with known values for bulk MgF 2 in the VUV. 7 

From the beam diagrams for reflectance and trans- 
mittance measurements given in Figs. 3 and 4 it follows 
that Tq is 

T* = T f T x [ 1 + R X RS + (fl,/tf) 1 2 + ... 1, (12) 


and again using the binomial expansion we obtain 


t f t } - /?,) 

l -R X RS l — R\R? 


(13) 


where R% is the calculated reflectance of a single film 
with a substrate as an incident medium and air as an 
emerging medium (substrate), and TV is the calculated 
transmittance of a single film sandwiched between the 
incident medium and semi-infinite substrate. 

Equation (13) gives the transmittance of the single 
absorbing film on the nonabsorbing substrate. Using 
the correction factor defined in Eq. (5) we express the 
transmittance of the substrate with a single film 


Tr = C 


T F (l-R ,) 


(14) 


1 - R X RS 

Using Eq. (14) and calculating /?£ from the initial 
values of n and k, we form the merit function F defined 


as 


F = W x iR f - Ri') 2 + WM7t, - 7?r\ 


(15) 


1 - v *i 

The correction factor given in Eq. (5) can be consid- 
ered as the intrinsic transmission of the substrate ma- 
terial. Using the known relation for light propagating 
through the absorbing media 

I{z) = 1(0) exp(-az), (8) 

where the coefficient of absorption a is given by 


where R{ and /?£ are the measured and calculated 
reflectances of the film on the wedged (semi-infinite) 
substrate; Ty and T h c are measured and calculated 
values of transmittance through the plane-parallel 
substrate with a single film; W\ and Wo are the weight- 
ing factors chosen according to the relative accuracy of 
the R and T measurements for each wavelength. 

Merit function F is then minimized using a damped 
least-squares approach. This is implemented as a 
subroutine in a computer program for thin film 
design. 21 Reflectances Rn and R £ and transmittance 
are calculated using values of n and k at the particu- 
lar wavelength. Merit function 



k t 


we get 



F = F[R^R[ y T^,T[(R^)} (16) 

depends on five variables of which three are dependent 
on n and k of the film. Thus, the minimum value of F 
has a high probability of providing accurate values for 
the optical constants of the thin film. 
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To minimize uncertainties in the optical constant 
determination, several single films of the same materi- 
al but with different thicknesses are deposited on sepa- 
rate substrates. The total merit function Ft is then 
given by 

>1 = 1 

where F n is the merit function of the nth film defined in 
Eq. ( 15) and L is the total number of single films with 
different thicknesses deposited either on the same 
type of substrate material or on several different types 
of substrate material. 

III. Experimental Techniques 

A. Sample Preparation and Handling 

All depositions are made on 12.7-mm diam by 2-mm 
thick magnesium fluoride substrates with root mean 
square roughness (usually referred as rms 22 ) <2.5 nm. 
To eliminate the contribution of the back side reflec- 
tion to the reflectance measurements, some of the sub- 
strates have a 3° bias. The substrates were cleaned by 
the supplier (Acton Research Corp., Acton, MA) using 
the following procedure: optical soap wash, water 

rinse, ethanol soak then ultrasonic bath, fresh ethanol 
rinse, and finally a Freon rinse. 

The substrates were shipped in Delrin holders 
wrapped with lens paper and were only removed im- 
mediately prior to deposition. After deposition, the 
substrates were allowed to cool to 40° C, and the vacu- 
um chamber was vented with dry nitrogen. The sub- 
strates were kept in the flow of dry nitrogen and placed 
in a stainless steel container to prevent contamination 
due to exposure. All depositions were made at the 
University of Alabama in Huntsville Optical Aerono- 
my Laboratory and the transmittance and reflectance 
measurements were made at Atomic Physics Branch of 
the Marshall Space Flight Center. 

B. Deposition of Thin Films 

The vacuum system consists of a cryo pump and a 
sorption pump giving an oil-free environment for all 
depositions and therefore providing a very low proba- 
bility for hydrocarbon contamination of the films. 

The film materials BaF 2 , CaF 2 , and LaF 3 were pre- 
pared for vacuum deposition by CERAC with a typical 
purity of 99.9%. AI9O3 (99.5%), Si0 2 (99.98%), and 
MgF 2 (99.95%) are standard Balzers coating materials 
while HfOo with a purity of 99.5% was prepared by EM 
Chemicals. 

The fluoride films were deposited with low imposi- 
tion rates on heated substrates. To reoxidize dissoci- 
ated molecules of the oxides ALO3, Hf0 2 , and Si0 2 , the 
films were deposited at a low deposition rate on heated 
substrates in an oxygen residual atmosphere. The 
temperature of the substrate was monitored with a 
Chromel-Alumel thermocouple attached to the alumi- 
num substrate ring holder. The substrate and its ring 
holder were placed in a 6-mm thick stainless steel plate 


Table I. Deposition Condition* 


Material 

Pa 

\ Torr) 

p 

(Torr) 

d r 

fnm/s) 

d 

(nm) 

BaF* 

9.5 X 10”' 

1.5 X 10"* 

0.16 

53.0 

CaF-j 

8.5 X 10" : 

2.5 X 10~* 

0.20 

99.0 

La Fi 

S.5 X 10" 7 

2.0 X 10“** 

0.14 

51.0 

MgF>> 

3.5 X 10~ 7 

2.5 X 10~ h 

0.23 

68.0 

Al»Oi a 

8.5 X 10’ 7 

l.OXlO- 4 

0.10 

112.0 

HfOj 

S.5X10- 7 

1.5 X 10 -4 

0.10 

30.0 

SiO? 

9.5 X10' 7 

1.2X10" 4 

0.10 

51.0 


Oxygen leaked into the chamber during deposition of the oxides. 


with 40-cm diameter. Further details about the con- 
ditions of depositions are given in Table I. 

The depositions were made with an electron gun. 
The gun had a fixed voltage of 10 kV and low power 
depositions were maintained by supplying low current 
to the gun. The source-to-substrate distance is 50 cm 
and the source-to-thickness monitor distance was 35 
cm. The thickness control and rate measurements 
during film depositions were done with the Kronos 
Digital Film Thickness Monitor QM-300 series and the 
Kronos Deposition Rate/Thickness Output Accessory 
RI-100/RO-200 series with FFT-300 transducer. 

The vacuum chamber geometry provided excellent 
calibration constants C c for all Film materials. The 
calibration constant C c is defined as the ratio of the 
film thickness expressed in transducer counts (hertz) 
and the measured thic. .ness in nanometers, i.e., the 
number of counts (hertz) needed to obtain a 1-nm 
thick film, ihe values of C c varied from 44.09 Hz/nm 
for MgF 2 up to 72.23 Hz/nm for LaF*. The physical 
thickness measurements of ihe films were made with a 
Talystep stylus profilometer. The stylus radius is 2 
/xm and the stylus loading is 1 mg. 

The physical thicknesses for each material are cho- 
sen according to the following criteria: 

(a) nonzero transmittance measurements over most 
of the considered spectrum; 

(b) no overlap of maxima and minima for R and T 
measurements; and 

(c) at least one measurement must be close to the 
possible physical thickness in the design of either a 
multilayer interference filter or an absorbing single 
layer coating. 

C. Photometric Measurements 

Transmittance and reflectance measurements were 
performed in a hydrocarbon-free vacuum system at 
pressures below 10 - ° Torr. A sealed deuterium lamp 
with a MgF 2 window was used in tandem with a 0.2-m 
monochromator producing a beam with 1-nm FWHM 
spectral resolution. Folding and collimating optics 
were used to produce a 1- X 0.75-cm reference beam 
which is incident on an eight-position filter wheel con- 
taining the substrates. Different detectors were used 
for transmittance and reflectance measurements. 
Each detector consists of a sodium salycilate-coated 
Pyrex window placed in front of a bialkali photometer. 
A schematic diagram of the system optics is shown in 
Fig. 5. 
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Deuterium Lamo 


100 



Fig. 5. Schematic diagram of the single beam VUV spectrophotom- 
eter. 


Absolute transmittance and reflectance were mea- 
sured by determining the ratio of the transmitted or 
reflected beam intensity and the unattenuated inci- 
dent beam intensity. For reflectance measurements, 
the unattenuated beam intensity was determined by 
the measured reflectance of a calibrated VUV-en- 
hanced aluminum mirror located in one of the filter 
wheel positions. The reflectance measurements were 
made at a 6° angle of incidence. 

The aluminum mirror was calibrated by Acton using 
the ACE-type 23 self-calibrating VUV photometer. 
The estimated error of this instrument for reflectance 
measurements is <2%. The uncertainties associated 
with the thin film thickness, reflectance, and transmit- 
tance measurements resulted in the total uncertainty 
for optical constants determination of the order of 
±5%. This uncertainty is derived from discrepancies 
between the theoretically and experimentally ob- 
tained spectral performance of deposited VUV multi- 
layer filters. In the theoretical calculation, films are 
assumed to be homogeneous, isotropic, and bounded 
by two infinite ideal planes. Thus, the discrepancy 
between theory and experiment is caused by both 
physical effects (such as film inhomogeneity, surface 
and volume scattering, diffusion, and possible con- 
tamination) neglected in the theoretical calculation 
and the uncertainty in the optical constant determina- 
tion. Taking all this into account we might say that 
the maximum uncertainty for optical constant deter- 
mination is <5%. 

IV. Optical Constants 

A. MgF 2 Substrate 

The reflectance and transmittance measurements 
shown in Fig. 6 were made on wedged and 2-mm thick 
parallel substrates, respectively. Optical constants n 
and k shown in Fig. 7 were determined using Eqs. (8) 
and (12). 

B. Fluoride Films 

The temperature of the MgF 2 substrate during de- 
position of the fluoride films is 250°C. The pressure 
before (P 0 ) and during deposition (P) as well as the 
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Fig. 6. Reflectance R\ and transmittance Tv/ measurements on the 
MgFo wedged and a 2-mm thick parallel substrate, respectively. 
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Fig. 7. Optical constants of the MgF> substrate determined using 
Eqs. (81 and ( 12). 


deposition rates (Dr) and physical thicknesses of the 
films ( d ) are listed in Table I. Single sets of the R and 
T measurements of BaF 2 , CaF >, LaF and MgF 2 films 
deposited on MgF? substrates are given in Figs. 8, 10, 

12, and 14, respectively. The corresponding optical 
constants determined from at least two independents 
and T measurements [Eq. (19)] are given in Figs. 9, 11, 

13, and 15. 

From the reflectance and transmittance curves ot a 
53-nm single film of BaF 2 , it follows that this coating 
material can be used for wavelengths longer than 135 
nm. Even a single film of BaF 2 could be used as a cut- 
on filter if wavelengths shorter than 130 nm are not 
desired. Refractive index n has values between 1.87 
and 2 for wavelengths from 125 to 135 nm and it is 
higher than 1.7 throughout the region from 140 to 210 
nm, increasing up to 1.98 at 230 nm. Extinction coeffi- 
cient k has values of the order of 10" 2 for wavelengths 
from 140 to 230 nm. 

The measured reflectance and transmittance of a 99- 
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Fig. 8. Reflectance /?f and transmittance Ty of a 53-nm thick BaF? 
film deposited on the MgF> substrate. 
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Fig. 10. Reflectance /?f and transmittance of a 99-nm thick 
CaF 2 film deposited on the MgF 2 substrate. 
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Fin 9. Optical constants of BaF > determined from R and T mea- 
surements of 53-, 45-, and 94.5-nm thick films deposited on the MgF > 
substrate. 
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Fig. 11. Optical constants of CaF> determined from R and T mea- 
surements of 99- and 54.5-nm thick films deposited on the MgF 2 
substrate. 


nm thick CaF 2 film indicate low values of refractive 
index n and a relatively small extinction coefficient k 
for wavelengths longer than 135 nm. The values of 
refractive index n of CaF 2 are lower than 1.3 for X ^ 180 
nm, and the h values are of the order of 10~ 2 for the 
same wavelength region. Low values of the refractive 
index make CaF 2 suitable for use as the alternative low 
index material for the longer wavelengths of the V UV. 

The material with the highest values of refractive 
index among all the fluorides n ^ 1.85 within the 135— 
230 nm range of the VUV wavelengths is LaF 2 . The 
extinction coefficient has values lower than 2.2 X 10 
for X > 145 nm. Compared with other fluoride and 
oxide coating materials LaF d seems to be the best 
choice for the high index material in the considered 
spectral region. 

Magnesium fluoride films have extinction coeffi- 
cient values lower than 10~ 4 throughout the entire 
140-230-nm wavelength region. This makes MgF 2 the 


most attractive low index material for the VUV. The 
n and k values of MgF 2 film deposited on the substrate 
heated up to 250°C reported by Wood et al. 1 at 121.6 
nm are n = 1.7 ± 0.1 and k close to 0.01 are in good 
agreement with values reported here. Unfortunately, 
we are unabL to compare these n and k values of MgF 2 
with the values obtained by some other authors be- 
cause either the MgF 2 films were prepared in different 
conditions or insufficient experimental data were pro- 
vided by the ether workers. 1 ' 10 

C. Oxide Films 

The temperature of the MgF 2 substrate during de- 
position of oxide films was 300°C. The deposition 
conditions Po, P, and Dr, as well as the physical thick- 
nesses of films d, are listed in Table I. The measured 
spectral curves of R and T for A1 2 0 ;, Hf0 2 , and Si0 2 
films deposited on MgF 2 substrates are shown in Figs. 
16, 18, and 20 while the corresponding optical con- 
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stants are given in Figs. 17. 19, and 21, respectively. 

The transmittance curve of a 112-nm thick ALO) 
single film suggests that this material could possibly be 
used for the design of absorption edge filters with cut- 
on wavelengths between 160 and 180 nm depending on 
the thickness of the AI 2 O 2 film. The values of refrac- 
tive index n are higher than 1.85 for almost the entire 
region from 120 to 230 nm while the k values are higher 
than 10 -1 for 120 nm < \ < 175 nm. 

The HfO > film has the highest extinction coefficient 
(shown in Fig. 19) and it does not seem that this mate- 
rial could be useful for the vacuum ultraviolet wave- 
length region, particularly for wavelengths below 200 
nm. 

The values of refractive index n of SiO? films are 
very close to the values of bulk SiCT, while the k values 
are an order of magnitude higher. 24 The refractive 
index n > 1.8 for 125 nm < \ < 150 nm and decreases 
gradually to 1.6 at 200 nm. Extinction coefficient k is 
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of the order of 10“ 1 for 120 nm > X > 135 nm which 
makes SiC >2 a useful material for the design of a single 
layer absorption edge filter with cut-on wavelengths 
within the interval from 125 to 135 nm depending on 
the relative thickness of SiO? film. 

V. Summary 

We demonstrated that iterative mathematical mod- 
eling of transmittances and reflectances measure- 
ments provides a reliable way for determining the opti- 
cal constants of thin films deposited on the weakly 
absorbing substrates. In part 2, proof of the validity of 
the approach is demonstrated by design and fabrica- 
tion of multilayer coatings. The model is used for the 
VUV wavelength region but its application extends 
over the whole visible and IR spectrum whenever a 
substrate has a low value of extinction coefficient k 
such that its effects on the substrate reflectance are 
negligible. 
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The fluoride film materials have generally lower 
extinction coefficients than the oxides. High values of 
refractive index n make LaF^ and BaFo useful materi- 
als for the VUV, particularly for constructing a high- 
low index pair with MgF2 being the most useful low 
index material. The designs of the VUV coatings such 
as narrowband pass, and narrowband reflection filters 
are possible with these materials. 

SiOo coating material, among other oxides measured 
in the VUV, seems to be the most applicable for design 
of multilayer stacks such as the narrowband reflection 
filters. AI2O3 coating material may be used for the 
design of absorption edge filters for A > 160 nm, while 
Hf02 becomes a useful high index material for wave- 
lengths longer than 230 nm. 
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Vacuum ultraviolet thin films. 2: Vacuum ultraviolet 
all-dielectric narrowband filters 


Muamer Zukic, Douglas G. Torr, James F. Spann, and Marsha R. Torr 


We report the design and performance of narrowband transmission filters employing the rapidly changing 
extinction coefficient that is characteristic of BaF? and Si02 films within certain wavelength intervals in the 
vacuum ultraviolet. We demonstrate the design concept for two filters centered at 135 nm for BaF 2 and at 141 
nm for SiOj. It is found that these filters provide excellent narrowband spectral performance when combined 
with narrowband reflection filters. The filter centered at 135 nm has a peak transmittance of 24% and a 
bandwidth of 4 nm at full width at half-maximum for collimated incident light. The transmittance for Ao < 
130 nm is <0.1% and for 138 < An < 230 nm the average transmittance is <3%. Another filter centered at 141 
nm has a peak transmittance of 25% and a bandwidth of 3.5 nm. 


I. Introduction 

b The design of all-dielectric multilayer interference 
filters for the vacuum ultraviolet (VUV) wavelength 
region from 120 to 230 nm is limited by the lack of film 
materials with suitable optical constants. Film mate- 
rials such as MgF 2 and LiF have low values of refrac- 
tive index n (the real part of the optical constant) and 
relatively low values of extinction coefficient k (the 
imaginary part of the optical constant) within the 
VUV wavelength region. The low resistance to high 
energy radiation damage in a space environment 1 
makes LiF not very useful as a material for interfer- 
ence optical filters intended for many space applica- 
tions, reducing the choice of low index film materials to 
just MgF>>. 

For their values of the extinction coefficient, both 
BaF ^ and LaF 2 may be used as the high index materi- 
als" to form a high-low (HL) index pair together with 
MgF 2 . A HL pair, made either with BaF 2 -MgF 2 or 
LaF }-MgF 2 , can provide the basic sequence of a multi- 
layer design which can be used for wavelengths as low 
as 130 nm. A brief review of the basic theory for such 
multilayers is given in Sec. II, while theoretical calcula- 
tions of so-called tuned multilayer interference filters 
with absorbing film materials are given in Sec. III. 

All-dielectric Fabry-Perot type filters employing 
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BaF 9 and MgF 2 are considered in Sec. IV along with an 
alternative design for narrowband transmission filters. 
Designs of narrowband reflection filters are also given 
in that section. Spectral performance is compared to 
transmission filters. It is found that, for a small cone 
of light (±5°) centered about some incident angle, 
reflection filters when combined with transmission fil- 
ters can provide excellent spectral performance in the 
VUV. A summary and conclusions are given in Sec. V. 

II. Basic Theory 

A. VUV Absorption of Dielectrics 

Absorption of dielectrics in the VUV can be treated 
with classical or quantum mechanical theory. Be- 
cause of the present accuracy of the reflectance and 
transmittance measurements there is no advantage to 
using the quantum mechanical approach; thus the 
classical model is used for the calculation of optical 
constants. In this model, bound electrons in a dielec- 
tric illumination with electromagnetic radiation are 
treated as damped harmonic oscillators. 3 

To determine the optical constants of metals, the 
classical theory uses the free electron gas model. 4 Ab- 
sorption of incident electromagnetic radiation in di- 
electrics and in metals are two different physical phe- 
nomena (bound electrons in dielectrics and free 
electrons in metals), but the final result is the same, 
i.e., loss of incident intensity. 

Both phenomena can be described in terms of a 
complex optical constant. However, the transition 
from the real optical constant of a lossless medium 
(dielectrics in the visible part of the spectrum) to the 
complex one of an absorbing medium requires a redefi- 
nition of phase velocity v and wavenumber k\ such that 
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( where N is the complex optical constant) are now the 
complex quantities. In addition, the angles between 
the direction of the light propagation and the normal 
to the film plane are also the complex quantities and 
they no longer represent just the refraction of the 
propagating light. 

B. Multilayers with Absorbing Film Materials 

The VUV coatings presented in this paper consist of 
an HL pair with BaF- 2 , LaF^ or SiO ? as the high index 
materials and MgF? as the low index material. Since 
extinction coefficient k of the MgFs films is <10“ 3 for 
wavelengths longer than 130 nm, it is neglected in the 
theoretical discussion concerning the angles of light 
propagation through the absorbing multilayers. 
However, the extinction coefficient of MgFo films is 
taken into account in the exact calculation of multi- 
layers. 

The amplitude reflection and transmission coeffi- 
cients for the plane electromagnetic wave incident on a 
multilayer are 3 

_ ~ + A 4y?} s ) ^ 

r= (M u + M yl ,) % + (M- n + M, u ) 


{M u + + Af?>v? 9 ) 

where 770 and v, are the effective optical constants of 
the incident medium and the substrate, respectively. 
They are defined as 

17,1 “ COSt^n. 

77 n = rj s COS#, 

for s-polarization; and for p- polarization 


The phase terms di , which are usually reierred to as the 
phase thicknesses of the films, are defined by 

0,. a — COSOn 

' ' ' 

where An is the vacuum wavelength of the incident 
light, Ni is the optical constant of the Ith layer defined 

as 

.V = n 1 + Ur)= n : + in,K l = n, + ik ; , (12) 

with K[ — htjnu and where fii is the refractive index, k / is 
the extinction coefficient. d ( is the physical thicknes- 
s.and 0/ is the complex angle of light within the Ith 
film. 

If the complex angle 0/ is written as 
sin0/ = Qi expU7j)» 

the generalized Snell's law applied at the interface 
between the nonabsorbing (i “ l)th medium and ab- 
sorbing Ith medium is given by 

•V, sine, * (n ! -I- ik,)q, expUT/) = n sin0,_[. (14) 

The right-hand side of Eq. (14) is real, thus 

n>qi cosy , - k,q ( sini^ = n -i sinO^,. (15) 

n,q t sin 7 / + k ; q t cosif = 0 

From Eqs. (15) and (16) it follows that q t and y / are 
given by 

n /-i sin0i_i q7) 

~ --- — -> • 

\ n l + 


7 / « COS 


v n/ + k\ 


If ki = 0, then 

n,_, sin0/_. 


where and n* are the refractive indices of the inci- 
dent medium and the substrate, respectively. It is 
assumed that both the substrate and incident medium 
have negligible extinction coefficients and. therefore, 
real optical constants. The terms A/ rM ij =1,2 are the 
elements of the multilayer characteristic matrix M 
which is defined as the product of the matrices of the 
individual layers M/ (/ = 1,2,. . .P), 


where P is the total number of layers. Matrices M; are 
given by 


it}, sin<), eos<5, 


and the generalized Snell’s law [Eq. (14)] should re- 
duce to the well-known law of refraction. 

The reflection and transmission coefficients are 
complex numbers of the form 

r = lr| expU0j. 

t = |fl expU0j, 

where <p r and </> t represent the phase changes on reflec- 
tion and transmission, respectively. The phase 
change of reflected light </> r is referred to the plane 
bou dary between the semi-infinite incident medium 
ana he front surface of the multilayer, while the phase 
ch; of transmitted light <$>> is referred to the plane 
bou- xry between the multilayer and the semi-infi- 
nite medium of the substrate. Reflectance R ♦ trans- 
mittance r, and absorptance A of the multilayers are 
given by 

R = rr*, < 23 > 
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r = 

.4 = 1 - + D. 


(24) 

(25) 


111. Tuned Multilayers 

Multilayers formed by high and low index materials 
alternating throughout the stack are. by analogy with 
electrical networks, called tuned filters or tuned multi- 
layers. Because of their importance in the design of 
VUV coatings presented in this paper, we investigate 
the properties of tuned multilayers with absorbing 
film materials in more detail here. 

Consider the multilayer (HL)^ with a total number 
of films P = 2p. Let the angle of incidence be zero, i.e., 
0o = 0. By denoting the angle within the low index 
material as 0 l and the angle within the high index 
material as 0 h we may write 


(26) 


Thus, phase thickness 6 h and <5 l are now given by 


“ T~ in H + t^H^H 

'Ml 


(27) 


for the high index material, and for the low index 
material 


. 2 7T , .. , 

5 L = — (n L + i* L )d L . 
An 


(28) 


If the optical thicknesses of both the H and L materi- 
als are quarterwave relative to the same reference 
wavelength A r , then 

n iA = T - 


n u dyj = — . 


Phase thicknesses <$h and are given by 


- /. , M 

0.1 = 1 + l I , 

- X 0 \ n H / 

X ^ 1 ‘ ( i u_ 

' 2A,\ n L ) 


Further, with 




Eqs. (31) and (32) can be written as 

T 'V 

- A 0 
7T X ' 

= 77“ <1 + lK 0- 

4. Afl 


(29) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


The matrices of such quarterwave films at A r = Ao 
become 


r.^G€ IS 
OF POOR QUAUTY 


M-- I 


M, =i 


I -i sinho H — coshcru \ 

•»H 

iiV H cosha H -i sinha H 


( — i sinhai — cosha L \ 

tiVj coshaj — r sinha L J 


where 


7T 7T 

T * 


2 n 


(37) 

(38) 

(39) 

(40) 


L 


At this point, we shall assume (just for theoretical 
consideration of the tuned multilayers with absorbing 
films and not for exact calculation) that the hyperbolic 
functions in Eqs. (37) and (38) may be approximated 
by their values at the origin, i.e., sinhx — 0 and cosh* 
— 1 for x 0. Matrices Mh and Ml within this 
approximation and for Ao — A r can be written as 


Mu = i 




iNu 



(41) 


(42) 


giving matrix Mi of the basic sequence (HL): 

o 


0 


M t = M h M l = | 


(43) 


The characteristic matrix of multilayer (HL) P is 
then given by 


Mp 


mM„M l ) p = 




ar. 


(44) 


From Eq. (3) the amplitude reflection coefficient fol- 
lows as 


no\Nj 
n o \N l / 


(45) 


The ratio of optical constants N h and N l can be 
written as 


,V H n H (l+u H ) n H 

— — (a - + b~) ~ exp(ta), 

&L ”l( 1 + 1 *l> "i. 


where 


1 + 


H a L 


b = 


1 + Kf 
*H ~~ «L 

i + *r 


(46) 


(47) 


(48) 
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.j = ran 




Thus. l\'hAV| )p can be written as 

/JV H V 

= Ffp) exp[iV(pl|. 

where 

'^YV 


Fip) 




2 (pi = pj. 

Now, Eq. (45) for r becomes 


l 


rip) = 


-Vi 

” 0 / 


(2p) exp(iV(2p)| 


1 + 




<49> 
1 50) 

(51) 

(52) 

(53) 


'(2 p) exp(t^(2p)j 


giving for reflectance Rip) of the (HL) p multilayer 
tuned at wavelength An = A r : 


Rip I 


i + ( — | rV2p) - i( -\ f(2 

\n 0 f 


F(2p) cos\P2p) 


1 + - 


*n/ 


(54) 


(2p) -I- 2| — )F(2p) cos\£t2p) 


If p = 0, i.e., if no films are present on the substrate, 
Eq. (54) reduces to the well-known Fresnel formula for 
reflectance at the boundary between two semi-infinite 
media with indices no and n v : 


fl(0> = 


'nr 


(55) 


From Eq. (54) it follows that the reflectance of the 
(HLV 1 stack will have maximum values if, for l = 
0 , 1 , 2 . . .. 


;c2p) = ' i + 1 ■' 1 ) t , 


(56) 


Thus, for / = 0 from Eqs. (47M49), (51), (52), and (56) 
it follows that the maximum reflectance is achieved if 


Equation (54) indicates that R(po) = 1 if p satisfies 
Eq. (57). This is caused by the high level of approxi- 
mation at the beginning of our derivation. Equation 
(57) is the final goal of this theoretical consideration 
for determining the maximum useful number of layers 
for achieving the highest possible reflectance with ab- 
soroing materials. 

From a certain number of (HL) pairs po, both the 
reflectance and absorptance of the stack will remain 
practically constant, adding to unity, and no signifi- 
cant improvement of the stack’s reflectance can be 
achieved by adding new (HL) pairs, i.e., for p > po: 

Rip) + Aip) — * 1. (58) 

The transmittance of the stack will be equal to zero for 
all practical purposes. 



The results of exact calculation of how reflectance 
R(p) and absorptance Aip) depend on the number of 
(HL) pairs p of multilayers formed with BaF 2 , LaF^, 
and SiO 2 as the high index materials and MgF 2 as the 
low index material are shown in Figs. L 2, and 3, 
respectively. The optical film thicknesses of all the 
films are quarterwave relative to wavelength Ao = A r = 
135 nm. At this wavelength the ratio k for the Si0 2 
film is kh = 0.05759, and for the MgF 2 *l = 0.00025, 
giving po = 13.7, i.e., po — 14, which agrees with the 
exact calculation shown in Fig. 1. Since the extinction 
coefficients of BaF 2 and LaF^ are an order of magni- 
tude lower than the extinction coefficient of Si0 2 ,Po 
attains higher values if these fluorides are used as the 
high index material. 

From the reflectance Rip) and absorptance A(p) 
curves shown in Figs. 1-3, it follows that BaF 2 has the 
highest potential as the high index material for wave- 
lengths close to 135 nm. At longer wavelengths, Ao > 
150 nm, the extinction coefficients of LaFj and BaF 2 
become similar while the refractive index of LaF .3 re- 
mains higher.- 

The Rip) and Aip) curves shown in Figs. 4-6 are 
calculated again for multilayers with BaF 2 , LaFi, and 
Si0 2 as the high index materials but with optical thick- 
nesses of 

n H d H = j , (59) 

Ml - 7 . < 60 > 

where is the refractive index of MgF 2 , and Ao - A r = 
135 nm. This type of tuned multilayer stack is some- 
times referred to as the third wave (TW) design, by 
analogy with quarterwave (QW) stacks. The thick- 
ness of the (absorbing) high index material is larger in 
the QW than in the TW multilayers, thus one may 
expect that the TW designs will have lower absorp- 
tance and therefore higher reflectance. That this ex- 
pectation is justified follows from comparison of the 
Rip) and Aip) curves for the QW multilayers shown in 
Figs. 1-3 and Rip) and Aip) curves for TW multilayers 
shown in Figs. 4-6. 

The maximum reflectance of the QW stack with 
BaF? is R = 85.4% (Fig. 1), while the TW multilayer 
with the same high index material has a maximum 
reflectance of R = 87.6% (Fig. 4). The largest relative 
difference between the maximum reflectances of :he 
TW and QW stacks is for a multilayer with Si0 2 as the 
high index material. The maximum reflectance of the 
QW stack is R = 60.9% (Fig. 3), while for the TW stack 
R = 66.3% (Fig. 6). 

The number po found for the quarterwave multi- 
layers using Eq. (57) seems to be useful even for the 
third wave stacks. It again represents the minimum 
number of (HL) pairs po needed to obtain the maxi- 
mum value of reflectance or zero transmittance. 

A more rigorous treatment of absorbing multilayers 
with the so-called potential transmittance and absorp- 
tance is given by KnittL and the original references 
listed therein. 
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Fig. 1. Reflectance and absorptance of the QYV stack as functions of 
the number of (HL) pairs: H = barium fluoride and L = magnesium 
fluoride. 
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Fig. 4. Reflectance and absorptance of the TVV stack as functions of 
the number of { HL) pairs: H = barium fluoride and L = magnesium 
fluoride. 



Fig. 2. Reflectance and absorptance of the QW stack as functions of Fig. 5. Reflectance and absorptance of the TVV stack as functions of 

the number of (HL) pairs: H = lanthanum fluoride and L = magne- the number of (HL) pairs: H = lanthanum fluoride and L = magne- 
sium fluoride. sium fluoride. 
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Fig. 3. Reflectance and absorptance of the QW stack as functions of Fig. 6. Reflectance and absorptance of the TVV 7 stack as functions of 

the number of ( HL) pairs: H = silicon dioxide and L = magnesium the number of I HL) pairs: H = silicon dioxide and L = magnesium 

fluoride. fluoride. 
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IV. Narrowband Filters 

A. Transmission Filters 


L Fabry -Perot Filters 

Transmission or bandpass filters are frequently de- 
signed using the basic structure of the Fabrv-Perot 
(F.P.) interferometer — the muitnayer in which one of 
the layers (spacing layer or spacer) is bounded by two 
partial reflectors. The transmittance of such an inter- 
ference filter is given by 6,7 


1 + F sirr'i' 


where 


T = ■ 


F - ■ 


r,7\ 

u- v 7 p;> 2 

^R\R, 


(1 - x 


* = h - 


01 + 0 .. 


<62) 

(63) 

(64) 


R[ and iZ 2 are the reflectances and T\ and 7T are the 
transmittances of the partial reflectors bounding the 
spacing layer, calculated with the spacer taken as the 
incident medium (seen from inside the spacer). The 
emerging medium for reflector R \ is air and for partial 
reflector Ro it is the substrate. 5 is the phase thickness 
of the spacing layer defined in Eq. (11), <p i and <t > 2 are 
the phase changes [Eq. (17)] associated with reflec- 
tions on partial reflectors R\ and R >. The maxima of 
transmission according to Eq. (61) are given by 

^ — h — 0 = m 7t , (65) 


where 

0 . + 0 ., 


( 66 ) 


and the bandwidth (A X)h full width at half-maximum 
(FWHM) is given by 7 


( AA)j, 



m 


d\ n 


1 I 



(67) 


of the mth order filter (m = 1,2,. . 

For BaFo, LaFi, and Si0 2 the values of the extinction 
coefficients are too high in the VUV- to be used as film 
materials for the spacing layer. Thus, in all our de- 
signs MgFj was used as the film material of the spacer. 
The possible designs of F.P. type interference filters 
are then given by 

air|H(LH) p 2L(HL) p H|substrate, (68) 


or 

air((LH) p 2L(HL) p )substrate. (69) 

where L denotes the quarterwave optical thickness of 
the MgF 2 , H is the quarterwave optical thickness of 
one of the high index materials (BaF2, LaF), or SiO;), 
and the substrate is bulk MgF 2 . The optical thick- 


nesses are quarterwave relative to the central pass 
wavelength of the F.P. filter. 

Reflectances R and R > and transmittances T i and T> 
appearing in the above basic equations for the F.P. 
filter are seen from inside the spacing layer. The R i 
and T\ are calculated for the boundary between the 
spacing layer and one side of the F.P. filter including 
the incident medium (air), while R- 2 and T> are calcu- 
lated for the boundary between the spacing layer and 
the other side of the F.P. filter including the substrate. 
The partial reflectors in both designs [(68) and (69)] 
are (HL) P QW tuned multilayers with an additional 
outer H layer in design (69). 

The dependences of reflectances R i and Ro on the 
number of (HL) pairs p for BaF 2 , LaFn, and Si0 2 are 
shown in Figs. 7, 8, and 9, respectively. The maximum 
values of transmittance T ma x calculated using Eq. (62) 
and the corresponding bandwidths (AX) b calculated 
according to Eq. (67) for the first-order filter (m = 1) 
are plotted vs p in Figs. 10-12 (for the high index 
materials listed above). The differentiation of the 
function 


^(Vl) — A-,0(Ao) 


70) 


in Eq. (67) was done numerically according to 


<I>(X ‘' ) ] V " A = + [ X ° 
r d 1 r*< x o 


d\, 

+ AX) - 1M A 0 ) 
AX 


0(X O ) 




(71) 

172) 


with AA < X 0 X 10" 6 and the central pass wavelength of 
the F.P. filter Ao = 135 nm. 

From Fig. 10 it follows that the F.P. filter centered at 
Aq = 135 nm should have T max — 30% and (AAo)fc ^ 10 
nm for H = BaF 2 and p = 6. The design of the F.P. 
filter becomes 


air((LH) 6 2L(HL) H ]substrate, I' 3 ) 

where the substrate is bulk MgF 2 . 

The experimentally obtained and theoretically cal- 
culated spectral performance of such a filter is shown 
in Fig. 13. The theoretical curve is not corrected for 
loss due to MgF 2 substrate absorption and back side 
reflection (23% at A 0 = 135 nm). The experimentally 
obtained bandwidth is smaller than the one predicted 
in Fig. 10. This can be explained by much higher 
values of the extinction coefficient of BaF‘ 2 for wave- 
lengths below 131 nm than tor Ao ^ 135 nm. 

To improve the transmittance at central wavelength 
Ao = 135 nm, a F.P. filter with two spacing layers may 
be designed such as 

air((LH) 4 2L( HL) 5 H2L( HL) -4 [substrate, (74) 

where H and L represent the quarterwave optical 
thicknesses of BaF 2 and MgF.3 respectively, and the 
substrate is bulk MgF 2 . Th- calculated and experi- 
mentally obtained transmittances of this filter are 
shown in Fig. 14. The effects of the MgF 2 substrate 
absorption and back side reflection are not taken into 
account in the theoretical curve. The bandwidth of 
the filter is again smaller than that calculated using a 
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NUMBER OF (HU PAIRS Fig. 10. Maximum transmittance and bandwidth of the Fabry- 

Fig. 7. Reflectances of the partial reflectors of the Fabrv-Perot Perot type filter calculated using Eqs. (62) and (6 /)♦ respectively. H 
type filter as functions of the number of (HL) pairs: H = barium = barium fluoride and L - magnesium fluoride, 

fluoride and L = magnesium fluoride. 
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Fig. 8. Reflectances of the partial reflectors of the Fabrv-Perot Fig. 11. Maximum transmittance and bandwidth of the Fabry- 

type filter as functions of the number of ( HL) pairs: H = lanthanum Perot type filter calculated using Eqs. (62) and (67), respectively: H 

fluoride and L = magnesium fluoride. = lanthanum fluoride and L = magnesium fluoride. 
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Fig. 9. Reflectances of the partial reflectors of the Fabrv-Perot Fig. 12. Maximum transmittance and bandwidth of the Fabry- 

tvpe filter as functions of the number of (HL) pairs: H = silicon Perot type filter calculated using Eqs. (62) and (67), respectively: H 

dioxide and L = magnesium fluoride. = silicon dioxide and L - magnesium fluoride. 
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Fig. 13. Twenty-five-layer Fabry-Perot type filter: H = barium 
fluoride and L * magnesium fluoride. 
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Fig. 14. Twenty-nine-layer Fabry-Perot filter with two spacing 
layers: H = barium fluoride and L = magnesium tluoride. 

formula similar to Eq. (67) for the bandwidth of the 
F.P. filter with two spacing layers. 0 This can also be 
attributed to the much higher k values of BaFj films 
for shorter wavelengths; k = 0.1 at X = 130 nm com- 
pared with k - 0.026 at A = 135 nm. 

2. Tuned Stack 

The fact that the extinction coefficient of a BaF> 
film is almost four times larger at 130 nm than at 135 
nm can be used for the design of a narrowband filter 
centered at 135 nm. The design of such a filter is a 
simple QW tuned stack with the high reflection zone 
centered at 140 nm. The theoretical and experimental 
spectral curves of the twenty-five-layer QW tuned fil- 
ter are shown in Fig. 15. The interference effects are 
predominant in the wavelength region above 135 nm, 
while the absorption of BaF? dominates for wave- 
lengths shorter than 135 nm. The peak value of the 
transmittance is T max = 39% at A 0 = 135 nm. The 
theoretical curve in Fig. 15 is not corrected for loss due 
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Fig. 15. Twenty-five-laver QW tuned filter with the high reflection 
zone centered at 140 nm: H = barium fluoride and L = magnesium 
fluoride. 
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Fig. 16. Twenty-five-laver TW tuneu filter with the high reflection 
zone centered at 147.5 nm: H = silicon dioxide and L = magnesium 
fluoride. 

to MgFo substrate absorption and back side reflection. 

The filter shown in Fig. 15, which combines absorp- 
tion effects of the film material ( BaF >) to reject shorter 
wavelengths and interference effects to reject longer 
wavelengths relative to the central wavelength Ao = 
135 nm, has a higher peak transmittance and at the 
same time provides better rejection of the longer wave- 
lengths. Another possible design of such a filter cen- 
tered at 141 nm with a twenty-five-layer TW tuned 
multilayer is shown in Fig. 16. The rapidly changing 
extinction coefficient k of SiO> film within 135 nm < A 
< 145 nm is utilized. 

All these filters suffer from pass windows at longer 
wavelengths. An edge filter is needed which will reject 
longer wavelengths. When the edge filter is combined 
with the narrowband transmission filter, the combina- 
tion should provide useful transmittance at Ao. For 
the filters centered at 135 nm, the edge filter might be 
required (for some applications) to reject longer wave- 
lengths up to at least 170 nm to better than 95%, and at 
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Fig. 17. Twenty- five-layer QVV tuned filter with the high reflection 
zone centered at 160 nm: H - lanthanum iluoride and L = magne- 
sium fluoride. 

the same time provide T > 50% at Ao = 135 nm. The 
basic design of such a filter is again a tuned multilayer 
with either QW or TW optical thicknesses. The QW 
stacks generally have a wider high reflection zone than 
TW stacks. The width of the high reflection zone of a 
QW multilayer with nonabsorbing film materials is 
given by 8 



where nn and hl are the refractive indices of the high 
and low index materials, respectively. If the presence 
of absorption in the films affects only the maximum 
value of the reflectance but not the width of the high 
reflection zone, (AA)h.r. = 21.6 nm for H = LaF; and L 
= MgFo at Ao = 160 nm. Because of the transmission 
requirements at 135 nm, the idea of coupling QW or 
TW multilayers centered at several wavelengths with- 
in the wavelength interval from 140 to 170 nm cannot 
be used. Hence, the edge filter which would provide 
the required transmission at 135-nm wavelength and 
at the same time reject longer wavelengths up to 170 
nm cannot be designed with the dielectrics known to 
us. The twenty-five layer QW stack shown in Fig. 17 
justifies calculations done using Eq. (75); the high 
reflection zone of the measured spectral curve ( AA)h.r. 
< 15 nm. 

B. Reflection Filters 

Designs of the dielectric cut-on filters which would 
have a useful range of transmittance (T > 50%) for 
wavelengths below 145 nm and at the same time reject 
the longer wavelengths up to 230 nm better than 95% 
do not seem to be feasible at the moment. Another 
possibility in solving this problem is to try to design a 
reflection filter at an incident angle of, say, 45°, which 
would reject a narrow spectral band with a reflectance 
of 50% or more at the desired wavelength. If this type 
of filter is then combined with one of the previously 
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Fig. IS. Twenty-five4ayer second-order QW tuned filter with the 
high reflection zone centered at 135 nm. The angle of incidence is 
45°: H = lanthanum fluoride and L = magnesium fluoride. 


presented transmission filters, the combination could 
provide excellent rejection for the shorter wave- 
lengths, high peak transmittance, and reasonably good 
blocking of the longer wavelengths. 

The design of the reflection filter is once again a QE 
or TW multilayer with the optical thicknesses of the 
films corrected for oblique incidence. A first-order 
QW stack generally has a wider high reflection zone 
and suffers from high side reflection ripples. The 
ripples can be reduced by introducing films with opti- 
cal thicknesses H72 or L72 (primes denote correction 
of the quarterwave optical thickness for the oblique 
incidence) at the first and last positions in the stack. 
The width of the high reflection zone decreases if the 
order of the QW multilayer is increased. The increase 
of the order by 1 means a change in the optical thick- 
ness of one of the film materials from Ao/4 to 3 Aq/ 4. 
Obviously, in the VUV range MgF 2 is a material whose 
optical thickness can be increased from L' to 3L' with- 
out affecting the maximum value of reflectance possi- 
ble with a first-order QW stack. 

The experimental and theoretical spectral curves of 
the second-order QW stack at an angle of incidence 0 O 
= 45° are shown in Fig. 18. The design of the filter is 
given by 

airj^L i H'3L') n H' ^substrate. (76) 

where H' = BaF 2 L' = MgFj, A 0 = 135 nm, and primes 
denote the correction for oblique incidence. The 
width of the high reflection zone (A A) h r. ^ 5 nm and 
the reflectance at the central wavelength is 60%. 

The measured overall transmittance of the filter 
presented in Fig. 15 combined with the reflection filter 
(Fig. 18) is shown in Fig. 19. The central wavelength 
of this filter is An = 135 nm, the peak value of the 
transmittance is T max = 24%, and the bandwidth (AA)^ 
= 4 nm. The average transmittance for the longer 
wavelengths is <3%, while for wavelengths shorter 
than 130 nm the transmittance is <0.1%. A cone angle 
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Fig. 19. ^asured transmittance of the filters shown in Figs. 15 and 

18 combined. 

of incident light of ±5° causes the bandwidth to in- 
crease up to 7 nm. 

The overall transmittance of the filter shown in Fig. 
16 combined with the QW multilayer reflector cen- 
tered at 1 1 nm is shown in Fig. 20. The bandwidth of 
the filter is 3.5 nm and peak transmittance is 25%. 
The film materials used for the reflection filter are 
LaFs and MgF2. 

V. Summary 

The idea of utilizing the natural absorption of one of 
the film materials to limit the transmission at shorter 
wavelengths and a combination of this filter with a 
reflection filter to control the transmission at the long 
wavelength end of the bandpass constitutes the basis 
of our design of narrowband filters. The two filters 
that we designed and evaluated to demonstrate the 
approach have bandwidths smaller than 5 nm and 
peak transmittances higher than 25% for collimated 
incident light. The average transmittance in the re- 
gion of longer wavelengths is <3% while the transmit- 
tance in the region of shorter wavelength is <0.1%. 

The values of the extinction coefficients of BaF2 and 
LaF 1ms are much smaller at longer wavelengths 

than 135 nm. ” hich makes the design of a narrow- 

band high reflectur much easier. Other film materials 
such as AloOu and HfO-i or suitable bulk materials 
(substrates) such as BaF?, fused silica, and AUO3 can 
be j d to reject different shorter wavelength ranges. 
Wht.i such materials are combined with reflection 
multilayers, narrowband filters with bandwidths 
smaller than 5 nm and overall transmittances higher 
than 25% can be made for the whole VUV region from 
120 to 230 nm. 
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Fig. 20. Combined filter centered ... 141 nm. 


The discrepancies between theory and experiment 
may be partially explained by the random errors which 
occur in the film thickness monitoring during deposi- 
tion. There are certainly some other effects to consid- 
er such as scattering, film thickness nonuniformities, 
and possible film inhomogeneities. Since the discrep- 
ancy between theory and experiment varied for differ- 
ent designs, it seems that these effects do not affect 
each of them equally. 

The first author would like to thank the Physics 
Department of the University of Alabama in Hunts- 
ville (UAH) for the graduate student assistantship 
received during this study, and the ZRAK Optics Co. in 
Sarajevo, Yugoslavia, for additional financial support. 

This work was supported bv NASA contracts NAS8- 
>7586 and NAS8-37576 and NASA grants NAG8-639 
and NAG8-086. 


References 

L. D. F. Heath and P. A. Sacher, ‘Effects of a Simulated High- 
Energy Space Environment on the Ultraviolet Transmittance of 
Optical Materials Between 1050 A and 3000 A.” Appi. Opt. 5, 937- 
943(1966). 

2. M. Zukic, D. G. Torr, J. F. Spann, and M. R. Torr, ’‘Vacuum 

Ultraviolet Thin Films. 1: Optical Constants of BaF 2 , CaF 2l 

LaFi, MgF 2 , AlsOs, Hf0 2 , and Si0 2 Thin Films,” Appi. Opt. 29, 
0000-0000 ( 1990). same issue. 

3. M. Born and E. Wolf, Principles of Optics t Pergamon, New York, 
1980), p. 60. 

4. Ref. 3, p. 616. 

5. Z. Knittl, Optics of Thin Films (Wiley, New York, 1976), Chaps. 5 
and 6. 

6. H. A. Macleod, Thin-Film Optical Filters (Adam Hilger, Bristol, 
1969), Chaps. 2 and 7. 

7. Ref. 3, p. 347. 

8. Ref. 6, p. 98. 


PASE ls> 

OF POOR QUALITY 


4302 


APPLIED OPTICS / Vol. 29, No. 28 / 1 October 1990 



APPENDIX D 



JOURNAL OF GEOPHYSICAL RESEARCH. VOL. 95. NO AT. PAGES 10.337-10.344. JULY L 1 9% 


Auroral Modeling of the 3371 A Emission Rate: 
Dependence on Characteristic Electron Energy 

P. G. Richards 

Computer Science Department and Center for Space Plasma and Aeronomtc Research 
University of Alabama in Huntsville 


D. G. Torr 

Physics Department and Center for Space Plasma and Aeronomtc Research 
University of Alabama tn Huntsville 


We have developed nn efficient two-strewn auroral electron model to study the deposition of 
auroral energy and the dependence of auroral emission rates on characteristic energy. This model 
incorporates the concept of average energy loss to reduce the computation time. Our simple two- 
stream model produces integrated emission rates that are in excellent agreement with the much 
more complex muitistream model of Strickland et ol. (1983) but disagrees with a recent study 
by Rees and Lmnmerzheim (1989) that indicates that the N 2 second positive emission rate is a 
strongly decreasing function of the characteristic energy of the precipitating flux. Our calculations 
reveal that a 10 keV electron will undergo approximately 160 ionizing collisions with an average 
energy loss per collision of 62 eV before thermalizing. The secondary electrons are created with 
an average energy of 42 eV When all processes including the backscattered escape fluxes are 
taken into account, the average energy loss per electron-ion pair is 35 eV in good agreement with 
laboratory results. 


1. Introduction 

There is currently renewed interest in the use of au- 
roral optical emission rates to deduce the characteris- 
tics of the precipitating particle fluxes, and ultimately, 
the global auroral energy input to the Earth’s upper 
atmosphere. Images from t lie Dynamics Explorer satel- 
lite have been used by Rees et al. [1988) to calculate 
the energetic electron flux and its characteristic energy. 
Imaging instruments planned for the ISTP mission will 
monitor key UV emissions on a global scale for the ex- 
press purpose of determining the global energy input. 

Early work in determining auroral particle character- 
istics from emissions concentrated on the use of the ra- 
tios of atomic oxygen emission rates (6300 A, 5577 A) to 
molecular nitrogen ion emission rates (3914 A, 4278 A) 
to deduce the incident auroral spectrum [Rees and 
Luckey, 1974; Vallance Jones, 1975; Shepherd et al., 
1980; Strickland et al., 1983). The higher energy auroral 
electrons penetrate deeper into t lie thermosphere where 
the relative proportion of atomic oxygen is smaller. 
Thus the ratio of atomic to molecular emission rates de- 
creases with increasing electron energy. Unfortunately, 
chemical processes play an important role in the atomic 
oxygen emissions and it is difficult to separate the ef- 
fects caused by the characteristics of the auroral energy 
flux from the effects caused by changes in the atmo- 
spheric composition. Therefore, it would be useful to 
find an emission rate ratio that is sensitive to the au- 
roral characteristics but which is not complicated by 
chemical factors. 
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Recently, Rees and Lummerzheim [1989] suggested 
that the ratio of the second positive to first negative 
emission rates could be used to determine the charac- 
teristic energy of the auroral electron flux. Using an 
auroral electron model developed by Lummerzheim et 
al. [1989], Rees and Lummerzheim [1989] found that 
the N 2 second positive (3371 A ) emission rate de- 
creases substantially with increasing characteristic en- 
ergy of the auroral electrons while the Nt emission 
rates are almost constant. This ratio would be an at- 
tractive alternative to those used previously because it 
would be independent of atmospheric composition and 
both emissions are prompt, thus eliminating chemical 
effects. Unfortunately, the calculations of Rees and 
Lummerzheim [1989] are in conflict with the earlier cal- 
culations by Darnell and Strickland [1986] who found 
that the 3371 A emission rate was nearly independent 
of the characteristic energy. 

The experimental evidence also seems to be in con- 
flict. Rees and Lummerzheim [1989] present data from 
high flying aircraft that support their theoretical calcu- 
lations. On the other hand, Solomon [1989] presented 
data from the visible airglow instrument on the At- 
mosphere Explorer C satellite showing that the ratio 
of the N 2 3371 A to N* 4278 A emission rates has 
only a small dependence on the characteristic energy, 
which can be accounted for by contamination of the 
3371 A second positive emission by the Vegard-Kaplan 
(0-9) band. The VAE data support the earlier calcula- 
tions of Damell and Strickland [1986] and Strickland et 
al. [1983]. Solomon was able to reproduce the observed 
ratios using his own two-stream auroral electron depo- 
sition code. We note that the experimental data pre- 
sented by Solomon [1989] for the ratio of N 2 3371 A to 
NT 4278 A is in excellent agreement with the ratio of 
N 2 3371 A to NT 3914 A that was measured on a 1974 
rocket flight by Sharp et al. [1979]. 
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The source of the discrepancy between the models is For each of the 3 main neutral atmospheric compo- 
difficuit to understand. Both Rees and Lummerzheim nents (O, O 2 , N 2 ) the total cross sections are made up 
[1989] and Strickland et ai [1983] use relatively com- of partial cross sections from the numerous electronic 
plete multi-stream models. The model used by Solomon states of each species each corresponding to a different 
[1989] is a simpler two-stream model and his calcula- energy state of the atom or molecule. In addition, each 
tio; s are in agreement with the calculations of Strick- electronic state of a molecule can be created in any one 
land et al. [1983] and Darnell id St: ckland [1986] of a large number of vibrational energy states. A com- 
wliich show little variation in the 3371 A to plete evaluation of the electron flux would require the 
3914 A emission rate ratio as a function of characteristic separate accounting of all these energy losses which we 
energy. In this paper we present calculations of a num- have replaced with a single, energy dependent average 
ber of important emission rates obtained from a two- energy loss for the excitation and ionization of each of 
stream auroral electron model that we have developed, the three main thermospheric species O, 0 2l and N 3 . 
The calculated emission rates for the 3371 A emission The treatment of excitation processes is relatively 
rate and a number of other emission rates are in accord straightforward; the average energy loss is specified and 


with the earlier calculations of Strickland et ai [1983] 
and Darnell and Strickland [1986] for the dependence 
on characteristic energy of the precipitating flux. How- 
ever, mainly due to the use of a revised cross section, 
our atomic oxygen 1356 A emission rate is a factor 
of 2.5 lower than that calculated by Strickland et ai 
[19831. 

2. Model 

2.1. General Principles 

The model that we have developed is based on the 
two-stream photoelectron flux model of Nagy and Banks 
[1970] that was subsequently extended to 500 eV and 
combined with a continuous energy loss model to calcu- 
late auroral electron fluxes by Banks et ai [1974]. Al- 
though both our model and the model of Solomon et ai 
[1988] have origins in the Nagy and Banks [1970] two- 
stream model, they have evolved substantially along en- 
tirely separate paths. Our model owes much to our ear- 
lier work with the ionospheric photoelectron flux 
[Richards and Torr , 1984; 1985a]. 

By incorporating a variable energy grid developed by 
Swartz [1985] and a variable altitude grid we have been 
able to extend the two-stream model up to energies 
greater than 20 keV. This lias eliminated problems en- 
countered by Banks et ai [1974] in matching the con- 
tinuous slowing down approach that they used above 
500 eV with the two-stream approach they used below 
500 eV. The continuously variable altitude grid allows 
altitude steps of less than a kilometer below 90 km up 
to 50 km at 500 km with a manageable number of grid 
points. 

Further economy in computer time is achieved by in- 
troducing the concept of an average energy loss for an 
excitation or an ionization event. This concept allows 
the use of only the total excitation and ionization cross 
sections instead of treating each partial excitation or 
ionization process separately. That is, the excitation 
(or ionization) is treated as arising from a single aver- 
age state for each thermospheric species. The average 
energy loss depends on the species and also on the en- 
ergy of the primary electron. The calculation of the 
emission rates then becomes a two stage process with 
just the total cross sections being needed to calculate 
the electron flux as a function of energy and altitude 
in the first stage. In the second stage, the partial cross 
sections are folded with the electron fluxes to produce 
the excitation rates. 


the electrons deposited in the correct lower energy bin. 
However, ionizing collisions are more complex because 
of the production of secondary electrons which may be 
produced with energies ranging from 0 up to E p — I{ 
where E v is the energy of the primary electron and U is 
the ionization potential of the state i being produced. 
We follow the approach of Banks et ai [1974] and des- 
ignate the higher energy electron as the degraded pri- 
mary and the lower energy electron as the secondary. 
This means that the maximum secondary energy is then 
(E v — /, )/2 as is the minimum energy of the degraded 
primary. We treat ionization as arising from a single 
state with an average ionization potential I and we use 
the measured secondary electron distributions of Opal 
et ai [1971] to determine the average energy E, of the 
secondary electrons produced by an electron of energy 
E v (we note that the measured secondary electron dis- 
tributions are in fact a sum of the contributions from 
all the ionization states). The average energy of the 
degraded primaries is then E v — I — E 9 , A separate 
ionization potential could be used for each species but 
a further improvement in computational efficiency can 
be made by observing that the ionization potentials and 
secondary electron distributions are similar enough that 
a single ionization potential and secondary electron dis- 
tribution for all 3 major species will suffice. Since N 2 is 
the most important species in aurora, we adopt the N 2 
ionization characteristics. At each electron energy, the 
total number of secondary electrons from O, 0 2l and N 3 
is calculated and then they are distributed according to 
the measured secondary electron distributions of Opal 
et ai [1971]. The model has been found to conserve 
energy to better than 5% for characteristic energies in 
the range .1 to 20 keV. 

2.2. Average Ionization Potentials 

The calculation of the average ionization potential 
follows from the knowledge of the ionization potentials 
and the partial cross sections for the various ionization 
states of each of the thermospheric species. For exam- 
ple, N 2 ionization results in the formation of the X, A, 
and B states of Nt and also N + which arises from sev- 
eral higher lying states that dissociate. The ionization 
potentials for the X, A, and B states are 15.6, 16.8, and 
18.8 eV respectively while the bulk of the N + probably 
arises from a state with a threshold near 37 eV [ Erd - 
man and Zipf \ 1986]. The actual energy lost by the 
primary electrons may be greater than threshold due 
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to the ions being vibrationally excited and because, in 
the case of dissociative ionization, the atomic fragments 
are observed to carry substantia] kinetic energy [Lil- 
iana and Stoc dale, 1975]. To complete the calculation 
of the average ionization energy for No we also need 
the ratios of the partk*. cross sections. We examine 
the high energy case first where the proportions of the 
various partial cross sections are practically constant 
and the calculation is straightforward. Above approxi- 
mately 100 e V\ dissociation accounts for — 20% [Rapp et 
a/., 1965] and the B state —10% [5or$t and Zipf } 1970] 
of the total cross section. According to Cartwright et 
a l. [1975], the X and A state cross sections are ap- 
proximately equal (that is, —35%). Using these per- 
centages for the partial cross sections and ionization 
potentials given above, the average ionization potential 
is —20 eV for high energy electron impact ionization of 
N 2 . We now examine the low energy average energy loss 
per collision. The average ionization potential begins 
at 15.6 eV at threshold when the X state is produced 
but increases to —16 eV above 17 eV when the A stat’e 
threshold is reached. There is only a marginal increase 
when the B state threshold is reached at 18.8 eV because 
it accounts for less than 10% of the total cross section 
at these energies and the energy loss of 18.8 eV is only 
marginally larger than 16.8 eV. After about 30 eV dis- 
sociative ionization becomes important. The ionization 
potential then rises steadily to approximately 19 eV at 
50 eV before leveling off toward 20 eV at high energies. 

A similar calculation yields an average ionization po- 
tential of approximately 18 eV for both O and Oo, at 
high energies. Because is the dominant constituent 
and all three ionization potentials are similar, we use 
he average N 2 ionization potential for all three species. 
This simplifies the calculation and reduces the compu- 
tation time without introducing significant errors. The 
adopted ionization and excitation energy losses as a 
function of primary electron energy are shown in Fig- 
ure 1. This figure also shows the average energy of the 
secondary electrons as a function of the primary en- 
ergy. Below 25 eV the ionization potential is assumed 
to be 16 eV while above 25 eV it is represented by 
10 (1 + (1 - 15 /E) l/2 ) eV. 


2.3. Average Secondary Electron Energy 

To obtain the secondary electron distribution we 
adopt the Banks et al. [1974] parameterization of the 
Opal et al. [1971] measured distributions. The proba- 
bility of a secondary electron of energy E. is given by 


P,(E.) 


A 

1 + (£./£) 2 


( 1 ) 


where A = [(E arctan( E m /E) -1 ] is a normalization fac- 
tor that ensures a total probability of unity when in- 
tegrated over all secondary energies from 0 to E in — 
( E v - J)/2 and E=14 is an empirical normalization fac- 
tor. The product of E , and Equation 1 integrated over 
energy yields the average secondary electron energy for 
a primary energy E v as 


( 2 ) 



Fig. 1. Average energy losses per collision nvernge secondary 
electron energies (E a ») as a function of primary electron energy. 
The average ionization potential is labelled I and the total en- 
ergy loss per ionization is labelled I+E«» . The average excitation 
potentials are indicated by an asterisk and the 0 3 excitation po- 
tential is set equal to that of N 2 . Note that below 5 eV, the N 3 
excitation potential is set at 1 eV . 


The dependence of the average secondary energy on 
the primary energy is shown in Figure 1 along with 
the average total energy loss of the primary (/ + E tlv ). 
The average energy of the secondaries increases steadily 
from 0 near threshold to 52 eV for 10 keV primary elec- 
trons. The total energy loss per ionizing collision for 
10 keV electrons is 72 eV. when the 20 eV ionization 
potential is included. The average energy of the secon- 
daries is approximately equal to the ionization energy 
for 200 eV primary electrons. 

A quantity of interest, in relation to energy degra- 
dation of high energy electrons, is the average energy 
required to produce each electron-ion pair. The aver- 
age energy required to produce eicli electron-ion pair 
is a quantity that is independent of electron energy 
and is also remarkably independent of the species be- 
ing ionized. The experimental value for the energy lost 
per electron-ion pair for high energy electron is 35 eV 
[Valentine and Curran , 1958]. This energy loss per 
electron-ion pair was used in early auroral electron de- 
position codes [Rees, 1963: Rees et al. 1969; Rees and 
Jones , 1973]. 

An approximate value for the energy loss per electron- 
ion pair can be deduced by using the average energy 
losses depicted in Figure 1, assuming that for electron 
energies above approximately 100 eV the energy lost 
to excitation collisions is small and can be neglected. 
With this assumption, it requires about 160 ionizing 
collisions to thermalize a single 10 keV electron. Thus, 
on the average, 62 eV is lost in the creation of each of 
the 160 electron-ion pairs. Since the ionization energy 
is 20 eV for electron energies above 100 eV, the average 
energy of the secondaries is —42 eV. This means that 


E„ = 0.5 A E- In (l + ( E,/E ) 2 ) 
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the 160 secondary electrons are sufficiently energetic to 
create one more electron-ion pair each. If all secon- 
daries created an additional electron-ion pair, the total 
number of pairs would be 320 and the average energy 
per electron-ion pair would decrease from 62 to 31 eV. 
In reality this does not happen because, below 100 eV, 
excitation processes begin to compete effectively with 
the ionization processes for the available electron en- 
ergy and the number of additional electron-ion pairs 
produced by the secondary electrons would be less than 
160. In fact, for a 42 eV electron the total ionization 
and excitation cross sections are approximately equal 
and only half the secondaries could be expected to pro- 
duce an additional electron-ion pair. This agrees with 
our previous calculations that show that a third rather 
than a half, of the total ionization is created by de- 
graded primaries and secondaries with energies below 
100 eV [ Richards and Torr , 1985b]. Thus the original 
10 keV electron would ultimately produce about 240 
electron ion pairs and yield an average energy loss per 
electron-ion pair of 42 eV. This is only an estimate of 
the energy lost per electron-ion pair and a more detailed 
calculation including transport is required to determine 
the actual value. It was pointed out by Banks et ai 
[1971] that escaping backscattered electrons will be lost 
to the system and act to increase the energy loss per 
electron-ion pair. 

We have summed the total ion production rate in our 
full auroral calculation, and we obtain an average en- 
ergy loss per electron-ion pair of 35 eV, which is smaller 
than our estimate but in agreement with the laboratory 
measured value. A slightly higher value of 37 eV was 
obtained by Fox and Victor [1988] using their discrete 
local energy loss method. The reason that the energy 
lost per electron-ion pair is not a strong function of 
electron energy has to do with the relationship between 
the average secondary electron energy and the primary 
electron energy. Electrons with higher initial energies 
suffer a greater energy loss per collision as they degrade 
but they produce higher energy secondaries which are 
more likely to generate secondary ionization. For ex- 
ample, a 1 keV electron will undergo only 22 ionizing 
collisions with an average energy loss of 45 eV before 
it thermalizes. Thus, the secondary electrons have an 
average energy of only 25 eV compared to the 42 eV for 
the 10 keV electrons and are much less likely to produce 
additional ions. 

2.4. Cross Sections 

In a number of previous studies of the ionospheric 
photoelectron flux we have chosen measured total exci- 
tation cross sections that produce good agreement be- 
tween theory and the photoelectron spectrometer mea- 
surements from the AE-E satellite [Lee et a/., 1980]. At 
low altitudes where N 2 is the dominant species, the to- 
tal cross section obtained from electron mobility studies 
by Pitchford and Phelps [1982] was found to be com- 
patible with the PES measurements. At high altitudes 
where atomic oxygen is the dominant species, the emis- 
sion cross sections measured by Zipf and co-workers pro- 
duced good agreement between theory and PES mea- 
surements. These total cross sections for energies be- 


low 100 eV have been published by Rickards and Torr 
[1988]. Basically the same cross sections have been used 
in this study but they have been reparameterized to ex- 
tend them to higher energies. 

Above 100 eV, the excitation cross sections decay 
rapidly with increasing energy and are much less im- 
portant than the ionization cross sections both because 
they are smaller and because the energy loss per colli- 
sion is smaller. The total ionization cross sections are 
better established than the total excitation cross sec- 
tions although there are some differences [Kieffer and 
Dunn , 1966]. We have adopted the N 2 and 0 2 total 
ionization cross sections of Rapp and Eng lander- Golden 
[1965] which have also been used by most other mod- 
elers. The total ionization cross section for O is from 
Brook et a/. [1978]. The elastic cross sections are very 
important because of their role in inhibiting transport. 
We have used the elastic cross sections of Solomon et 
al. [1988] and also their elastic backscatter coefficients. 
The total cross sections used in our auroral model are 
shown in Figure 2. Our N 2 total excitation cross section 
is comparable to that of Solomon et al. [1988] below 
25 eV but is smaller at higher energies. The differences 
at high energies have little effect on the calculated fluxes 
because the excitation cross section is smaller than the 
ionization cross section. However, differences in cross 
sections below 30 eV produce comparable differences in 
fluxes. The N 2 total excitation cross section of Strick- 
land et al. [1983] is almost a factor of two larger than 
ours at all energies and their fluxes would be a fac- 
tor of two lower below 30 eV, at least below 200 km 
where N 2 is the dominant species. The cross sections 
of Solomon et al. [1988] and Strickland et al. [1983] 
were obtained by summing the partial cross sections 
and there is the possibility of double counting some 
cross sections; for example, those that lead to dissocia- 
tion. Moreover, Strickland et al. [1983] included large 



Fig. 2. Total elastic, excitation, and ionisation cross sections em- 
ployed in the model. The ionization cross sections are indicated 
by a plus while the excitation cross sections are indicated by an 
asterisk. 
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Rydberg cross sections from Green and Stolarski [1972] 
which have been revised sharply downward by Porter 
et ai [1976]. There is now good agreement between 
the total cross section of Pitchford and Phelps [1982] 
and the sum of the partial excitation cross sections of 
Cartwright et al. [I977a,b] as revised by Trajmar et ai 
[1983], below 20 eV. Above 20 eV, the ionization cross 
section becomes an increasingly important component 
of the total inelastic cross section and it is not easy to 
compare the two cross sections. 

The total atomic oxygen excitation cross section em- 
ployed by Solomon et al. [1988] is a factor of 2 larger 
than ours above 15 eV and will produce a similar dif- 
ference in flux above 250 km where O is the dominant 
species but the atomic oxygen cross section has little 
effect on the integrated emission rates. Our atomic oxy- 
gen excitation cross section was obtained by summing 
the measured emission cross sections for 1304, 1356, 
and 1027 A [ Zipf and Erdman , 1985], the 989 cross 
section from Gulcicek and Doering [1988], and the the- 
oretical 1 D and L S from Henry et ai [1969]. Implicit in 
this procedure is the assumption that the higher lying 
triplet and quintet states are included in the 1304 and 
1356 emission cross sections via cascade. We have left 
out some theoretical Rydberg cross sections proposed 
by Jackman et al. [1977] and some minor states that 
radiate directly to the ground state but for which there 
is no experimental data. Thus, our cross section must 
be regarded as a lower limit. 

Figure 3 shows the excitation cross sections for the 
second positive (C 3 7 t„) and Lyman- Birge-Hopfleld 
(a 1 7r f/ ) systems of N 3 . Also shown is the cross section 
used for calculation of the 0( 5 S) 1356 A emission rate. 
We obtained t his cross section by reducing the measured 
cross section of Stone and Zipf [ 1974] by the factor 3.1 
which is the same factor that the 1304 A emission cross 
section of Stone and Zipf{ 1974] was reduced by Zipf and 
Erdman [1985]. The (C 3 tt (< ) cross section was obtained 



Fig. 3. Cross sections for the three excited states giving rise to 
the emissions studied in this paper. 


by multiplying the 3371 A cross section of Imami and 
Borst [1974] by 4 and the (a l x,j) cross section is from 
Ajello and Shemansky [1985]. 

We have examined the sensitivity of the emission rate 
ratios to cross sections and this will be discussed later. 
In ail these calculations we have used the 1 erg cm'V 1 
Gaussian incident flux distribution, and the neutral at- 
mosphere employed by Strickland et al. [1983]. 

3. Results 

3.1. Comparison With Previous Work 

We have calculated the N 3 3371 A , 3914 A , N 2 

4278 A , O 1356 A , and several N? LBH band emission 
rates as a function of energy and these are shown in 
Figure 4. This figure shows that both the 3371 A and 
3914 A emission rates are independent of the char- 
acteristic energy of the precipitating flux for energies 
above 2 keV in agreement with the results of Strickland 
et ai [1983] and Damell and Strickland [1986]. Not 
only is the shape in good agreement but, except for the 
1356 A emission rates, the magnitudes are also in good 
agreement. Although the shape of the 1356 A curve is 
in good agreement with that of Strickland et al. [1983], 
the magnitude is a factor of 2.5 lower owing to the use 
of the revised cross section of Zipf and Erdman [1985]. 



ENERGY (keV) 


Fig. 4. Calculated emission rotes ns n function of the chnracteris- 
tic energy for n Gaussian energy distribution with o total incident 
energy flux of 1 ergem^s -1 . When differences in cross sections 
are tnken into account, there is excellent agreement with the cal- 
culations in Figure 8 of Strickland et al. (1983) and Figure 11 of 
Damell and Strickland (1986). 
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We have included 0 2 Schuinan-Runge absorption which 
affects both the 1356 A and the 1200-1600 A LBH 
bauds when the characteristic energies are high and the 
electrons penetrate to lower altitudes [Strickland et al. y 
1983]. The ratios of the LBH bands are taken from 
Ajello and Shemansky [1985]. 

Our calculated second positive to first negative emis- 
sion rate ratios are within 20% of the measured val- 
ues for the Nt emissions. At 10 keV, our 3371 A to 
3914 A ratio is 0.3 compared to 0.25 from Sharp et al. 
[1979] and our 3371 A to 4278 A ratio is 0.98 compared 
to 0.8 from Solomon [1989]. Solomon obtained better 
agreement between his model ratios and the measured 
ratios but there is sufficient uncertainty in the input 
parameters to account for the differences. 

3.2. Sensitivity of Ratios 



We have performed some parameter studies to char- 
acterize the sensitivity of the ratio of the second positive 
to first negative integrated emission rates to possible 
errors in the model inputs. Obviously, a reduction of 
20% in the 3371 A emission rate cross section would 
bring the calculated and measured values into excellent 
agreement but a 30% increase in the Ns total excitation 
cross section has a similar effect by decreasing the low 
energy electron flux which is responsible for most of the 
3371 A emission. Likewise, a 30% increase in the N 2 
total ionization cross section above 100 eV reduces the 
ratio from 0.3 to 0.25 by increasing the 3914 A pro- 
duct 1 n rate. The integrated ratio is not sensitive to 
changes of up to a factor of two in most other param- 
eters including: the atomic and molecular oxygen in- 
elastic cross sections, the O, 0 2 , and N 2 elastic cross 
sections and backscatter coefficients, and the relative 
concentrations of the species. We estimate a possible 
error of 10% in our computed average excitation and 
ionization potentials but t his has negligible effect on 
our computed ratios. 


ENERGY (ktV) 

Fig. 5. The initial partitioning of the incident 1 erg cm -3 s" 1 
energy flux between ionisation, excitation, thermal electron heat- 
ing, and backscatter as n function of characteristic energy. The 
largest proportion of the energy (—35%) goes inr.iaUy into the 
ionization potential of the while f — 20% ) goes into Nj excita- 
tion. Only ( — 16%) is backscattered out of the thermosphere. O 
is an important absorber of energy at the lowest energies while 
Oj becomes increasingly important as the characteristic energy 
ncreases. 


3.4. Electron Flux Spectra 

Downward moving fluxes at 120, 174, 223, and 326 
km are shown in Figure 6 for a 5 keV incident Gaussian 
flux with an energy flux of 1 erg cm“ 2 s“ 1 . The incident 
flux can be seen centered at 5 keV in Figure 6. At the 
two highest altitudes, there is very little degradation 
of this initial flux but the degradation is noticeable at 
174 km and pronounced at 120 km. Because there is so 


3.3. Energy Budget 

T!:*» incident electron energy flux is initially parti- 
tioned into a large number of excitation and ionization 
processes before it finally emerges as heat for the *!ier- 
mosphere or is radiated into space. Figure 5 shows the 
gross energy partition amongst ionizations, excitations, 
thermal iectron heating and backscattered escape flux 
as a function of characteristic energy. N 2 ions capture 
the greatest share of the available energy (35%). Exci- 
tation of N 2 is next with (20%) while only 15% is lost 
through lie escaping backscattered flux. This escape 
flux is cli smaller than the 45% obtained by Banks 

et al. , ~4], possibly as a result of the use of different 

cross sec: ions and backscatter coefficients. Below 1 keV, 
* uization and excitation of atomic oxygen absorption 

e important energy sinks for the electron energy, but 
they become small for high energy incident fluxes. Ab- 
sorption in molecular oxygen shows the opposite trend, 
becoming more important with the deeper penetration 
of the higher energy fluxes. Thermal electrons capture 
a greater proportion of the available energy, the lower 
the characteristic en rgy. 



Fig. 6. Downward moving electron flux spectra at several alti- 
tudes for a 5 keV Gaussian incident flux. The incident energy 
flux is 1 erg cm" 2 s -1 
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little interaction with the thermosphere for high energy 
electrons at the high altitudes, there are few degraded 
primaries to fill in the region between 300 eV and 5 keV. 
However, at the lowest altitudes this intermediate en- 
ergy range is filled in. Comparison of the downward 
fluxes in Figure 6 with the upward fluxes in Figure 7 
reveals that, below 225 km, where transport is inhib- 
ited, the electron flux is isotropic for energies less than 
300 eV. At 326 km, the upward (escape) flux is a factor 
of 2 larger than the downward flux at low energies and 
orders of magnitude larger at intermediate energies. 

4. Conclusions 

We have developed an efficient two-stream auroral 
electron model that incorporates the concept of aver- 
age energy loss. This model produces integrated emis- 
sion rates that are in excellent agreement with the more 
sophisticated multi-stream model of Strickland et ai 
[1983] but is in disagreement with the model of Rees and 
Lummerzheim [1989] with regards to the energy depen- 
dence of the N 2 3371 A second positive emission rate. 
Our calculations give a value of 35 eV for the average 
energy lost per electron-ion pair produced independent 
of primary electron energy and we have explained this 
behavior in terms of the variation in the energy of the 
secondary electrons. We find that more than 30% of 
the initial energy flux is stored initially as ionization 
energy of Nt while about 20% goes into excited states 
of N 2 while only 15% is backscattered out of the ther- 
mosphere. All other processes are minor except at low 
incident energies where 20% of the energy is stored in 
atomic oxygen ions. 
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Fig. 7. Upward moving electron flux spectra for the 5 keV Gaus- 
sian incident flux shown in Figure 6. 
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Images of the entire auroral oval at carefully selected wavelengths contain information on the global energy 
influx due to energetic panicles and some information on the characteristic energy ot the precipitating panicles. 
In this paper we investigate the sensitivity ot selected auroral emissions to changes in the neutral atmosphere. In 
particular, we examine the behavior of 01 1356 A and two Lyman Birge Hoptield (LBH) bands and their ratios 
to each other with changing atmospheric composition The two LBH bands are selected so that one lies in the 
region ot strong O? absorption 1 1464 Ai and one lies at a wavelength where O absorption is effectively negli- 
gible ( 1838 A) We find that for anticipated average uncertainties in the neutral atmosphere t factor ot 2 at aurorai 
altitudes), the resultant chance in the modeled intensities is comparable to or less than the uncertainty in the 
neutral atmosphere. The smallest variations, tor example, arc tor I 1 83 8 (approximately iO to 20#) while the 
largest variation is seen in the Ol 1356 A emission which is linear with [O] to within 20# We have also 
investigated the dependence ot these intensities, and their ratios, to much larger changes in the composition tie.. 
[0|/|N^]) such as mieht be encountered in large magnetic storms, or over seasonal or solar cycle extremes. We 
find that the variation in the 1 1356/1 1838 ratio over the equivalent ot a solar cycle is less than 50#. The 
summer-to-winter changes are approximately a lactor ot 2. The 1 1 356/1 1 838 ratio is a very sensitive indicator ol 
the characteristic energv. showing a change ot 13 over the energy range 200 eV to 10 keY The corresponding 
change in ihe LBH long-to-short wavelength ratio is much less (about a lactor ot 3). However, the latter is 
insensitive to changes in the neutral atmosphere t < 20# changes in LBH emission ratio tor large changes in NO. 
The three emissions therefore potentially provide a most valuable diagnostic ot particle characteristic energy and 
energy flux. 


I Introduction 

While in situ observations of energetic particles provide 
accurate information on the panicle characteristics at the point of 
measurement, imaging from space of the enure auroral oval holds 
the potential for providing details on total auroral energy influx, 
estimates of the characteristic energy of the auroral particles, and 
the capability to map and relate the footprint of this derived 
information hack along the magnetic field lines to various regions 
of the magnetosphere Auroral imaging in the vacuum ultraviolet 
permits observations of the regions of interest under both day and 
night conditions. Work by Rees and Luckey 1 19741 on the ratios 
of visible emissions. UV emission intensity calculations by 
Strickland et al . [!983|. and analysis of UV auroral spectra by 
Isiumoti) et ai. |1988| all indicate the potential value of using 
ratios of emission intensities to study auroral processes. A major 
focus of work in this area at the present time is to establish the 
quantitative footing on which such determinations can be placed. 

With the exception of HI Lya. the Ol muluplets at 1304 A and 
1356 A and the N : Lyman Birge Hopfield (LBH) bands arc the 
most prominent vacuum ultraviolet auroral emissions. The Ol 
1304 A emission has a high efficiency for multiple scattering. As 
a result, it has limited use for actual auroral imaging, although it 
does have potential value as an indicator of the O concentration. 
While the 1356 A emission does undergo multiple scattering, the 
efficiency is relatively small [Strickland and Anderson. 19831 and 
we ignore multiple scattering for 1 1356 for this study. Similar 
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considerations allow us to also ignore multiple scattering for the 
Na LBH emissions that are also considered in this study. The Ol 
1356 A emission is absorbed increasingly by O;* with decreasing 
altitude. Thus its intensity varies strongly (inversely) with 
increasing depth of penetration of the incident aurorai electrons 
and hcncc with increasing energy. The NU LBH transitions are 
electric dipole forbidden and the only prominent excitation mech- 
anism ts electron impact. The LBH emission may therefore serve 
as a direct measure of the total energy tlux of charged particles 
into the atmosphere. The longer wavelength LBH bands, which 
lie outside the region of substantial O? absorption, are useful 
indicators of the total energy influx, while the long-to-short 
wavelength LBH intensity ratio provides information on the 
and thus also some information on energy. These are the 
emissions (Ol 1356. long and short wavelength LBH) on which 
wc shall concentrate in this study. 

The purpose of this paper is to examine the sensitivity of these 
emissions to both likely uncertainties and anticipated changes in 
the neutral atmosphere. This is just one step tn the process of mak- 
ing quantitative interpretations of auroral images, but an 
important one. We will consider other aspects (energy spectral 
characteristics and wavelength spectral extraction) elsewhere. In 
this paper w'e conduct a series of sensitivity studies using an 
auroral emission code that has been developed by our group 
[Richards and Torr, 1990], The results are discussed below. 

2. Description of Auroral Code 

The behavior of auroral Ol 1356 and N : LBH emissions has 
been studied with the use of an aurorai computer model. The 
model is a two-stream auroral electron energy loss code that deter- 
mines the energy degradation of the primary spectrum as a func- 
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non ot energy and altitude and determines the production rates of 
prominent auroral emissions. A more complete description ot the 
program, the selected cross sections, and comparison with other 
auroral models is given by Richar d and Torr \ 1990). The 
N : ( a 1 tt , ) cross section, as well as tl * mhos ot the individual LBH 
bands are taken from A/ello and Sm manskx 1 1 9S5 ] . Attenuation 
due to molecular oxygen absorption is explicitly computed with 
0 2 absorption cross sections taken from Ogawa and Ogawa 
1 1975] and Hudson 1 1971 ). The model currently assumes that 0 2 
absorption at wavelengths beyond 1750 A can be ignored. 

The model is optimized by incorporating variable energy bins 
|7Y;rr et al. 1974; Swartz , 1985] for the energy grid. To prevent 
numerical instabilities and violation of energy conservation, the 
altitude grid is variable to allow small grid steps (less than 1 kmi 
at low altitudes. As a result, energy is typically conserved 
to within 10%. The code utilizes either the MSIS-86 neutral 
atmosphere \Hedm, 1987] or a user-supplied atmosphere. Either 
monoenergetic fluxes or a specified energy spectrum may be 
used. The incident energy spectrum may be modeled as a 
Gaussian or Maxwellian distribution, after Strickland et al. 

1 1983| (hereafter SJW). or a user-supplied distribution may be 
used. All simulations reported below employed a 1 erg cm * s 1 
Gaussian incident flux distribution. The Gaussian scale parame- 
ter. labeled W in 57VV. has been set equai to 0.25 E char which 
yields a full width at half maximum of 0.5 [ 1 n 2 ] 1 ~ E vhar , where 
E^har is the characteristic energy. 

3. Sensitivity Studies 

The emission studies reported here involved modeling auroral 
emissions at local midnight at 60 degrees north latitude. Table I 
lists the MSI S parameters used in this study as well as detailing 
the range of solar activity investigated in the latter part of the 
study. 


TABLE I MS1S Model Parameters 


Solar Activiiv 


Minimum 

Moderate 

Maximum 

I'm - ci 11 Flux Index 


i ID 

:do 

\\erage F,,, < an Index 

75 

1 ID 

:oo 

\p Magnetic Activity Index 

4 

:o 

UK) 


Geographic lamude. 60 degrees; geographic longitude. 0 degrees, solar 
jpparent time. 0.0 hours; days 173. 356 


Three emission ratios were studied. The first ratio was Ol 1356/ 
LBHiony. where LBH (fmj , designates an N : LBH emission not 
strongly dominated by G : absorption. Specifically, the (2.8) band 
at 1838 A was chosen for this purpose. Second, the ratio OI 1356/ 
LBFf hon was also modeled to investigate the relative influence of 
absorption by molecular oxygen. Here. LBH shor1 is represented by 
the ( 1 , 1) band at 1464 A. By analogy with the previous definition 
LBFf hon is an LBH emission which is strongly absorbed by 0 : . 
The final ratio studied was LBH (on|r LBH, hnn . The volume emis- 
sion rates integrated over altitude give the surface brightness or 
column intensity of the emission which we shall designate 1 1356, 

I 1464. and I 1 838. These are the intensities that would he seen by 
a nadir viewing instrument from above the emission layer. 

3.1. Sensitivity to the Uncertainty of a Single Constituent 

The first question we chose to investigate was the dependence 
of the selected emissions and emission ratios on the uncertainties 
at any given time in our knowledge of the neutral atmosphere. We 
have assumed for this purpose that if we base our calculations on 
the MSIS-86 model atmosphere, the concentrations of O, 0 2 , or 
N 2 at auroral altitudes may on the average be uncertain by as much 


as a factor of 2. There will be occasions on which the uncertainty 
will exceed a factor of 2. but typically it will be less. This is 
similar to studies performed in SJW , but extends the investigation 
to study the dependence of auroral emissions to each of the major 
atmospheric constituents. In addition, in section 3.2 below, we 
further extend the study to include larger compositional variations 
due to seasonal and solar cyclic variations. 

The unperturbed, or reference, atmosphere is an MSIS model 
for hieh solar activity at summer solstice. Figure 1. Figure 2 



Number Oensity (cm'*) 

Fie. 1. MSIS-86 reference atmosphere used as the standard case in this 
study (Day = 173. F,<> 7 - -00. Ap - UK)) 



0 2 4 6 8 10 

Energy (kev) 

Fie. 2. Nadir viewing column brightnesses (I 1356. ( 1464. and I 1838) 
calculated using the reference atmosphere shown in Figure I The 
diamonds here and in the remainder of the figures show the selected 
Gaussian characteristic energies lor the incident electron energy 
distribution. 
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shows the computed i 1356. 1 1464, and 1 1 838 intensities tor this 
model atmosphere as a function of energy over the range 200 eV 
to 10 keV as determined by our auroral code. In what follows we 
shall compare these results with those obtained when each ot the 
atmospheric constituents (O. CT. NM is. in turn, individually 
multiplied by 2 at all altitudes, while the other two are held 
constant. 

Before examining the results ot these atmospheric changes, let 
us consider the possible impact of the changes. Doubling the con- 
centration of a constituent might at first be expected to double the 
effect ot that specie on the column brightness ot the monitored 
auroral emissions. In reality, however, there are a number or 
possible options in the interaction of the penetrating electrons and 
the atmospheric gases that render the situation more complex. An 
electron of a given initial energy will undergo a fixed number ot 
collisions in a particular gas before thermalizing. Thus in the very 
simple case of a single constituent atmosphere, changing the con- 
centration simply raises or lowers the altitude ot the peak energy 
loss (and peak emission). This example (single-constituent 
atmosphere) is representative ot those altitudes in which the con- 
centration of one atmospheric constituent dominates. From Figure 
1 it can be seen that (for the conditions chosen) O tends to 
dominate above 400 km and N? tends to dominate below 500 km. 
In a mixture of gases, increasing the concentration ot one specie 
relative to the others may also have the effect ot raising the 
penetration depth. However, the gases will compete tor incoming 
electrons in proportion to their mixing ratios and collision cross 
sections, and the ratio of the resuiting emissions changes 
accordingly. 

Figure 3 shows the volume emission rate profiles for 01 1356 A 
photons for the reference atmosphere case for selected energies. 
Only the very soft electrons (<200 eV) lose their energy in the 
altitude region where O is the major species. All other energies lie 
in altitude regions where the various gases can compete tor 
collisons with the precipitating electrons. Thus, lor example, 
increasing the concentration ot Ni will result in a decrease in the 
production ot OI I 356 A photons, because electrons that would 
have collided with O atoms now have an increased probability ot 
colliding with molecules. The results ot changing the O, CT. 
and N-* concentrations individually by a factor ot 2 are shown in 
Figures 4 and 5. 



Fig. 3 Volume emission rate profiles lor O! 1356 A tor ihc reference 
jimosphere shown in Figure l 


Figures 4<j anil 5a show the impact on the computed OI 1 356 A 
and LBH surface brightness or ihe doubling of |N,). The 1356 A 
intensity at 3 keV drops to b59c of the reference model, and then 
rises back to 86 c /r of the reference model by 10 keV. The I 1 356 is 
reduced at all modeled energies due to the fact that there are now 
more collisions with N-. with a corresponding decrease in the 
production ot O emissions. The decrease is not a full factor of 2 
because of the abundance of atomic oxygen at the higher altitudes. 
For the very lowest energies (highest altitudes), where |N,| is 
much less than |()|. the 1 1356 should tend to an intensity level 
unchanged from the standard case, as the emission is simply 
raised in altitude. The modeled emissions do not include initial 
energies less than 200 eV which would lose their energy above 
400 km. but the lowest energy emissions do show this trend. At 
the higher energies, the emission is produced primarily at alti- 
tudes where N, is the major constituent and increased N, does not 
result in a significant change in the competition between O and 
N,. 

For the LBH 1464 A emission, doubling the N, reduces the 
relative concentration of the dominant absorber. O,. Thus tor the 
higher energies which penetrate to greater depths, the emission 
from the increased N, overwhelms the 0 : absorption. Absorption 
by O, is not significant for the LBH 1 838 A emission: there is thus 
little dependence on the energy ot the incident electrons. 

Doublina the O, density (Figures 46 and 56) increases the 
absorption 'of I 1356 and I 1464 at the higher energies (lower 
altitudes) resulting in reduced column brightnesses. The LBH 
1 838 A emission is relatively unaffected by O, absorption and is 
influenced onlv by increased competition for collisions ot the 
energetic panicles with CL molecules. However, since IN,] 
remains the major specie relative to 1 0,1 , 1 1838 shows only small 
changes. 

The effect of doubling the ( O ] is show-n tor I 1356 in Figure 4c 
and for the LBH emissions in Figure 5c. For I 1356 the et feet is 
close to a factor of 2 increase m the emission at all energies, while 
for the LBH emissions there is almost no effect at all energies. For 
the very low energies (not modeled) where O is the major con- 
stituent. the effect of doubling the O is simply to raise the altitude 
of the 1356 A emission. The trend to an unchanged emission can 
be seen at the lower energies. For the LBH emissions, the etfect 
of doubling the O concentration is only seen at the very low 
energies (high altitudes) where the competition with N : is further 
increased. In the altitude regimes where N : is a larger component, 
the O does not play a significant role. 

From the results shown in Figures 4 and 5 it can be seen that 
uncertainties of a factor of 2 in any of the principal neutral atmos- 
pheric species translate into uncertainties ot less than 20% for I 
1838. The LBH 1464 A emission shows variations up to 70% for 
factor of 2 uncertainties in |N : ) and less than 40% variation due to 
other constituents. OI 1356 is weakly sensitive to changes in 0 2 
and N-». but varies almost in direct proportion to changes in O. 


3.2. Sensitivity to Larger Compositional Changes 

The neutral atmosphere exhibits relatively large compositional 
changes in the course of the seasonal, solar cyclic and magnetic 
storm variations. In this section we report the results of our 
assessment of the dependence on the computed emissions chosen 
for this study to changes of this magnitude. In order to simulate 
changes we have varied the input parameters ( F,„ 7 * A P- ^ of 
year) to the MSIS-86 model atmosphere, yielding neutral atmos- 
pheres at summer and winter solstice for conditions corresponding 
to low, moderate, and high solar activity (Table 1). The relative 
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Fie. 4. 01 1 356 dependence on | O) , [0>]. and | N : ] Each curve shows I 1 356 for the doubled constituent case divided by l 1356 
tram the standard (unchanged) case. 


compositional changes produced by these cases are illustrated in 
Figure 6. 

The variations in the emission ratio of 01 1356 A to LBH 1838 
A due to such composition changes are illustrated in Figure 7. 
This particular ratio is very sensitive to incident energy. For any 
given atmospheric conditions, the ratio vanes by a factor of 13 
over the energy range shown in the figure. What is interesting to 
note in Figure s that the variation due to the compositional 
changes produced by solar activity variations (low. moderate, and 
high F M ) 7 cm flux) are small (=£30%). while the vanations result- 
ing from compositional changes of the type produced by seasonal 
vanations are much larger (about a factor of 2). In Figure 8 we 
show the energy dependence for these various composition cases 
tor the individual 1 1 356 and 1 1838 intensities. The LBH intensity 


is relatively insensitive to the changing atmosphenc conditions, 
while the OI 1356 is found to be pnmanly responsible for the 
vanations shown in Figure 7. 

The reason for this can be seen in Figure 9 which shows the 
ratio of the individual concentration changes relative to the 
standard case. Figure 9a shows that the atomic oxygen concentra- 
tion (for altitudes below 300 km. which correspond to the initial 
energies modeled here) is significantly higher for the winter cases 
than the summer cases. 

Figure 10 shows the modeled volume emission rate altitude 
profiles for the solar minimum and solar maximum summer cases. 
As the atmosphere expands under the influence of increased solar 
activity, the production rate for a given energy peaks at higher 
altitudes. note is the fact that the behavior of the production 
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Fig. 5. LBH dependence on |Oj. [CKJ. and(N : ]. Each curve shows the LBH intensity for the doubled constituent case divided bv 
the LBH intensity from the standard (unchanged) case. 
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Fig. 6 . Relative altitude density profile variations as a function of solar 
activity and seasonal variation. 



Fig. 7 Sensitivity of the OI 1356/LBH 1838 intensity ratio to solar 
activity and seasonal variation. Solid lines indicate summer conditions, 
dashed lines are for winter. Within each seasonal set of curves, the largest 
intensity ratios are obtained at solar minimum activity, the smallest ratios 
occur at solar maximum. 


rate profiles for electron energies below 2 keV is markedly dif- 
ferent from that above 2 keV. The altitude at which this energy 
loss peaks is approximately 140 km. This is the altitude below 
which Ot becomes a competitive constituent (see Figure 1). 

Figure l i illustrates the effect of local 0 2 absorption. As would 
be expected, for emissions lying outside the region of 0 2 absorp- 
tion (LBH 1838), the production rates are unchanged by local 0 2 
absorption. Foremissions within the Schumann-Runge absorption 
continuum, however, the shape of the emission rate profiles is 
changed significantly due to local 0 2 absorption. 

The final ratio modeled was LBH| ong ; LBH s h orl (Figure 12). The 
ratio of LBH 1838 to LBH 1464 shows a dependence on the 
incident electron energy that vanes only slightly with solar 
activity. As above, this variation can be explained by the relative 
densities of N 2 and 0 2 . The observed variability from solar 
minimum to solar maximum is due to changes in the 0 2 column 
density and hence in the 0 2 absorption. 

4. Discussion 

We have shown the intensity ratio OI 1356/LBHi ong to be a 
useful diagnostic for determining the charactenstic energy of the 
auroral particles using LBH )ong to be LBH 1838. The I 1356/ 
I 1838 ratio is a very sensitive indicator of characteristic energy, 
changing by a factor of 13 or more over the range 200 eV to 10 
keV, but this ratio can vary by up to factors of 2 with changes in 
the neutral atmosphere. Almost all the change is due to variations 
in I 1356. In addition, the I 1838/1 1464 ratio shown in Figure 12 







'"30 


Germany f.t al Auroral Emission Dependence on Composition 




Fig. 8. Dependence of LBH 1838 and Of 1356 on solar activity and seasonal variation. The curves have the same interpretation as 
in Figure 7 
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Fig. 9 Variation of altitude density profiles as a function of solar cyclic 
and seasonal variations. The profiles have been normalized to the standard 
case (summer. F (0 7 — 200). 


is found to be another potentially useful determinant of 
characteristic energy. The N 2 LBH long-to-short wavelength ratio 
vanes by about a factor of 3 between 0.2 and 10 keV but, unlike 
the OI 1356 emission, is almost insensitive to changing atmos- 
pheric composition. The longer wavelength LBH bands are also 
useful indicators of the total energy influx, while the long-to-short 
wavelength LBH intensity ratio provides information on the 0 2 
Other workers have also attempted to use optical emissions to 
charactenze the auroral electron prectpitatio or example, 
Ishimoto et al. [1988] studied a similar ratio (OI 1j56/LBH 1928) 
using a recent version of the Strickland model. Their modeled 
emission ratio shows a sensitivity of 9 between 1 ana 10 keV, in 
good agreement with our results. In their description of the auroral 
code used in this study. Richards and Torr [1990] conduct a 
comparison between the two-stream model used here and the 
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Fig 10. Emission production rate altitude profiles tor Ol 1356 and LBH 


more sophisticated multi-stream model of SJW and find good 
agreement in the shape of the OI 1356 emission curve. (Dif- 
ferences in magnitude are due to the use of revised Ol cross 
sections in our model.) In addition to these studies. Rees and 
Lummerzheim ( 1989) have suggested the N? 3371/N : * 4278 emis- 
sion ratio as a determinant of the incident auroral energy. Their 
results, however, disagree with calculations by SJW and with our 
model, which shows the 3371 emission to be independent of 
characteristic enerev above 0.5 keV {Richards and Torr. 19901. 


We have investigated the sensitivity of OI 1356 A. LBH 1464 
A. and LBH 1838 A auroral emissions to changes in the neutral 
atmosphere. Our studies show that OI 1356 vanes linearly with 
f O) to within 209c and shows much less variation with other 
atmospheric constituents. The LBH 1838 A intensity is relatively 
insensitive to typical uncertainties in the neutral atmosphere (fac- 
tor of 2 at auroral altitudes). 1 1464 shows larger variations 
because of its additional interaction with O?. Our results are in 
aood acreement with similar sensitivity studies performed by SJW 



r.hKMANY hi \l Auroral Emission Dependence on Composition 




Production Rot* (cm'* j") 



Production Rot* (cm' J »**) 


Fig. I ! Same as Figure 10 (solar maximum!, but illustrating the effects 
of local Q : absorption The dashed curves are with local 0 ? absorption. 


who used a Jacchia model atmosphere [Jacchia . 1977] to model 
Ol 1356 dependence on [O! The dependence of these intensities 
on much larger changes in the composition such as might be 
encountered over seasonal or solar cycle extremes has also been 
investigated. It is found that the OI 1356 A intensity is sensitive to 
compositional changes while :e N? LBH long wavelength emis- 
sion is relatively insensiti\ ■ such changes. 
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Fig. 12 Sensitivity of LBH| onF /LBH, hof1 intensity ratio to solar activity 
and seasonal variation. The curves have the same interpretation as in 
Figure 7 
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The capability has been developed lo model thermosphenc airglow emissions on a semiglobal scale \L < 5). 
This model produces volume emission rates as a function of altitude, latitude, longitude, and local time tor anv 
selected date, and solar and magnetic conditions. The model can thus be used to provide three-dimensional maps 
ot ihe selected emission tor comparison with data obtained trom orbiting vehicles. As such it becomes an 
essential tool in the planning and interpretation ot airglow observations. A unique teature ot the model is that it 
incorporates lull interhemispheric coupling by solving all the appropriate coupled equations along the magnetic 
tlux lubes Irom the mesosphere in one hemisphere to the mesosphere in the other hemisphere. As a result the 
effects ot conjugate photoelcctrons land heat fluxes) can be fully explored. In this paper we select two thermo- 
spheric emissions with which to demonstrate the capability. The first is the 7320-A emission Irom the metastable 
0~( : P) The second is the permitted emission at 3371 A from the N : second postive 0-0 band. These two 
emissions, lor which the photochemistry is relatively well understood, are used to show the seasonal, diurnal and 
solar cyclic variations on a scale that covers mid- and low-latitudes, and the ettects of interhemispheric coupling 
(conjugate photoelectrons ». 


Introduction 

Airglow emissions are important indicators ot atmospheric 
composition and the mechanisms responsible tor the production 
and loss of the particular excited state trom which the airglow is 
radiated. For example, the 0-4) band ot the N; second positive 
svstem. which radiates at 3371 A. is excited in the airglow' by 
photoelectron impact, and ts lost only by radiation As a result, 
this emission is an excellent indicator ot the photoelectron excita- 
tion rate \Kopp ct al.. I977|. The O ( P) metastable state, which 
radiates at 7320 A. is excited both by photoelcctrons and by 
extreme ultraviolet photons. However, because it is a long-lived 
state, it is lost bv quenching by O. NU. and electrons in addition to 
radiation Thus at high altitudes, where radiation is the dominant 
loss mechanism, the 7320-A emission can be used to inter either 
the atomic oxvcen concentration or the solar ultraviolet flux it the 
other ts known [Meriwether et a!.. 1978; Rusch et al .. 19761. 

In the past, a number of detailed studies ot these emissions have 
been made using measurements of the surface brightness altitude 
profiles obtained from orbiting spacecraft (see. for example. 
Walker et al. 1 1 9751 and Rusch et al. j I977|>. In this study we 
have globally modeled these two emissions tor a variety ot 
conditions. 

The model we have used here is one that we have steadily 
developed over the years. We solve the coupled time dependent 
energy, momentum, continuity, and photoelectron transport 
equations from 80 km in one hemisphere, along a field line to 80 
km in the other hemisphere (Figure h {Young et al .. 1980 a.b: 
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Richards and Torr. 1985«. 19881. The equations that are solved 
arc summarized below: 

1 . Ion continuity equation for major ions is given by 


= Q, — L,N, — V*tb f 
dr 


(I) 


where A is the concentration of the /th major ton. Q ( and L , are its 
production and loss frequency, respectively, and is the ion tlux 
defined below. The electron density is assumed to be equal to the 
sum of the ion densities. 

2. Momentum equation is given by 

4), = N,U, (2) 


where 


.^ + J_vr, + 2LZA\v r 
at, t, "<■ 


‘'■-(szh-H*™' 

<3> 


where U„ is the neutral wind velocity and where the subscripts / 
and j are applied first to 0 + and H * respectively, and then 
to»H + and He* respectively. O*. H + , and He* are coupled 
through collisions and the polarization electric fields. Since the 
influence of 0 + and He* on H + becomes significant in dif- 
ferent altitude regimes and since the effect ot He + and O is small 
due to their mass ratio, this decoupling ot a system of three major 
ions into 2 pairs of major ions considerably simplifies the numer- 
ical calculation of the major ton densities. Here a f/ and a , / are 
thermal diffusion coefficients and D, is the ordinary ion diffusion 
coefficient of species /. as in the work by St. Maurice and Schunk 
1 1 977 1 . 
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3. The thermal electron energy equation is 


and 


i MQc 'lL = - N,MT e X-U r - ^ .V,J W, ?T, - V-I ir + IQ, - lL r 

4 The ion energy equation is 
l ,V,* 'lL = - H,kT,V U , - 1 \,kU,^T, - V £/, + IQ, - 1L, 

^ * lit ‘ ' *> ^ 


3 


(4) 


5. The ion heat flow equation is 


q , = -L (-K.VT,- 

1-4 \ N/’W ' 7 


(5) 


( 6 ) 


6. The electron heat flow equation is 

q e = -A..VT, (7) 

The thermal conductivity coefficients A, and A,., and the term 

where v,,' , u,/. v,' , and v/ arc the effective collision frequencies, 
are given by St. Maurice and Schunk 1 19771. The ion-neutral 
collision frequencies v ltl are from Schunk and Na^y |I980|. 

7. The photoclcciron Liouville equation is given by 


R d <t>* 

B JI¥ 


"" 7^2^ ^ f 1^ 4- 


2<cos 0> <co.s 0>' 


b £j 


2<cos 0> <cos 0> 


( 8 ) 

(9) 


where 


<t>UE.s) 

4T(£.i> 

q{E,s) 

q- 


<cos 4>> 
B 
Ti 


photoelectron flux outward along s: 
photoelectron flux inward along s ; 
photoelectron production rate in the range E to 
E + dE due to direct ionization processes; 
photoclcciron production in the range E to E + dE 
due to cascading from higher energy photoelcctrons 
undergoing inelastic collisions; 
average cosine of pitch angle, 
magnetic field strength: 

k 

T? = X #i t |a/+/» r Vl; 

k 



Fig. 1 . Illustration of the mterhemisphenc nature of the code in which the 
coupled and time dependent equations are solved from 80 km in one hemi- 
sphere. along a field line, to 80 km in the conjugate hemisphere. 


// , At h species number density; 

photoelectron backscatter probability for elastic con- 
ditions with the Ath species: 

t< k photoelectron total scattering cross section for elastic 
conditions with the Ath species; 
r r / inelastic cross section for excitation of the Ath pan- 
icle species. 


The model includes an option to increase the 0"~0 collision 
frequency as recommended by Burnside vi al. [1987]. The use of 
the Schunk and Nayy |I980| values here does not significantly 
influence the results ot this paper. 

The above formulation corresponds to conditions where the dif- 
ferences between species temperatures and flow velocities are 
assumed to be small, i.e.. stress and nonlinear acceleration terms 
are neglected. In addition, density and temperature gradients 
normal to the geomagnetic field lines arc neglected and we 
assume that the electron and ion temperature distributions are 
isotropic. 

The continuity equation is solved using a rather unique 
approach. We can rewrite ( I) in terms of a function E: 

dF = — + V ■ ( N,U,)-Q, + L,N , (10) 

Bl 


Then, using a Newton iterative procedure to find the minimum of 
F, we solve for the density at the grid point, j. Figure 2 shows 
how the field line is divided into elements about the actual grid 
point j. The lower limit of the element (f) is placed midway 
between the grid points j and y — I , and the upper lo ot (u) is pla- 
ced midway between y and j + I . The tower limn of one element is 
the upper limit of the preceding element. We then integrate (10) 
between u and C, and obtain the densities by solving 



0 (II) 


where B is the amplitude of the magnetic field The values of tiNl 
fit. Q , and L arc obtained at the limits of integration, u and J. by 
interpolation between the actual grid points. 

In the past, many models have encountered numerical problems 
above about 3000 km. due to the large diffusion coefficient in this 
region which results in small density changes producing large 
changes in velocity. At lower altitudes, both ion-ion and lon- 
neutral collisions are important, while at greater altitudes, colli- 
sions become less important and the plasma can be accurately 



Fig. 2. The numerical grid scheme used for the solution ot the continuity 
equation. 
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described using a diffusive equilibrium approach. The solution ot 
<11) distinguishes this model from previous methods which eval- 
uate the terms of the integral at only one point. As can be seen 
trom ill) and Figure 2. the flux at the upper limit ot one element 
becomes the flux at the lower limit ot the next element. These 
fluxes must be identical and hence this method has been called the 
'flux preserving scheme. “ Furthermore, the flux at any grid point 
is closely tied to the flux at neighboring grid points, allowing 
stable solutions even for regions where large changes would give 
rise to unstable solutions with other numerical methods. 

The model starts at noon with ‘best guess ‘ initial values. It is 
allowed to run for 12 hours tin local time! before results arc used 
in order to reduce dependence on initial conditions. It has been 
found that the ionospheric densities over a midnight to midnight 
diurnal cycle typically repeat with only small differences due. tor 
example, to plasmaspheric refilling. The plasmasphcnc H" and 
He * contents arc initially low. and the flux tubes are allowed to till 
continuously. The transport equations for the three major tons. 
He * . H ~ . and O". arc carried out in two steps. As 
mentioned above, the coupied O" and H~ equations are solved 
tirst. followed by the He* and H~ equations, where the latter 
use the O " -H * results. Thus the equations tor the three major 
ions are essentially solved in a simultaneous manner. Below 
approximately 1X0 km. NO" and O?" become major ions, but 
arc obtained from photochemical equilibrium calculations 

The numerical solution of these equations, boundary conditions 
and other details are further discussed by Yonntt et al. \ 1980a,/>|. 

The full interhcmisphcric coupling is difficult to handle, but 
once incorporated it imposes no artificial upper boundary con- 
ditions for both thermal and photoclcctron fluxes. This is most 
important for thermal coupling and the proper treatment ot con- 
jugate photoclcctrons [Richards and Ton . 1985/?| w hich can be 
significant in the calculation of airglow emissions. Typically, 
attenuation of the conjugate photoclcctron flux by Coulomb colli- 
sions results in approximately 5C energy loss The model 
includes the option to specify loss duo to pitch angle scattering 
which we assumed to be zero for this paper. A lilted dipole 
approximation is used for the Earth's magnetic tield | Rulum/s 
and Torr. I986a|. 

The concentrations of the major neutrai species are provided by 
the MSIS-S6 j Hedin. 1 987 1 to the model which (hen computes the 
concentrations of minor and excited slate species and major ions 
j Torr. 1985: Richards ci al., 1982a. 1986/fl. In this paper a 
simple model giving daytime poleward winds and equatorward 
neutral winds was employed. This behavior is consistent with the 
results obtained using the method of Richards and Torr | I986/?| 
and Miller el id. (I986|. The model includes the option to use 
model winds of Hedin el al. |I988|. Killeen el al. |I987| and 
Killeen (private communication 1989) Note, if the Hedin et al. 
j I988| model is used, the O " -O collision frequency of Burnside cl 
al j 1 987 1 should be used in order to produce the observed h,„F2 at 
night. 

The chemistry' of all significant emitting species is incorporated 
in detail, including the excitation of the metastable states | M. R 
Torr and Torr. 1982] and vibrational states [Richards et al., 
1986a; Richards and Torr. 1986c|, and the odd nitrogen 
chemistry' [Richards et al., 1981. 1982/?; Richards. I986j. The 
vibrational population distributions ot 6N are determined, an 
important factor in calculating the ionospheric O" and ;V ( 
concentrations. The model also includes calculation of the vibra- 
tional populations of N : ^ . but this does not significantly affect the 
results reported here. The photochemistry is that described by 


Torr ( 1 985 j as updated and shown here in Table 1 and illustrated 
m Figure 3. Al equatorial latitudes, where electric fields play an 
important role, the electron concentrations are obtained from the 
fully analytical ionospheric model of Anderson et al. |I989|. 
Hlsewhere. the electron densities are computed self-consistentlv 
by this model. The transition occurs between L - 1.5 and 1.8 
w ith interpolation between these L shells. The major elements of 
the code arc showrn in Figures 4a and 4 h. The solar EUV flux is 
obtained in the following wav. For solar minimum (FI0.7 = 71) 
the model utilizes the F74I 13 reference spectrum from Torr et al. 

[ 1 979 1 , with the fluxes below 250 A doubled as recommended by 
Richards and Torr | 1984} and supported by Ottawa and Judtte 
[1986|. For other levels of solar activity, each of the 37 
wavelength intervals is scaled linearly as a function of FI 0.7 
using the solar maximum measured flux at FI0.7 = 206 given by 
Torr and Torr 1 19851. For further details, see Richards and Torr 
| 1 988 1 . which also provides the cross sections used. 

The main outputs of the model include ion densities (O". 
0~r*S). O '(->). H". He". N". NO^. N^. 

N\' M. neutral densities (N( 4 5). N (~D). N i‘P). NO. 0(‘D). 
Ol’.S). N-{A'22‘). NyM. electron and ion temperatures and 
flow velocities, the photoclcctron flux, and a large number ot 
emissions (see Figure 4a). 

Dunne the Atmosphere Explorer C. D. and E program, 
numerous studies were conducted which compared the photo- 
chemistry of the code with in situ measurements of species con- 
centrations. The photochemistry yields results consistent with the 
data base taken over the lifetime of the AE satellites. Generally, 
verv good agreement with measurements has been obtained with 
re card to all parameters with the exception of high altitude 
electron temperatures in the plasmaspherc and ionosphere [New- 
berry cr al.. 1 989 1 . The model has been extensively tested against 
comprehensive satellite and incoherent scatter radar data bases 
j Yountt cr al.. 1980a./?: Chandler et al.. 1983: Richards and 
Torr. 1988: Newberry et al.. 1989: Horwitz el al.. 1990: Richards 
et al.. 1 989 1 Apart from the input parameters (such as data and 
location) the only tree parameter in the code is the pitch angle 
mattering of photoelectrons m the plasmaspherc. In its present 
torm. the model is ideally suited tor studies ot the airglow 
emissions. 

The code is run on the Marshall Space Flight Center CRAY 
XMP computer. Values of output parameters arc provided on a 
elobal end of points, providing results in a latitude, longitude, 
altitude, and local time mesh for any selected dale, or solar or 
magnetic conditions. For the cases discussed in this paper, we 
have run the model for 144 flux tubes, which corresponds to 144 
northern and 144 southern hemisphere locations Figure 5 shows 
the locations of the field lines along which the equations arc 
solved. These are constrained to L ^ 5. 

We have chosen the 0"( : P) emission at 7320 A and the N : 
second positive 0-0 band emission at 3371 A tor the initial global 
modeling. The calculations have been made tor November 28. 
1983 for which the FI 0.7 cm flux was 89. and the Ap index was 
23. Thus the November 1983 calculations correspond to a period 
of relatively low' solar activity. An earlier example of the results 
for the 7320 A case has been shown by Torr et al. \ I990|. The 
results shown here represent a significant improvement over the 
Torr et al. \ I99()| case, in that we have added approximately 50 
more flux tubes at low latitudes ( 144 versus 96). In addition we 
have used much smaller time steps through the twilight conditions 
(5 minutes versus 20 minutes). For comparison, wc have also run 
(he calculations for the same day of year, but tor conditions 
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TABLE 1 Summary of Photochemistry Used in the interhemisphenc Model 


Reaction 

Number 


Reaction 


Rate Coefficient (cnr s _( ) 
or Rate (s _l ) 


Reference 


-4 x I0“* - (77300)° 


Torr [ 1 985] 


1 

o" 

+ r -40 + hv 

-> 

or 

4 r -40 + o 

3. 

cr 

+■ o : +or 4 o 

4 

O" 

4 N ? h»NO " + N 


5. 

NO" + e + N 4 O 

6. 

Nr + O+NO" 4 N 

7. 

Nr 4 e’-^N 4- N 

8. 

Nr 4 o>0" 4 n 2 

9 

Nr 4 0,-40," 4 n. 

10. 

N" 4 0 } 40r 4 N 

1 1 

N" 4 0.-4NO* + O 

12. 

N" 4 0 : -4 NO" 4 o CO) 

13. 

or 4 N -4 NO * 4 o 

14. 

or 4 NO+NO" - o. 

15. 

0"( 2 D) 4 N,-4N : " 4 o 

16. 

0"( 2 D) 4 0-40"( 4 S) 4 0 

17. 

o"( 2 d) 4 0,-40," 4 o 

18. 

O "( 2 D) 4- r-40"( 4 S) 4 e~ 

19 

0"( 2 P) -40 "( 2 D) 4 hv 

20. 

0"( 2 P) -*-0 " ( 4 5) 4 hv 

21. 

O "( 2 P) 4 f>0 "( 2 D) 4 e~ 

22. 

O "< 2 P) 4 e+O "( 2 D) 4 e~ 

23. 

O *( 2 P) 4 N,- 4 products 

24. 

0"( 2 P) 4 0-4 products 

25. 

He 4 b-fHe' e~ 

26. 

He" 4 N,-4N" *- N 4 He 

27. 

He" 4 N : 4Nr - He 

28 

O" 4 H -4 H " 4 () 

29. 

H" 4 0-40" 4 H 

30. 

O" 4 N( 2 D)-4N" + O 

31. 

NO" 4 e' + N C-D) 4 O 

32. 

Nr 4 f>N ( : D) + N 

33. 

Nr 4 0-4N( 2 D) 4 NO" 

34. 

N" 4 o 3 - 4N( 2 D) 4 or 

35. 

N( 2 D) 4 0-4N( 4 S) 4 0 

36. 

N( 2 D) 4 0,-4 NO 4 0 

37. 

N( 2 D) 4 e‘-4N( 4 5) 4 e~ 

38. 

N( 2 D) 4 O/ 4 NO" 4 O 

39. 

N( 4 5) f 0,-4 NO 4 O 

40. 

N 4 NO-4N, 4 O 

41. 

O, 4 h v -4 0( 1 D ) 4 0 

42. 

or 4 f'40(‘D} 4 0 

43. 

O i'D) 4 N ; -40 < 3 P) 4 N, 

44. 

O CD) 4 0,-40 OP) 4 O, 

45. 

0(‘D)-4CH 3 P) 4 hv 

46. 

N ( 2 D) 4 0,-40('D) 4 NO 

47. 

or 4 r -40 VS) 4 0 

48. 

or 4 n^o vs) 4 no- 

49. 

O i'S) 4 0( 3 P) - 4 O 4 0 

50. 

O CD) 4 0-40 4 O 


1.6 x I (T 7 (300/7*,)° ” for T e ^ 1200 K 
2 x 10- 7 (300/7,)° 7 for 7, < 1 200 K 

2. 1 x 10*' 1 {7„ 4 2 7/3 x 300}“° 763 

1 533 X 10- 42 - 5.92 x 1CT* ' (7 cr /300) 

4 8.60 x I0" 14 (7 eff <300) 2 
for 300 ^ 7 eff ^ 1700 K 
2.73 x |0"‘ 2 - 1.155X I0" ,: (7 <# /300) 
+ 1.483X lO" 1 ' {7 ef y300) 2 
for 1700 < T eif < 6000 K 

4.3 x 10“ 7 (7/300r‘ 

1 .4 x I0“‘° (T^OOr 44 for T, < 1500 K 

2.7 x |(T 7 

0.07 kt, (7/300)° 21 

9. 1 x 10"' ‘ exp ( -0.002 7 eft ) 

4 x 10"'° 

2 x I O' 10 

3 = 0.7 

1.2 x I0“*° 

4.4 x I0“‘° 

8 x 10-'° 

5 x kt 12 

7 x 10*'° 

6.6 x 10"* (300/7,) 5 
0.173 s" 

0 047 s" 

1 .5 X I (T 7 (300/7,) ■ 

4.7 x 10^ ■ (300/7,) 3 

4.8 x I0“'° 

5.2 x 10-" 

4 Ox I0~* - 1.2 x 10 7 <c> 

1 x l(T 

6.5 x IO- tu 

2.2 x 10"" (T,) s 

2.5 x 10"" ■ <T„) s 

5 x 10"" 

3 • kt, where (3 = 0 76 
(3 /r 7 where (3=19 

(3 • U, where [3=10 
[3 Ac| 0 where (3=10 
-7 x I0"' 3 

6 x ICT 12 

5 x l(T i0 (7/300) 5 

1 x 10"* 

4.4 x 10'* 3 exp ( -3220/7) 

3.4 x I0‘" 

3 = 1; J x (0,)sr = (1.5- 2.8) x 10^ 
3 Jt, where 3 = 1.2 
2. OX 10*" exp (107.8 T„) 

2.9 x 10*" exp (67 5/7J 
0.00934 

—5 x I0* 12 

3 k 2 where 3 = 0.08 
-2x 10*" 

2 x I0“ 14 

8 x lO"' 2 


Torr and Torr (1981): 

Mehr and Biondi ( 1 969] 
Chen et al. [ 1978]*' 

St. Maurtce and Torr j 1978): 
Albritton (1978]: 

Chen et al. (1978] 


Torr and Torr ( 1979] 

McFarland et al. [ 1974]; Torr ( 1979] 

Abdou et al. ( 1984] 

McFarland et al. (1974] 

Lindinger et al. ( 1974] 

Huntress and Anrcrch ( 19761 
Huntress and Anicich ( 1976] 

Langford et al. ( 1985] 

Fehsenfeld (1977] 

Lindinger et al. ( 1974] 

Rowe et al. ( 19801; 

Johnsen and Biondi ( 1980] 

Abdou et al. [ 1984] 

Johnsen and Biondi (1980] 

Henry et al. ( 19691 
Seaton and Osterbrock [1957] 

Seaton and Osterbrock ( 1 957] 

Henrx et al. [ 1969] 

Henr\ et al. [ 1969] 

Rusch et al. (1977] 

Rusch et al. [ 1977] 

Torr and Torr ( 19851 

Adams and Smith [19761 

Adams and Smith ( 1976] 

derived from Banks and Kockarts ( 19731 

derived from Banks and Kockarts { 1973] 

Torr et al. ( 1979] 

Kiev et at. [\911] 

Queffelec et al. ( 1985] 

Frederick and Rusch j 1977] 
assumed 

Richards et al. [1981] 

Lin and Kaufman |197l] 

Frederick and Rusch ( 1977] 

Daigarno ( 1970] 

Becker et al. ( 1969] 

Lee et at. (19781 
Torr et al. ( 1980] 

Abreu et al. [ 1986] 

Streit et al. (1976] 

Streit et al. (1976] 

Fischer and Saha ( 1983] 

Rusch etal. (1978]; D. G. Torretal. [1981] 
Bates and Z/p/11981 ]; Abreu et a I. (19861 
Frederick et al. [1976] 

Slanger and Black [1981] 

Abreu et al. (1986] 
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TABLE I (continued) 


Reaction 

Number 


Reaction 

0( '5) + 0( ] D +■ hv 
+ O (V*) +■ hv 
N ( : D) + NO + N > i- O 
0 ( ' 5 ) + 0 : + 0i'P) + 0 > 
N : (A^,r) + O + products 
SdA^,;} + O+O < l S) r : 

+ 0+0" + Ni 


O + hv^O~(*S) +■ e 
+ 0"{‘D) + e 
+ 0 "f 2 P) + e 
+ 0"( 4 J p ) + e 
+ 0*( 2 P") + e 


Rale Coefficient tcm s" 1 ) 
or Rate is -1 ) 


3 ■ Ac <; 4 where [3 = 0 37 

— 2 x lQ- iOS 

(1.78 - 5 14) x KT'* 

(3 06 - 8.82) x 10 7 
(0.98 - 2.81) x I(T 7 
(0.79 - 2.34) x 10" 7 
(0.45 - ! 38) x I0" 7 
(1.04 - 3.43) x 10“* 

(0 46 - 1 .42) x 10“** 


Reference 


Kernahan and Pan g 11975] 

Un and Kaufman 1 1971] 

Zip/ (1979] 

Piper (1982) 

Piper 1 1982] 

Abdou et al. [1984]; Torr (1985] 
Torr and Torr j 1985]‘ 

Torr and Torr ( I985)‘ 


“Since the results ol Chen eial. 1 19781 stop at 700°K. »e normalize the converted drift tube data parameterized bv Sr Maurice and Torr ( 1 9781 at this 

"The model computes this rate coefficient as a function ot N- * vibrational temperature i7l. The value reduces to that given bv equation (bl when T, = T„. 
Inclusion of this process is an option available [see Abdou et al.. 1984], 

The ranges given for the ionization frequencies indicate the variation over a solar cycle. 



Fig. 3. Schematic of the thermospheric and ionospheric chemistry that is incorporated in the model. 
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corresponding to solar maximum \h 10 7 = 194. Ap = 23). This 
study deals only with the airgiow. and we do not show results tor 
latitudes higher than L = 5. 

Details ot the resuits are given in the following section. 

Mode-ling of the 7320-A Airgi.ow Emission 

The 7320-A emission arises from the 0*( P) state, which is 
produced in the thermosphere by two mechanisms; photoiomza- 
tion and photoeiectron ionization: 

O ■+■ fw ->-0 ~ (~P) c 
O + e* + 2e 

where the photoelectrons may be those produced locally, or those 
produced in the conjugate hemisphere and transported along the 
magnetic Meld lines. The O'D/ 3 ) is lost by several mechanisms: 

idiation 

O ~( : P)+ O VS) + hv ^ 2 o v 
collisional deactivation 

0*i 2 P) 4 - N : - O 

Cr< : P> + O -> O * * O 
CT( : P> - e ^G*( : D) + e 
CT< : P) + e + OVS) + *' 


This photochemistry has been reviewed bv Torr and Torr { I982|. 
It should be noted that since publication of the acronomicaily 
derived values ot the rate coefficients by Rusch et al. [1977]. the 
solar EL Y Mux below ^2^0 A was doubled. We estimate that 
when quenching dominates, the reported intensities may be about 
30# too high. 

Plate I shows the 7320-A volume emission rate at the peak of 
the layer as a function ot latitude, local time (longitude) and alti- 
tude. The upper plot shows the peak volume emission rate as a 
function o: '.ititude and local time. Because the model has long- 
itudinal variability resulting from both the MSIS | Hedin 1987) 
model atmosphere, and the mterhemispheric coupling, this par- 
ticular plot is shown for (X) UT. This UT is equivalent to placing 
midnight at 0° longitude, the Greenwich meridian, and noon at 
180° longitude, i.c, over the Pacitic Ocean. Where the volume 
emission rate becomes so small as to be effectively zero, the 
values arc not plotted. These regions can be seen near midnight 
for equatorial latitudes, and represent no production at these 
times. The graphics tend to smear the nonhem and southern lati- 
tude boundaries by a tew degrees to the north and south, respec- 
tively. The solutions arc only valid, however, within the L ^ 5 
region illustrated in Figure 5. 

The ba e features shown in Plate I have been discussed by 
Torr ct a! , j 1990| but we will review them here briefly as this 
plate w ill be used for the comparison with other cases. Summer is 



(a) 



Fig 4i</>: Schematic illustration ot the major input and output elements of the held line mterhemispheric plasma (FLIP) code. 

Illustration ot the Mow ot ihe solutions 
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Fig. 4. (continued > 
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about 23' ; Hie maior hieii altitude loss mechanism is radiation, 
winie quenchinn b\ V dominates at lower aitttudes. It is interest- 
ing to note in Plate 1 that the midday peak volume emission rates 
m the winter i northern ) hemisphere are somewhat larger than in 
die summer hemisphere. I he reason tor this can he seen in Fissure 
fv While the V peak quenching rate is aimosi identical in both 
hemispheres, the V quenching tails ott with altitude with a 
smaller scale heudn in the winter hemisphere, so that radiation 
dominates io a lower altitude with the resulting higher emission 
rate. 

At and bevond the terminators i tor solar zenith angles larger 
than dm sunlight illuminates increasingly higher altitudes, and 
the peak altitude rises white the volume emission rate fails 
sharply. In Future 7 we show the production and loss rates for a 
solar zenith aneie of 105 3 in the evening southern hemisphere. 
The peak LA' photoionization has risen to above 6(K) km. A 
second peak is formed near 300 km. The latter is a result of photo- 
electrons which arc produced in the ionization process that caused 
the upper peak. The photoelectrons arc transported downward and 
lose their encruv near 300 km in the impact ionization of atomic 
oxygen. Plate 2 shows the same information as was given in Plate 
1 . but versus solar zenith angle instead ot local time. Phis plate 
shows ihe extent to which seasonal hemispheric asymmetries are 
reduced when the data are plotted m [his formal. 



Plate I The 7320* A \ oiuine emission rate (photons cm" ‘ s' 1 > at the peak ol the layer as a I unction ol altitude, latitude, and loeai time 
loneitudc i lor DO [ 'T The upper plot shows the peak volume emission rate only as a function ol I at nude and local time I he volume 
amission rate is on a ioe s v alc. Phe results are only valid within the L '-■■■■ 3 boundaries shown in Figure .v 
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An interesting feature of Plate l is the fact that there is signifi- 
cant production (-0.2 cm' 3 s' 1 ) at high altitudes at midnight in 
the winter hemisphere. This is due entirely to photoelectrons 
transported from the summer hemisphere. This conjugate electron 
production rate together with the loss processes are shown in 


Fieure 8. Without the intcrhcrmspheric capability this production 
source would not be modeled. 

Because of the longitudinal asymmetries of the Earth s mag- 
netic field, the conjugate photoelcctron production rate varies for 
different maps generated for different universal times. Plate 3 
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Fig. 7. Production and loss rates tor 54°S tor a solar zenith angle ot 105 (LT IK)). 


shows the global plot for the same conditions as Plate l . but for 
LT = 0600. This is equivalent to placing midnight at 90°W. over 
the eastern United States. In this longitudinal sector, higher 
southern latitudes (solar illuminated) map to lower northern lati- 
tudes. Above 50°S the Sun does not set and conjugate electrons 


arc present all through the night near 50°N. Thus the conjugate 
source is seen to be significantly more pronounced than tor 00 
LT. Plate 4 shows the integrated volume emission rate, or surface 
brightness un ravicighs). with and without the conjugate photo- 
ciectron source. The scale has been adjusted u nhancc the winter 
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Piute -in With conjueuie photoeleciron production. 



Piute 4 V crticui 7120- A lritcn’Ots iravieiizhsf lor soiar minimum ut LT - Ob The solution is onh valid within ihe L * "'boundaries 
'iiown in Injure 5 
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Plate 5. The 7320- A volume emission rate t photons cm ’ ' s" \ at the peak of the laser as a tunction of altitude, latitude . and local time 
i or IX) IT. The results are lor the same dav as Plate l i Novcmocr 28) but lor solar maximum conditions. The volume emission rate is 
shown on a ioa scale 


nienmme conjugate production. This interhcmisphertc source 
results in 2 to 3 R of 7320- A airglow 

Plate 5 shows the same day of year as Plate 1 . but for conditions 
representing solar maximum. The principal differences between 
these two are in the overall increase in emission rate and the 
overall increase in altitude of the emission peak. Both ot these 
effects are to be expected on the basis of the solar flux and neutral 
atmosphere changes. 

Plate b shows the results of integrating the volume emission 
rates ot the solar minimum and maximum cases (Plates 1 and 5) 
over altitude to obtain the vertical surface brightness. These 
intensities (in raylcighs) arc what an orbiting instrument would 
see looking vertically downward on the atmosphere. 

Modeling of the N : 337 1 -A Airglow 

The production mechanism for the excitation ot the N : second 
positive system in the dayglow' is photoelcctron impact: 

ViX'l,. ' ) * v N>(( "tt ,i - c 

E he 337 I - A emission arises from the iMi band of the second posi- 
tive s\ stem: 

-r ■ - 'O' 


Since the transition is permitted, radiation is the only loss process. 
The excitation cross sections for the 0-0 band are those ot hmimi 
and Borst [1974| which are 25* of the total cross section. 

Plate 7 shows the results of the semiglobal solution ot the 337 1 - 
A volume emission rates. The production rates tor noon at mid- 
latitudes arc shown in Figure 9. The results shown in Plate 7 
follow the behavior anticipated from a simple photoelcctron 
source and radiative loss. The behavior changes alter sunset in 
that the peak production rises in altitude and tails in magnitude. 
A earn, the high latitude conjugate photoelcctron source can be 
seen during the night in the winter hemisphere. While we do not 
show the 06 UT ease here, the conjugate source again becomes 
more pronounced for the different magnetic field orientation, as it 
did for the 7320-A case, because of the more favorable alignment 
of the magnetic field lines tor this purpose. 

Plate 8 shows the global model of the peak 337 1 -A volume 
emission rate for solar maximum. What is interesting to note in 
this case is that, apart from an overall increase in the altitude of 
the emission peak, the peak emission rate distribution is very 
similar at solar maximum and solar minimum. However, the dif- 
ference in the vertical intensity is significant. 1 his is seen in Plate 
4 which shows the integrated column surface brightness for solar 
minimum and maximum. The difference is largely due to the dif- 
fering scale heiehts. as can be seen in Figure 9 
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Plate bo. Solar mintmum case shown in Plate I. 




Plate tb. Solar maximum case shown in Plate 5 


Plate 6 \ ertical column integrated surface brightness t m ravleighsi. The results are only valid w ithin the L > boundaries show n in 

I njure > 
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! late i \ ntumc emission rale (photons cm * s ' 1 1 at ihe peak ol the 337 1 - A layer v ersus attitude, latitude, and local time < longitude) 
nr near soiar minimum conditions The volume emission rate is plotted on a log scaie The results are on;\ valid within the L * s 
^‘lindanes shown in Future 3 
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In this paper we have presented the first mterhem i spherically 
coupled models ol the and Na'O airgiow emission at 

OO A and 337! A. respectively, covering mid and low latitudes. 
These are two emissions tor which the mechanisms arc relatively 
well understood, allowing us to examine here the moronoioey. 
The semi global graphical maps clearly illustrate the diurnal, 
seasonal and solar cyclic variability ot these emissions in both 
cases, conjugate photoelectrons result in excitation at mid to meh 
nighttime latitudes in the winter hemisphere. The m* tu- 

spnenc coupling capability ot the model allows us to he 


extent and the longitudinal variabihtv ol the comugate photoclec- 
tion e fleets. 

This ; odd prov ides a capability w hich vv lii greatly enhance the 
ability to interpret airgiow observations made from spacecraft. 
Any line of sight geometry can now be protected through the 
three-dimensional model solution, and bv integrating along the 
projected line tit sight, surface brightness (vertical or slant path) 
can be obtained tor comparison with the observations. Further- 
more. because ot the complex variability exhibited by these 
emissions over the L ' 5 range, this model also provides a 

valuable tool tor planning of such observations. 


Plate ^ \cmcal intensity of the 3a7!-A airgiow on ravteighs) for solar minimum and solar maximum conditions. The solar 
maximum plot is shown on two scales: one lo show the intensity variation and one to allow comparison with ihc solar minimum plot 
'•hown in Piaie The results are only valid within the *.5 boundaries shown in Hi cure 3 
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Plate 8. Volume emission rate (photons cm * s~') at the peak of the 3371 -A layer for solar maximum conditions. The volume 
emission rate is plotted on a log scale. 


Plate 'hi 
Solar minimum 
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PI ale 9/> Solar maximum, scaled to brightest intensities. 
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Plate yc. Solar maximum, same scale as soiar minimum 
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Abstract. The results are reported of the first interhemispheric 
model ot the 7320 A airglow covering mid- and low latitudes, at all 
longitudes. A comprehensive model of the ionosphere and thermo- 
sphere is used to compute volume emission rates as a function of 
altitude, latitude, longitude, and local time. Selected results are 
shown here to illustrate the computational capability. In particular 
we discuss the diumai and seasonal variability and interhemispheric 
coupling ot conjugate photoeiectrons. The model is particuiariv well 
NUited for airglow studies, and provides a valuable tool for the 
comparison, interpretation, and planning of spectroscopic observa- 
tions made Irom orbiting platforms for all significant thermospheric 
emissions 


Introduction 


The O CP) airglow emission at 7320 A provides a means of deter- 
mining the O CP) concentration. Measurements of this emission can 
be used to determine the atomic oxygen concentration or the solar 
UV ionization frequency [Men wether et aL, 1978: Rusch et al,. 
1976], The photochemistry of this species has been established 
largely on the basis of the comparison of steady-state altitude profile 
calculations w'lth surface brightness measurements made from satel- 
lites [Walker et al. . 1975; Rusch et al., 1977], Observations of limb 
brightness are inverted to yield altitude profiles of volume emission 
rates. The principal sources and sinks of O ^( 2 P) have been discussed 
in detail by Torr and Torr [ 1982]. 

Over the years, we have developed a comprehensive model of the 
ionosphere, thermosphere, and plasmasphere. This model solves the 
coupled time-dependent energy, momentum, continuity, and photo- 
electron transport equations from 80 km in one hemisphere, alone 
the field line, to 80 km in the conjugate hemisphere (Youne et al.. 
1980], The full interhemispheric capability allows for the proper 
treatment of thermal coupling [Richards and Ton. 1986] and con- 
jugate photoeiectrons [Richards and Ton, 19851, which is important 
in the calculation of airglow emissions. The concentrations of the 
major neutral species are provided as input from MSIS-86 (Hedin, 
1987] to the model which then provides the minor and excited state 
species and ions (Ton. 1985]. The chemistry of all the emitting 
species is comprehensively included [Ton. 1985). Ab initio calcu- 
lations ot the excitation and loss rates are performed for the metas- 
table species [Ton and Ton, 1982] and the vibrational states of 
molecules and ions [Richards et al., 1986], The concentrations of 
odd nitrogen species are also computed together with those of other 
minor constituents [Tonetal., 1980]. The photochemistry cunentjy 
used in the model is that defined bv Tonf 1985] {Tables A-8 through 
A- 1 2 and A- 1 9 through A-21) with the updates, conections. and 
additions summarized in Table I. 

The three-dimensional capability is achieved by running the code 
for approximately 100 magnetic flux tubes. Because the model is 
interhemispheric, solutions are obtained simultaneously for the con- 
jugate hemisphere. This yields a total of approximately 200 gnd 
points. The flux tubes are selected along various L-shells for L i 5 
and the grid on which the code is run is shown in Figure I in geo- 
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graphic space. Because of the large computational requirements, the 
MSFC CRAY II computer is used. Results are output on a four- 
dimensional global grid, compnsing latitude, longitude, altitude, 
and local time. Input parameters are date and solar and magnetic 
indices. 



Fig. 1. Location of the conjugate points along L - I 08 to L = 5 
which mark the field lines along which the computations were made. 


In this paper we have selected the 7320 A emission to illustrate the 
global airglow modeling capability. The calculations have been 
made for conditions corresponding to November 28. 1983, for which 
the FI 0.7 cm flux was 89 and the Ap index was 23. Thus these calcu- 
lations correspond to a period of relatively low solar activity. 

Results 

Figure 2 shows the results of the global solution of the 7320 A 
emission. The volume emission rate at the peak of the layer is plotted 
as a function of altitude, latitude, and local time (longitude). This 
particular piot is shown for 00 UT. which is equivalent to placing 
midnight at 0° longitude, and noon at 180° longitude. 

Several interesting features appear in the results shown in Figure 
2. Summer is in the southern hemisphere and the more extensive 
solar illumination in this hemisphere is immediately evident in the 
longer duration of the daytime peak values (volume emission rates of 
the order of 10 photons-cm ’-s' 1 ). The model results are valid for the 
latitudinal regime indicated by the envelope of the points shown in 
Figure I . The interpolation routine tends to smear the high latitude 
boundary in Figure 2 by about 3°. Figure 3 shows the production and 
loss rate profiles for noon at southern midlatitudes. The peak produc- 
tion of 0^( : P) occurs near 170 km with photoionization the major 
source. The dominant loss mechanism above 280 km is radiation, 
and below 280 km it is collisional deactivation by The combina- 
tion of these processes results in the emission peak being formed at 
approximately 260 km. 
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TABLE l. Updates to Photochemical Model of Torr 1 I985[ 


Reaction 

Rate Coefficient tenr s 1 > or s ' > 

Reference 

1 > O ' - O;- <),• - ( ) 

2.1 \ 10 

,, 1 t„ - : t, | ■" 

Chen et al. [1978] 



I * i 


2i N, - O -» O.' - \ ; 

9 1 x m 

M exp i-O 002*T CII ) 

Lindincer et al. 1 1974) 

^ O' - NrDl-N ' - () 

7 x 10 H 


Torr et al. 1 1 979| 

4i O' - H - H * O 

2.2 x 10 

M T. 5 

Banks and Kockarts ( 1973] 

5i H' - O-O’ - H 

2.5 x 10 

n T rt 5 

Banks and Kockarts [ 1973] 

am- 3“ , 




6| N, - e- - N(-D) + N 

a = 2.7 

x 10' 

Abdou et al. ( I984| 


(3=19 


Queffeiec et ai. (1985] 

?) O('D) + O, -O('P) * O. 

2.9 x 10 1 

1 1^67 < T, 

Strett et al. f 1976] 

Si Ol'D) -OOP) hi> 

\ t = i) 00934 

Fischer and Saha { 1983] 

l * * - O — products 

2 x 10 11 


Piper et al. [ 1981 } 

— Ot 'Si - S'. 

(3 = 0.37 


Piper [1982] 

IDl Nl Dl - NO — N, - () 

7 x II) 1 ' 


Lin and Kaufman 1 1971 1 

1 1 > Ot'D) 0-0 ~ O 

8 x I0' 1 ’ 


Abreu et ai. [ )986| 

I-) O(’S) - Ot 1 D) - hu 

1.07 


Kemahan and Pang [1975] 

O('P) -*■ hv 

0.0444 




“ Three-dimensional plot of the 7320 A volume emission rate at the peak of the layer versus latitude and local time, for 00 
LT. The calculations were done tor conditions appropriate tor November 28. 1983. The upper panel shows the same results, 
hut without the altitude coordinate, i.e.. peak volume emission rate versus latitude and local time. The color code is on a log 
scale, and the units are photons -cm T s -' 
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Fig. 3 Sources (ai and sinks (b) of 0'( : P) for the location 
corresponding to 48 S at 11:16 LT in Figure 2. Sources are photo- 
lonization and photoelectron impact ionization Sinks are radiation 
and quenching by O. N :> and electrons 


At twilight the situation changes rapidly. The photoionization 
peak rises rapidly in altitude and the production rate falls sharply. 
Figure 4 shows the production rates for 48°S at 22:23 LT. At this 
time the solar zenith angle is 107°. so that the high altitudes are 
illuminated. There is a second peak at about 300 km due to trans- 
ported local photoelectrons which are produced during the ionization 
process which created the upper peak. Radiative ioss dominates most 
of the protile. Ai low latitudes at midnight, the peak production rate 
tends to zero at very high altitudes. When the production rate effec- 
tively drops to zero these values are omitted m Figure 2. so as not to 
obscure the main results. 



Fig. 4 Sources (a) and sinks (b) of 0~( : P) for the location 
corresponding to 48°S at 22:23 LT in Figure 2. 



Fig. 5. V ertical column integrated 
surface brightness (in Rayleighs) for 
the same case shown in Figure 2. 


^0 


Torr ei al . : Midlatitude Imerhemispheric Model of 7320 A Airgiow 


In the winter ( nonhem j hemisphere the twilights draw closer 
together but the photochemistry ot the daytime and twilight 7320 A 
emission is essentially the same as that described above. What ts 
interesting to note here is that there is some proo iction of OM 2 P) at 
hich latitudes even at midnight. The source is n.ipact ionization by 
photoelectrons from the sunlit conjugate hemisphere. The conjugate 
photoelectron source produces a column integrated surface bright- 
ness of 2 to 3 R of 7320 A emission. 

Figure 5 shows the global map of the 7320 A intensity for the same 
case as Figure 2. This is what would be seen by an instrument look- 
ing straight dow n at the Earth for a fixed LT. The midday intensities 
in the winter hemisphere are slightly larger than those in the summer 
hemisphere. This is because the winter hemisphere is somewhat 
colder, and radiative loss dominates to a lower altitude. 


Summary 

We have reported the first results of a midlatitude interhemispher- 
coupled model of the 7320 A airgiow. These results show the 
1 ;nal and seasonal variation for conditions representative of near 
solar minimum. The magnitude and extent of the conjugate photo- 
electron source is quantified, together with the sources and sinks at 
selected times. In another paper (Tonr et al.. 1989] we show further 
details ot the longitudinal and solar cycle variations. This model 
allows three-dimensional maps of any airgiow emission to be gener- 
ated. Any line of sight geometry can be projected through these maps 
and the volume emission rate integrated along the viewing direction 
lable comparisons with observation. A three-dimensional 
perspective of airgiow emission rates or brightnesses is of relevance 
to the interpretation of measurements taken from an orbiting vehicle. 
A satellite tlving in a circular orbit through the topology repre- 
sented in Figure 2 would observe large variations in emission 
through the structured regions. Semi-global maps of the type pro- 
duced here provide rapid insight into the sources of the variability. 
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ABSTRACT 


Emissions in the vacuum ultraviolet of the Lyman-Birge-Hopf ield 
(LBH) bands of N 2 have been observed at night from the S3-4 spacecraft 
and on Spacelab 1. Conway et al. have reported that the intensity of 
the LBH emission observed on S3-4 in the nadir varied approximately as 
C N 2 ] 3 with altitude, indicating a vehicle-atmosphere interaction. 
Observations were made at right angles to and into the velocity vector 
(ram) on Spacelab 1. The LBH intensity in the ram direction shows a 
factor of 3 to 4 decrease with respect to the 90° case which were taken 
two days earlier. The latter observations also suggest that the N 2 
Vegard Kaplan bands occur simultaneously with the LBH. The Shuttle 
LBH intensities are brighter than those observed on S3-4 by several 
of magnitude. The following model is proposed to account for 
the above observations. 

N 2 striking the vehicle surface is excited to vibrational levels 
near v = 13 , allowing the exothermic surface reaction N 2 (X) * + o -» NO + 

N to proceed rapidly, thereby producing N in abundance on the surface. 
Surface recombination of the atomic nitrogen leads to the formation of 
vibrationally excited and electronic states of N 2 including N 2 (A 3 l* ) 

V r 9 to 13 (7.6 to 8.2 eV) which we shall designate N 2 (A)*. After 
thermal accomodation on the surface, the N 2 (A) * desorbs into the local 
Shuttle environment. The long-lived metastable N 2 (A) * component is 
then collisionally excited to the a l H g state which radiates in the LBH 
bands. In the ram direction there is a large enhancement in gas 
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concentration which results in attenuation of inflowing ambient N 2 
and 0 and possible collisional deactivation of the N 2 (A) * reducing the 
LBH intensity. For collisional excitation of N 2 (A) * by ambient n 2 
the net altitude dependence is - [N 2 ] 3 above 200 km altitude, changing 
to - [N 2 ] 2 [0] at lower altitudes in agreement with this behavior seen 
m the S3-4 data. The S3-4 spacecraft intensities are smaller than 
the shuttle intensities because of the smaller spacecraft size, and 
because of increased attenuation at the lower altitudes, and possibly 
due to different reaction rates for different surface materials. 
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1 . INTRODUCTION 


Conway et al. (1987) have reported observations of vehicle 
induced emission in the Lyman-Birge-Hopf ield (LBH) bands of N 2 on the 
S3-4 spacecraft. The N 2 (a 1 ^) state, from which the LBH bands arise, 
requires 8.6 eV for excitation of the v = 0 level and up to 9.75 eV for 
excitation of the higher lying levels below the dissociation 
threshold. There is no identified source for these bands in the non- 
auroral nightglow. Figure 1 reproduces the results of Conway et al . 

t-h® altitude variation of the integrated band intensities between 
1400 and 1700 A for nadir viewing. They report that the intensity 
var i- es with altitude as [N 2 ] 3 (or [N 2 ] 2 [0]) . They demonstrated that 
fhe source mechanism is a three— stage excitation process involving 
the surface, then the expected LBH intensity would be *16 orders of 
magnitude short of explaining the observed intensities. 

Observations of vacuum ultraviolet (VUV) emissions were also 
made by Torr et al. (1985) on the Shuttle with the Imaging 
Spectrometric Observatory (ISO) . Observations were made in the 
nadir and at 90° to the velocity vector at night, the ram at twilight, 
and the wake in daylight. For the nadir and 90® observations, the 
Shuttle was flying in an upside down airplane mode, i.e. with the nose 
(x axis) pointing into the velocity vector and the bay (-z axis) 
pointing in the nadir. Under these conditions the ISO views out along 
the -z axis when the scan mirrors are stowed out of the field of view. 
(The open mirrors are stowed > 180® to the closed position) . The 90® 
to ram data were made looking under the — y wing using the scan mirrors 
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to point the f ield-of-view of the instrument between 90 and 250 km 
tangent ray height. The ram data were taken with the scan mirrors 
stowed and the line-of-sight of the instrument pointing directly into 
ram, i.e. with the -z axis of the Shuttle into ram. 

Intensities were observed to be brightest at 90° to ram, 
decreasing in the ram and wake directions. 

Figure 2 shows an example from Torr et al (1985) of the vacuum 
ultraviolet spectrum observed on Spacelab 1 at 250 km at night on 
December 5, 1983 at 90° to ram. Curve a shows the synthesized LBH band 
system for comparison. The maximum allowable intensities are 
constrained by the lower values observed at the shorter wavelength 
side. Thus, it is apparent that the VUV glow observed on Spacelab 1 
cannot be attributed entirely to the LBH bands. Curve b shows the 
effect of including the Vegard Kaplan bands of N 2 , where the maximum 
allowable intensities are constrained by the lower values at both the 
short and long wavelength sides of the measured spectrum. 

Torr et al. (1985) also identified the presence of strong NO 7 
bands at night at wavelengths longward of those shown in Figure 2. 
Figure 3 shows the effect of including the NO f and 3 bands which have 
transition characteristics which allow the system intensities to peak 
near 1650 A. Thus it is evident that a plausible synthetic spectral 
fit can be achieved with a composite system of bands comprising N 2 LBH, 
VK and NO e and 5 bands . 

It is clear that the spectral identification of the shuttle VUV 
emissions may not be unique, because there is no one-to-one match with 
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Figure 2 : An example of the vacuum ultraviolet spectrum observed on 

Spacelab 1 at 2 50 km on December 5, 1983 at ' 130° W, at 21 hours local 
time, solar zenith angle = 107 at mid-latitudes. A mirror was used 
to view the 90° direction across the -Y wing of the Shuttle. Curves a 
and b are synthetic spectra of the N 2 Lyman-Birge-Hopf ield and Vegard 
Kaplan band systems respectively. Curve c shows a composite spectrum 
of these two systems. 
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Fiqure 3 : A comparison of the observations shown in Figure 2 with a 

composite synthetic spectrum comprising the N 2 LBH, VK and NO e an 


bands . 



discrete spectral features, but only with the general envelope of the 
observed VUV pseudo continuum glow. Thus the identification of the 
VK emission is not definitive. Detailed fitting may not be possible 
because it requires several grating steps to cover the full VUV range, 
and temporal variations which are severe may distort the relative 
emission ratios . Given the [N 2 ] 3 dependence of the glow indicated by 
the S3-4 data it is evident that relatively small changes in the N 2 
concentration due to gravity waves for example, could cause 
significant changes in the VUV emission because of the strong non- 
linear dependence. 

If the presence of the VK bands is accepted, one is led to the 
conclusion that there could be a significant far-field component to 
the glow. This conclusion is based on the fact that the ISO field of 
view was directed away from shuttle surfaces, and for the viewing 
geometry in question the instrument itself was shielded from direct 
interactions with the ram flux by the Spacelab 1 module and the aft 
bulkhead. Slanger (1986) has presented evidence which strongly 
suggests that shadowing of surfaces from exposure to ram 
significantly reduces visible surface glow brightness. 

For comparison, Figure 4 shows the VUV emission observed in the 
ram direction at 250 km on December 7, 1983. The most striking 
feature of these data in comparison with that shown in Figures 2 and 3 
is a reduction in intensity of the N 2 and NO bands by of a factor of 3 to 
4. 

Given the fact that the inferred intensities for the VK system 
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(*3kR) can be used to roughly obtain the column abundance of N z (A) , 

which we find to be large ('6 x 10 9 cm" 2 ), we investigate the 

plausibility of collisional excitation of desorbed N 2 (A)* by the 
influx of ambient 0 and N 2 as a source of the LBH emission. We show 
that the apparent decrease in N 2 LBH (and VK) emission in the ram 
direction, may be due to the attenuation of the incident ambient flux 
of N 2 and 0 and quenching of N 2 (A)* by O and N 2 re-emitted by the 

surface. The concentration of 0 and N 2 is believed to be 

significantly enhanced on the ram side of the Shuttle (Rantanen and 
Gordon (1987) . We consider only excitation of N 2 emissions in this 
paper, and will deal with the source of NO excitation in a future 

paper. 

2 . THE MODEL 

We investigate the production of LBH emission via the reactions 
0 + N 2 (A 3 E + u ) shut -» N 2 (a 1 TI g ) + o (1) 

N 2 + N 2 (A 3 I + U ) shut -> N 2 (a 1 n g ) + N 2 (2) 

followed by: 

N 2 (a 1 n g ) -» N 2 (X 1 E + g ) + h k ( LBH ) ( 3a > 

N 2 (A 3 £ + U ) shu t - N 2 ( X 1 S + g ) + h v (VK) (3b) 

where N 2 ( A 3 l + u ) s h u t refers to metastable N 2 desorbed from the vehicle 
surface and the usual notation (0,N 2 ) refers to the ambient gases. 
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3 . THE PROPOSED LBH GLOW MECHANISM 

The desorbed vibrationally excited N 2 (A)* internal energy peaks 
between 7.6 and 8.2 eV. Both S3-4 and Spacelab 1 observations 
indicate LBH emissions arising from v' = 0 to 6, that is from 
vibrational levels lying between about 8.6 and 9.7 eV. The mean 
relative energy in the center-of-mass system available for processes 
(1) and (2) is -3.9 ± 1 and *4.6 ± 1 eV respectively. Thus, both 
processes are energetically capable of exciting N 2 states above the 
a^g state. 

The volume emission rate for the LBH system, ■<) , at a distance s 
from the surface is given by 

-n (s) = eF; (s) [N 2 (A)*] s (4) 

where Fj (s) represents the total incident flux of ambient neutral 
species j, <r jE represents an excitation cross-section for collisions 
with the j th particle species where j = 0 and 1 designates 0 and N 2 
respectively, and [N 2 (A)*] S is the concentration of desorbed 
vibrationally excited N 2 (A)* at a distance s-s 0 from the 
vehicle/ instrument surface at s 0 . The total LBH irradiance is given 

by 


7 


CO 


00 


{F(N 2 ) s <r 1E + F (O) s o e ) t^2 (A) ] S ds 


(5) 


LBH 


i)(s)ds = 


where <r 1E and (r 0E are cross-sections for collisional excitation of 
N 2 (A)* to N 2 (a) by N 2 and 0 respectively. The flux F(x) s is given by 


— T a x “ T aj 

F(x) s = F (x) ^e or F j = F joj e (6) 

where F refers to the unattenuated flux, x refers to O ( j =0) or N 2 

00 

(j = 1) , and r a j represents the attenuation depth at a point along a 
given line-of -sight s defined by 
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a j 


oo 


f 


in A (s)ds' 

/ 1 

s 


(?) 


where s' is in a direction antiparallel to the velocity vector, c j ^ is 
the collision cross-section for species j and i and n j (s) is the 
induced gas concentration of the i t ^ constituent in the environment of 
the vehicle at distance s from the surface. Equation (6) is 
formulated under the assumption that after one collision , incident 0 
and N 2 will no longer be capable of exciting N 2 (A) * to the a 1 n g state. 
To evaluate the viability of the proposed excitation mechanism we 
first consider the case where the attenuation of the incident flux by 
the gases comprising the vehicle environment is negligible, i.e., 
where << 1. Under these conditions, (5) simplifies to 

^ j 
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00 


[N 2 (A)*]ds 


( 8 ) 


LBH 


= l Fi 


j oo“ j E 
S 


= ^F jo3 <r jE [N 2 (A)*] C0l ( 9 ) 

where [N 2 (A)*] col is the column concentration of vibrationally 
excited N 2 (A) , and subscript j® refers to the unattenuated flux of the 
j th species. 

We use the Spacelab 1 measurements of the VUV VK band systems to 
evaluate [N 2 (A)*] col . Given that the LBH and VK nadir intensities 
are "3.5 kR and 3 kR respectively and using the relationship 


where I is the total intensity of a band system, Aj is the Einstein 
coefficient for all transitions from the i th vibrational level of the 
emitting electronic state X, and [Xj ] is the concentration of the i* h 
vibrational level, we evaluate the total column abundance of N 2 (A) 
from the intensities for the VK system, assuming an average lifetime 
of 2 seconds based on the results of Shemansky et al. (1971) for low- 
lying vibrational levels. Thus, we find that 

[N 2 (A)*] co1 r 6 x 10 9 cm" 2 (11) 

We use (9) and (11) to estimate the collision excitation cross- 
section needed to account for the observed LBH emission. 
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F„[N 2 (A)] col 

Where is an intensity weighted cross-section, and F^ is the total 
ambient neutral flux. Using the MSIS model atmosphere (Hedin, 1987) 
for the Spacelab 1 conditions, i.e. November/December 1983 at 250 km, 

IF: = 1.3 x 10 15 cm" 2 s _1 . 

j J 05 

Thus, 


5 x 10 4 * * * * 9 

1.3 x 10 15 x 6 x 10 9 


= 4.5 x 10" 16 cm 2 


which does not appear to be unreasonable. 

The obvious questions that arise are the following: 1) What is 
the source of the N 2 (A) molecules? 2) Why are the LBH and VK 
intensities weaker in the ram direction? 3) How does the [N 2 ] 3 (or 
[N 2 ] 2 [0]) altitude dependence arise? 


4. SOURCE OF THE N 2 (A) GI/DW 

In principle, there are several possible ways of generating 

electronically excited N 2 (A) . However, all but the surface 

recombination of atomic nitrogen can be effectively eliminated by 

energy arguments (Kofsky, 1988). Surface catalyzed inverse 

predissociation of N 2 appears to be the only process capable of 
accounting for the highly efficient excitation of N 2 (A) in high 
vibrational levels. Reduction of available energy in the center of 
mass system of the reactants, and the fact that a significant fraction 
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of the kinetic energy must be shared with the surface eliminates the 
direct excitation of N 2 investigated by Conway et al. (1987) or impact 
dissociation of N 2 (suggested by Green, 1984) as sources of 
electronically excited N 2 or of atomic nitrogen. The available flux 
of ambient thermospheric nitrogen would require unrealistic 
efficiencies for N-N recombination. 

Meyerott and Swenson (1990) have noted the importance of being 
able to account for the observed anomalous vibrational distribution 
of the LBH emission. They point out that direct recombination of N 
into the a 1 ^ state via the A 5 £ + g surface should result in predominant 
population of vibrational levels near v = 5 where the crossing of these 
states occur. 

In order to provide a comparison with our model we briefly 
describe their approach. To explain the far-field vehicle glow they 
propose resonance fluorescence scattering of EUV photons which will 
give rise to a 1 Ilg by cascade from higher lying singlet Rydberg states 
of N 2 . To initially populate these singlet states they invoke a 
Rideal surface recombination mechanism in which one N is attached to 
the surface and the second is ambient atomic nitrogen. The relative 
kinetic energy of the latter provides the energy source needed to 
access the singlet Rydberg states. It is presumed that the 
transition probabilities from these states will be consistent with 
dominant cascade to the v = 0 and 1 levels of the a 1 Ilg state. 
Recombination can occur via either one or both of the A 5 £ + g or 7 E + g 
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surfaces . 


Swenson and Meyerott (1988) have also proposed the gas phase 
reaction 

N 2 + 0 -» NO + N (13) 

as the main source of surface N, and their calculations suggest that 
(13) could provide a sufficient source of N to account for the LBH glow 
observed on Spacelab 1 and the S3-4 satellite. Because the branching 
ratio from the singlet states to the ground state of N 2 is about 30 to 
100 times larger than to the a 1 ^, strong EUV emission is predicted 
together with the Herman Gaydon near UV bands associated with the 

cascade process. 

It should also be noted that the 5 E + g and 7 S + g paths require one 
more step than the A 3 E + u in the recombination process . In the case of 
direct recombination into N 2 (A), one would expect high lying 
vibrational levels to be populated, with possible partial surface 
induced vibrational relaxation. According to Green (1984) some 
fraction of the molecules will leave the surface in the N 2 (A) state. 
Currently it is believed that recombination should be preferred 
through the 3 E + u channel, because this state correlates directly with 
the ground state atoms, and it is thought that laboratory observations 
of emissions from singlet states are due to quenching of high lying 
vibrational levels of the A 3 £\ state (B. D. Green private 
communication, 1989) . 
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As mentioned above, the N 2 (A) vibrational distribution observed 
on Spacelab 1 is constrained rather well by the envelope of the pseudo 
VUV continuum. The distribution required to fit the data shown in 
Figure 2 is given in Figure 5. These values show that the N 2 (A) must 
be populated predominantly in v = 7 to 13. Also measurements of N 2 (A) 
emissions from v' = 0 in the near UV preclude low lying vibrational 
levels as a possible precursor. Although collisional excitation of 
N 2 (A) (v" » 0) to a : H g (v' » 0) is preferred from the Franck-Condon 

perspective, transitions which violate the Franck-Condon principal 
are not precluded for heavy particle collisions. Thus collisional 
excitation from high lying vibrational levels of N 2 (A) could populate 

low lying vibrational levels of N 2 (a). 

As stated earlier, the only mechanism capable of producing 
N 2 (A)* is surface catalyzed recombination of N. Thus, the next step 
is to identify a source of N that will yield the correct altitude 
dependence of the LBH emission. We consider the surface reaction of 0 
and N 2 * where the following sources of O and N 2 * are proposed. 

0 ( * 5eV) + surface — % O surf ( 14 ) 

N 2 (*9 eV) + surface » N 2 (X)„2i3 + - 5 " 5 eV ( 15 ) 

followed by the exothermic surface reactions 

N 2 (X ) V>13 + °surf NO surf + N surf (16) 
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Figure 5: The N 2 (A 3 E U + ) vibrational distribution required to 

produce the synthetic spectral fit shown in Figure 2 . The results are 
nomalized to unity at the peak value. 



N surf + N surf * N 2 (A) surf ( 17 ) 

followed by desorption of the N 2 (A) * and Reactions (1) and (2) . The 
species is probably long lived with a radiative lifetime of the order 

of a second. 

In this model, the N 2 (A) * could be formed on the surfaces of the 
shuttle wings. Then the desorbed N 2 (A) * could flow through the field 
of view of the instrument when viewing over the wing in a direction 90 
to ram. 

Alternatively, it is possible that direct impact of surface O and 
N 2 by the flux of ambient O and N 2 could also result in the production 
of N and NO, however, this will result in a different scale height for 

the observed LBH and VK emissions. 

Processes (15) and (16) need to be highly efficient in order to 
result in significant production of N. The ram flux of 0 will be the 
major source of surface 0 needed in (16) . As pointed out by Kofsky 
(1988) , the probability for N recombination ranges between 1 and 10 6 
depending on surface material (Evenson and Burch, 1966; Halpern and 

Rosner, 1978, 1982). 

The production of N, Q (N) is given by 

Q(N) = kj [N 2 (X) * surf ] [0 surf ] (18) 

It will become evident that the processes which control the 
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concentration of N 2 (X)* surf and 0 surf will determine the dependence 
of the LBH glow on altitude. In Section 5 we demonstrate that an [N 2 ] 3 
dependence requires that [N 2 (X)* surf ] scale linearly with ambient 
[N 2 ] and that [0 surf ] scale as [0] 1/2 . Here we show that these 
dependences can be provided in terms of likely processes. If we 
consider the steady state production of N 2 (X) surf by ambient inflow 
and loss by (16) and desorption, then we can write: 

e 1 F(N 2 ) = (k 1 [0 surf ] + Ji) [N 2 (X) * slir f ] ( 19 ) 
where J 1 represents a desorption coefficient (units: s l ) 


If >> t^sur f 1 f 

(which is possible because of the internal excitation energy of N 2 (X) 
which may be available to overcome the surface bonding energy) 


[N 2 (X) S u r f] 


e ! F (N 2 ) 


( 20 ) 


which satisfies the reguired dependence. This condition implies 
that Reaction (16) is a major source of N but a minor loss for N 2 (X)*. 

If we assume that the primary loss mechanism for surface 0 is 
recombination (at least for Shuttle surfaces), either directly or 
catalytically such that the net effective process can be represented 

by 


^surf + ^surf 


0 


2 s u r f 


( 21 ) 
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and if the source of 0 s u r f , 


Q (0) , is the flux of ambient 0 , F ( 0) , then , 


The assumption that surface recombination of 0 is a major loss 
process for 0 is supported by the well known fact that satellite borne 
closed source mass spectrometer measurements of 0 yield both 0 and 0 2 
signatures when 0 is the dominant ambient species present. Thus the 
atomic nitrogen production rate (18) can be written as 


<iF(N 2 ) 

Q(N) = k x 

Ji 


(e 0 F(O)) uz 


-3 j 


a [N 2 ][0] 1/2 


(23) 


where 

-(?! + t 0 /2) 

a = *1^7 k 1 v 3/2 /(J 1 -Jk7)e 


(24) 


where r j 


00 


/ 

S 


1 < j j in jds' 


(25) 


from (7), and s' is always anti-parallel to the velocity vector. 
Subscript j refers to the incident particle and i to the attenuating 
constituents. 

If Reaction (16) is sufficiently exothermic desorption of N may 
be the dominant loss process and assuming state steady state 
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conditions apply, 

Q (N) = Jj [N] 

Then 


Q[N] 
[N] = 

J 2 


( 26 ) 


and the production rate of desorbed N 2 (A) from (17) is given by 


F (NjA) * 


e 2 k 2 


r Q (N) ' 

(, J 2 * 


-(2T! + T 0 ) 

7[N 2 ] 2 [0]e 


where, using (24) 
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e 2 k 2 




f O f l 2f 2 2 ^2 V 

j i 2j 2 Z k 3 


(27) 


(28) 


e 2 is the yield of desorbed N 2 (A)*. 

Once the N 2 (A)* desorbs from the surface it may be quenched by 

constituents of the gas cloud. If quenching of N 2 (A) * is negligible 
then the concentration of N 2 (A)* close to the surface is given by: 


- ( 2 T J + T o ) 

[N 2 (A)*] a [N 2 ] 2 [0]e 


(29) 


It is known from laboratory and aeronomical data that N 2 (A) v > 0 is 
quenched by 0, ando, (see Torr and Torr, 1982; Piper etal. 1981: sharp 
and Torr, 1979) . Quenching of N 2 (A) * by HO is also rapid (6 x 10‘ “ 
cm 3 s" 1 ) . Thus, it is reasonable to assume that N 2 (A) V27 will be 
quenched both vibrationally and electronically. It is likely that 
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all major constituents of the vehicle gas cloud could quench the 
N 2 (A)*. The fraction of the N 2 (A)* flux that is quenched will be 
determined by the relative magnitudes of the transport and chemical 
(or quenching) lifetimes. 

Having computed the desorbed flux of N 2 (A) it remains only to 
include both the effects of transport and quenching on the 
distribution. To evaluate the effects of both transport and 
quenching on N 2 (A) * we solve the continuity equation (and dropping the 
superscript *) 

3(N 2 A) 1 

= q - B[N 2 A] v(0G) (30) 

at g 

where q = gas phase production rate of N 2 (A) = 0 

6 = N 2 (A) loss frequency = + A v k 

kj = quenching rate coefficient for species i = ^iV; 

A v k = Einstein coefficient for the VK system 

n j = concentration of the i th constituent of the vehicle 
gas "cloud" 

^ = v t [N 2 A] = N 2 (A) flux at point s 

v, = is the bulk velocity of desorbed species which is 
characterized by the temperature of the surface at 
which desorption occurs 

G = geometrical function which characterizes the effects 
of the radial outflow of gas 

t = time 

For purposes of deriving analytical solutions to (30) to 
facilitate discussion, we assume steady state conditions, 
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(31a) 


1 

6[N 2 A] = (GV0 + *VG) 

G 

or 

7[N 2 A] vG 

“6/v t = + — 

[N 2 A] G 

If 7 = d/ds for a given viewing direction and 


& = 


+ A 


VK 


- V T ? ff q i n i 


+ A 


VK - 


<r q nv T 


l VK 


where 7 q is a density weighted quenching cross-section 
n = Sn i 

Upon integrating (31) and adding attenuation of the 
flux to 6 we obtain 


[N 2 A] 


~(T q + 2T X + T 0 ) 

[N 2 A] so G 0 e 


G s 


where r q 



+ A VK /v T )ds 


and Tl and r 0 are the attenuation depths for the incident 
respectively . 


(31b) 


(32) 


incident 


(33) 


(34) 


N 2 and 0 
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5 . 


THE ANGULAR DEPENDENCE OF THE SHUTTLE VUV GLOW 
Using (34) in (4) the LBH volume emission rate is given by 


L B H ( S ) 


-T j (S) 

?<7j E^j so G O e 

G S 


(35) 


which can be written as 


Ti (s) = 


Effj ^ds' 


+ 


<r q nds 


+ 2 


Effjjnjds 7 + 


? c o i n i 


(36) 


s o 


S 0 


S 0 


-Tj (s) 

? , lj(®o) G 0® (■ 

^ L B H = 5 ” 

G s 

The LBH integrated emission rate can be readily evaluated 


B H 


CD 


(-')<■) 

e ds 



' G s 

s o 


(38) 


where, as mentioned previously, s is along the line of sight, and s' is 
antiparallel to ram. The attenuation terms represent the following 
in order of appearance: 

1) Attenuation of the incident ambient flux of 0 or N 2 (j 
0 or 1 respectively) between ® and point s 
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2) Quenching of N 2 (A)* between the surface and s 

3) Attenuation of the ambient N 2 flux to the surface which 

produces N 2 (vib.) 

4) Attenuation of the ambient 0 flux to the surface which 
produces surface 0. 

in practice, (38) is evaluated numerically for any given viewing 
direction. However, since it is relatively straight forward to 
evaluate the case analytically for the line of sight directed into 
ram, and since this case is the most relevant for our discussion, we 
will utilize this case as a point of departure for further discussion. 

First we define 

(39) 

<r j n Scj^n^ 

For ram viewing (36) then becomes: 



f 

f 

r j (ram) = <rj 

J 

nds + ff q 

) 

nds + ( 2 <t j + <r 0 ) 

t 


nds 


(40) 


S 0 


S 0 


i. e. 


Tj (ram) = r a j(s) + T a (s) + 2 t x + t 0 


(41) 


where 


nds = N, the total ram column density of the 
vehicle gas cloud 


(42) 


S 0 


Thus (38) becomes 
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I LBH ( raltl ) = ? T| oj G o 


f -Tj (ram) 


ds 


S 0 


where -o 0 j = T i 0 j( s o) " F j co 0 ? e [ n 2 A 1 s ° 
If G s = s 2 and <rj r c q 


(43) 

(44) 

(45) 


I LBH( raltl ) - ? 1, 0j S 0 e 


— t j (ram) 


(46) 


where s 0 is 


of the order of the vehicle dimensions (L) , and 


T = T + t 0 + 2Ti where r 2 = T aj + T q - <r aj 


nds 


(47) 


S 0 


The attenuation depths for the incident 0 and N 2 fluxes depend 
critically on the surfaces underlying the line-of-sight of the 
instrument. For example, if the line-of-sight runs parallel and 
close to a large surface such as the Shuttle wing, the flux of desorbed 
N 2 (A) * will be largest if the wing surface is normal to ram. However, 
this will also result in large values for all the attenuation depths , 
yielding a large value for Tj . On the other hand, if the wing edge is 
directed into ram, as in the airplane mode, the flux of desorbed N 2 (A) 
and the attenuation factors 7 aj , t 0 and ?! will be reduced roughly in 
proportion to the area exposed directly to ram. The attenuation of 
incident ambients will decrease as exp [-(? aj + T 0 + 2r x )], that is 
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exponentially with respect to the decrease in the flux of N 2 (A) . 

Similarly, when the line of sight is directed into ram, there is a 
large buildup of quenching/ attenuating constituents which may cause 
the attenuation and quenching terms to overwhelm the enhancement in 
N 2 (A)*. 

In what follows we report the results of a preliminary assessment 
of the magnitude of these effects for the databases discussed in the 
introduction. 

From (38) we see that the angular dependence of the emission will 
depend on the attenuation depth Tj (9) which varies as the column 
concentration N of desorbed species (see (40) for example) which can 
be obtained via detailed modeling of the vehicle "gas cloud" dynamics, 
which is currently underway and will be reported at a later time. 

In order to estimate values of Tj (0) for 0 = 0 and 90° we 
use the calculations of Rantanen and Gordon (1987) . Their results 
typically indicate a decrease in column density from ram to 90° (for 
the vehicles studied) of about an order of magnitude. 

If a o = 7.i - *o = i 2.5 x 10- 15 cm 2 for example, the 

N 2 (A) quenching rate coefficient by N 2 is given by k N 2 r v T <r R : 4 x 10 x 
2. 5 x 10“ 1 5 = 1 x 10' 1 0 cm 3 s" 1 which may not be unreasonable for high 

lying vibrational levels. 

Thus t- , . - 10' 14 N = 1.2 and the column density of the gas 

cloud needed to account for the attenuation of the ram flux is 1.2 x 
10 14 cm' 2 . If L: 30m, (see(46)) the near surface gas density is 4 x 
10 10 cm" 3 which is reasonable and corresponds to a ram enhancement of a 
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factor of 27. 


(Note the Rantanen and Gordon (1987) ISEM model predicts a non- 
linear increase in gas density with a few meters scale length near the 
surface, which can severely affect the attenuation of molecules en 
route to the surface.) 

Table 1 provides a summary list of approximate numerical values 
reguired for relevant parameters if the proposed model is to be 
viable . 


6. ALTITUDE DEPENDENCE OF THE S3-4 GI/DW 
ANALYTICAL FORMULATION 

Differentiating (46) with respect to altitude h, for a single 
species and simplifying notation, we obtain 


dl dl n -r 


-r dr 


dh dh 


e - 1 0 e — 

dh 


(48) 


where I 0 = ri 0 L, the unattenuated LBH intensity from which we deduce the 
height H for the LBH emissions, namely, 

-1 1 dl 1 dl 0 dr 

H I dh I 0 dh dh 


(49) 


i.e. 


Ill 
H H 0 H t 


(50) 


where H 0 is the scale height associated with i» 0 defined by (51) 
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TABLE 1 


SPACE LAB 1 PARAMETERS 
90° DATA 


F(N 2 ) 
F (O) 


= 4.5 x 10 1 4 
= 8.5 X 10 14 


cm" 2 s" 1 

) MSIS 1986 


AT 

cm" 2 s" 1 

J 250 km 


[N 2 A]* co1 = 6 x 10 9 cm" 2 (measured) 

F(N 2 A) = 8.7 x 10 10 cm" 2 s" 1 

N 2 (A)/N 2 = 1.9 x 10“ 4 

N 2 (X) V > 13/N 2 = 0.08 ASSUMED 

N /N 2 * =0.1 ASSUMED 

N 2 (A)/N RECOMB. =0.02 

N 2 (a) /N 2 (A) =0.04 

LBH PHOTONS/N 2 = 7.6 X 10" 6 

VK PHOTONS/N 2 = 6.6 X 10" 6 



From (35) and (29) 


ri 0 a [N 2 ] 2 [0] n j (n j = [O] or [N 2 ] 

Thus 


(51) 


1 



and 


3 1 

+ 

H N2 H 0X 


12 2 
Ho 2 H ox 


for nj = [N 2 ] 


for nj = [O] 


(52a) 


(52b) 


where H n2 and H ox are the scale heights of atmospheric N 2 and 0 
respectively . We now evaluate H T . 

From (40) and (41) we approximate Tj by 


T 


and 



dT <rLdn 0 
dh dh 


<rn 0 L 


where 

n 0 = fn r T r e 


T 


h 

( dh 



(53) 


(54) 


(55) 


( 56 ) 
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TABLE 2 


Height Dependence for N 2 Collisions (Reaction 2) 


LBH 


Ht 

(km) 

[ 0 3 3 

(cm 3 ) 

[N 2 ] 

(cm -3 ) 

a 0 

n 

l-a 0 

a l 

m 

3-a^ 

ALT DEPENDENCE 
[N 2 ] m [0] n 

180 

5 . 50x10 9 

5 . 60x10 9 

1.10 

-0.10 

1.12 

1.9 

[N 2 

] 1 ‘ 9 [ 0 ] ” 0 ' 1 

190 

4 . 34x10 9 

3 . 70xl0 9 
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[Nj 
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0.68 

0.32 

0.50 
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[N 2 

] 2 . 5 [O] 0 • 3 
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0.56 

0.24 

2.8 

[N 2 

K> 

00 
r 1 
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240 

1. 50X10 9 

6.00x10 s 

0.30 
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2.9 

[N 2 

] 2 • 9 [0] 0 • 7 
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Ht [N 2 ] m [Q] n 

180 [N 2 ]°- 9 [0]°- 9 

190 [N 2 ] 1 ‘ 3 [O] 1 ' 1 

200 [N 2 ] 1 * 5 [O] 1 ‘ 3 

220 [N 2 ] 1 ’ 8 [O] 1 ' 6 

[N 2 ] 1 * 9 [O] 1 * 7 


240 



n r = the total ambient density at reference height h r 

H n = the mean density weighted scale height of the ambient 
gas 

T = atmospheric temperature 
f = ram density enhancement factor at s 0 

We assume that the atmosphere is isothermal (i.e. T — T r ) and 
obtain 


1 dT -fffLn 0 dz r 
H t dh dh 

= -fffLn 0 


(57) 


(58) 


From Rantanen and Gordon (1987) for L = 6m (assumed for the S3 4 
spacecraft), f z 30. 

Recall 7 = E<Ti ~ 10" 14 cm 2 ( 59 ) 

i 

Thus f <rL z 2 x 10" 10 cm 3 for the S3-4 spacecraft. 


Table 2 gives average concentrations for O and N 2 from MSIS 86 as 
a function of altitude for the S3-4 spacecraft. 

Equation (57) can be written in terms of the individual number 

densities of 0 and N 2 , i.e. 


= -f crL 


H. 


[O] [N 2 ]) 

+ 

Ih 


OX 


H 


(60) 


n2' 
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= - a o " a l 

H 0X H M2 

Using (61) with (52) in (50) 

1 (3 - a t ) + (1 - a 0 ) 

H H K2 H ox 

where a 0 = f<rL[0] 
a x = fffL[N 2 ] 

Equation (62) allows the altitude dependence of the LBH glow to 
be expressed in the form 

i L B H [N 2 ] to] 

where 

m = 3 ol x and n 1 o 

Table 2 also lists values for a 0 , a lf m and n as a function of 
height for the S3-4 spacecraft. Figure 6 shows a comparison of the 
intensity of the computed LBH emission compared with the S3-4 
observations. Here the intensities have been scaled to the 
measurements at *200 km, in order to evaluate the altitude dependence 
predicted by the model, which is reasonably good. If the quantity f<rL 
has been overestimated, i.e. a 0 and a a are too large, then the 
calculated scale height will approach the high altitude limit of 
[N 2 ] 3 [0] more rapidly than observed. If on the other hand <x 0 and a i 
have been underestimated, then there would be indications of a more 


( 61 ) 


(62) 


27 



SPACECRAFT ALTITUDE (km) 




rapid increase in scale height below '200 km, which would ultimately 
lead to the formation of a peak in glow brightness, because of the 

effects of increased attenuation. 

It should be noted that these calculations do not take into 
account possible surface saturation effects. This will be done 
during the next phase of the work when the emission will be 
comprehensively numerically modeled. 

7 . SURFACE GENERATED LBH GLOW 
Equation (17), namely 

XT * ( 63 ) 

N surf + N surf N 2 surf 

also allows for the direct production of N 2 (a 1 ^) on the surface. The 
radiative lifetime for the a x n g state is -100 us which yields an e 
folding distance of 4.4 cm, which would not be observable by the ISO, 
but could be the mechanism responsible for the S3-4 glow. The scale 
height dependence of the glow is then given by (27) as modified for 

N 2 (a) , namely 

-(2t 1 +t 0 ) 

F(N 2 a) S0 = e 3 k 2 (Q<N)/ J 2 > 2 = ^a[ N 2] 2 t°] e (64) 


where 


^ a 
f 3 

and 


7f 3 (65) 

e 2 

efficiency for the production of N 2 (a 1 II g ) by (17) 
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-(2t 1 +t 0 ) 

^ L B H = ^ a ] [°] e (66) 

V T 

Thus, the LBH glow will vary as [N 2 ] 2 [0] at heights when 
attenuation of the incident 0 and N 2 ambient flux is negligible, 
changing to a larger scale height when attenuation becomes 
significant at lower altitudes, which is also qualitatively 
consistent with the S3-4 observations shown in Figure 1. It is 
possible that surface materials on S3-4 result in negligible 
production of N 2 (A) and with a small production of N 2 (a) accounting 
for the lower intensities than those observed on Spacelab 1. 
CONCLUSIONS 

Observations of vacuum ultraviolet emissions (VUV) made by the 
Imaging Spectrometric Observatory (ISO) on Spacelab 1 in 1983 and the 
S3— 4 spacecraft are interpreted in the light of the following proposed 
model. Key features of the observations include a decrease in VUV 
emissions in the ram direction on Spacelab 1 compared to 90° to ram, 
and an [N 2 ] 3 altitude dependence on S3-4. 

It is argued that since the ISO optical surfaces were shielded 
from the direct ambient neutral flux, the observed emissions must be 
of far-field origin. The Spacelab 1 observations could be 
synthetically fitted with the following bands: N 2 LBH and Vegard 
Kaplan, NO £ and 5. Because the LBH radiative lifetime is short 
(~100us) and because of energy considerations, it appears that 
collisional excitation by the inflowing ambient N 2 or O of a long lived 
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excited state of N 2 to the a 1 !^ state is needed to account for the LBH 
component of the glow. The metastable N 2 (A 3 E + U ) species not only 
meets this need, but it also accounts for the presence of the inferred 
VK bands. Interpretation of the NO emissions is deferred for a later 
paper. Based on the measured LBH and VK intensities, the cross- 
section for excitation of N 2 (A) to N 2 (a) is estimated to be - 4 . 5 x 

10' 16 cm 2 . 

A mechanism is proposed for the surface production of N 2 (A) by 
inverse predissociation of surface N. To provide an adequate source 
of N which yields the correct altitude dependence observed on S3-4 it 
is suggested that vibrationally excited N 2 recombines on the surface 
with O yielding NO and N as a product. The vibrational excitation 
which is needed to render the reaction of N 2 with O exothomic is 
generated by the impact of N 2 with the surface. The N 2 (A) formed on 
the surface then desorbs into the cloud surrounding the vehicle. The 
N 2 (A) flux is estimated to be -lO^cm 2 s 1 for the Spacelab 1 

conditions . 

To explain the angular dependence of the glow observed on 
Spacelab 1, it is argued that attenuation of the ambient flux of N 2 and 
O though the gas cloud plays a significant role. First, the N 2 and 0 
that must reach the surface to produce surface N is attenuated, and the 
attenuation of each of these streams enters multiplicatively into the 
surface production of N 2 (A) . Second, the N 2 (A) in the gas cloud may 
be quenched. Third, the incident N 2 and O impacting the N 2 (A) in the 
gas cloud is also attenuated. The net result is an accumulation of 
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several attenuating terms which become significant in determining the 
angular dependence of the glow. When the Shuttle flies with large 
surface areas such as the wings facing into the ram direction, there is 
a large buildup in the concentration of desorbed gases around the 
vehicle, which could result in significant attenuation. When the 
Shuttle flies in an airplane mode, there is an exponential reduction 

in attenuation (-an order of magnitude) . 

On Spacelab 1, the ISO instrument viewed directly out of the bay 
when the -z axis was directed into ram, and over the -y wing when the 
shuttle flew in the airplane mode. The larger intensities observed 
in the latter case are consistent with the exponential decrease m 
attenuation resulting in larger intensities at 90 than in the ram 
direction. 

This model was also applied to an analysis of the S3-4 data which 
covered an altitude range from -180 to 240 km. It is demonstrated 
that at high altitudes, when attenuation is less significant, the 
altitude dependence of the LBH glow reduces to [N 2 ] [O] or N 2 [0] , 

and at low altitudes this changes to an [N 2 ] or [N 2 ] [0] scale height 
depending on whether excitation of the a 1 ^ state is by N 2 or 0 
collisions respectively. A simple model calculation, with 
theoretical intensities normalized to the S3-4 observations at 200km, 
yields reasonably good agreement for the variation of the LBH glow 

with altitude. 
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